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Foundations of Control 


tive control, many recent papers have been 
devoted to the philosophy of automatic con- 
trol. For example, in the Dallas Control Confer- 
ence there were papers entitled “Fundamental 
Theory of Automatic Linear Feedback Control 
Systems,” “Impact of Information Conversion on 
Control,” “General Approach to Control Theory 
Based on the Methods of Lyapunov,” and “On the 
General Theory of Control System.” In the 1960 
Moscow Congress Program, there are papers en- 
titled “Information Flow Criteria for Feedback 
Control Systems,” “Theory of Nonlinear Control,” 
“Adaptive Control System Philosophy,” and “The 
Philosophy of Control.” Many of these papers dis- 
cuss the advantages and disadvantages of feed- 
back, the importance of the sensitivity concept, 
the control of information flow or energy conver- 
sion, the relationship between reliability, efficiency, 
and control loop configurations, the fundamental 
structure of control loops, and the deeper thoughts 
pertaining to memory, adaption, and learning. 
What is the significance of the return to funda- 
mental explanations of a field that has been grow- 
ing almost exponentially in the last twenty years? 
Certainly, basic control principles have been un- 
derstood during this period, but at the same time, 
they have been occasionally misinterpreted or mis- 
used, and it is a healthy sign, a sign of maturity 


| Daeewerte instigated by the concept of adap- 


perhaps, to note that control engineers are pausing 
to reflect on the fundamental knowledge that gov- 
erns their sphere of activity. Some of the reflections 
are made verbally and others are expressed mathe- 
matically. Although fundamental in nature and 
sometimes inapplicable, thinking of this nature 
often leads to more useful concepts and broader 
applications. 

The advance of control theory has been based 
on a combination of mathematical theory and its 
practical interpretation in constructing physical 
equipment. Future advances will undoubtedly be 
derived from applications of different mathemati- 
cal innovations, and this issue of the TRANSAC- 
TIONS includes several papers which feature new 
applications of mathematics to the control field, 
Some of them may have no immediate application, 
but enterprising control engineers may find them 
useful in the near future. The excellent, thought- 
provoking paper by Bellman and Kalaba is a par- 
ticular example. It stresses the underlying princi- 
ples of control including feedback and adaptive 
processes, and it attempts to define them with 
mathematical terms as a foundation for future 
work in adaptive control. It is hoped that future 
TRANSACTIONS will feature some papers which 
utilize these mathematical tools in the develop- 


ment of new control systems. 
—The Editor. 
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The Issue in Brief 


As noted in the editorial, many of the papers in this issue are 
concerned with basic mathematical tools that may be applied to 
the field of automatic control. Also included is a discussion of the 
International Federation of Automatic Control, a second installment 
of a tutorial paper on vector quantities in control, and another bibli- 
ography on Time Lag Systems. 


The International Federation of Automatic Control, Harold Chestnut 


As explained in the previous issue,! the American Automatic 
Control Council (AACC) represents the United States in the Inter- 
national Federation of Automatic Control (IFAC), and the PGAC 
represents the IRE in the AACC. We are fortunate to have the first 
president of IFAC, Harold Chestnut, discuss the policies, progress, 
and significance of this international organization. 


Dynamic Programming and Adaptive Processes: Mathematical 


Foundation, R. Bellman and R. Kalaba 


This paper discusses fundamental aspects of control systems, the 
principles of feedback, and the nature of adaptive processes. The 
central purpose of the paper is to provide a foundation for the mathe- 
matical treatment of broad classes of the adaptive processes through 
the use of the concepts of dynamic programming. 


The Properties and Methods for Computation of Exponentially- 
Mapped-Past Statistical Variables, Joseph Otterman 


Various tools are used in statistical analyses and this paper dise 
cusses the properties and methods of computation with an ex- 
ponential weighting function decreasing into the past. This weighting 
function results in ease of computation, and it provides simple 
theoretical relationships when the interest is focused on the recent 
behavior of the process involved. Analog computer circuits and digital 
computer flow diagrams which serve to compute the exponentially- 
mapped-past statistical variables are presented. 


Generalized Weighting Function and Restricted Stability of a Linear 
Pulse-Modulated Error Feedback System, William A. Janos 


An extensive mathematical development is used to produce a 
generalized weighting function, for a linear feedback system with a 
pulse-modulated error signal, in the form of a matrix operator acting 
on an input vector consisting of the first R-1 input derivatives where 
R is the order of the unmodulated closed-loop system. Since the 
system operator has the form of a finite dimensional matrix, it 
becomes possible to make and to realize more stringent conditions 
on the transient stability in terms of a preassigned bounded output 
after a preassigned time. Two examples are given to illustrate the 
procedure. 


1 R, Oldenburger, “The American Automatic Control Council,” IRE TRANs. 
on AUTOMATIC CONTROL, vol. AC-4, pp. 5-7; May, 1959. 


Optimization Based on a Square-Error Criterion with an Arbitrary 
Weighting Function, G. J. Murphy and N. T. Bold 


In recent years many error criteria have been devised to measure 
system performance. This paper briefly reviews and compares most 
of these error criteria, then a new, generalized version of several 
methods is introduced, and it is referred to as the mean-weighted- 
square-error criterion. An optimum, realizable weighting function 
is derived by minimizing the error defined by the new criteria. 


Multiple-Rate Sampled Data System, Lester A. Gimpelson 


Large classes of digital systems employing computers have sev” 
eral pulse rates for economy or because of the nature of the com7 
putation. This paper discusses techniques for analyzing these sys” 
tems containing both discrete and continuous variables with flow 
graphs and submultiple rate Z transforms. Tables and examples are 
given to illustrate the techniques involved. 


Automatic Control of Three-Dimensional Vector Quantities—Part 
II, Allen S. Lange 


This is the second part of a three-part tutorial article on the 
use of vector algebra in control systems. The first part considered 
the position vector problems encountered in spherical trigonometry 
and this part introduces the angular velocity vector for the purpose 
of analyzing and designing geometric stabilization systems. 


Root Locus Properties and Sensi tivi ty Relations in Control Systems, 
Hanoch Ur 


Properties of the well-known root locus including pole sensitivity, 
angle of slope, and curvature are investigated in this paper, and a 
relationship between the sensitivity function and pole sensitivity is 
established. It is shown that the sensitivity determines variations in 
the transfer function due to large as well as small variations in the 
loop gain. Additional properties of the locus are also developed. 


Time Lag Systems—A Bibliography, N. H. Choksy 


Previous TRANSACTIONS have included several bibliographies: 
This particular one is a large supplement to the one given in the 
May, 1959, issue? and it includes a thorough discussion of the 
mathematical characterization of time lags. 


Correspondence 


An interesting property of the root locus, often overlooked, is 
clearly explained by C. S. Lorens and R. C. Titsworth in their dis- 
cussion on the “Properties of Root Locus Asymptotes.” 


PGAC News 


Details of the progress being made in organizing the first IFAC 
technical meeting in Moscow, as reported by Professor Letov of the 
U.S.S.R., are given to indicate the complexity and difficulties that 
attend an international meeting of this size. 


2 R,. Weiss, “Transportation lag—an annotated bibliography,” IRE TRANS, ON 
AUTOMATIC CONTROL, vol. AC-4, pp. 56-68; May, 1959. 
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The International Federation of Automatic Control” 
HAROLD CHESTNUT 


Summary—The initials IFAC, standing for the International 
Federation of Automatic Control, are appearing more and more in the 
control systems literature these days. It is appropriate that readers 
of these TRANSACTIONS be more fully aware of this organization 
in which they participate through the IRE’s membership in the 
American Automatic Control Council. What are IFAC’s purposes, 
who are its members, how does it operate, what has it accomplished? 
Now, after two years of existence for IFAC, these data are better 
known and have more significance. 


PURPOSE 
AN STATED in its Constitution, the purpose of 


IFAC is “to promote the science of automatic 

control.” Automatic control is deemed to cover 
the field of open and closed loop (feedback) control of 
physical systems in theoretical and applied aspects. 

“TFAC is to promote the science of automatic control 
by: 1) The interchange and circulation of information on 
automatic control activities in cooperation with national 
and other international organizations; 2) International 
Congresses; 3) Such other means as may be considered 
desirable as for instance publications.” 

Its purpose is to strengthen the national organizations 
and technical societies working in the field of automatic 
control by providing them with a way of working to- 
gether on an international basis. It also provides a 
means for individual members of these societies to learn 
of their counterparts’ activities in other countries and 
to establish contact and exchange technical information 
more directly. 


ORIGIN AND MEMBERS 


In a recent article,! Rufus Oldenburger has described 
the significant facts concerned with the founding of 
IFAC in Paris, France in September, 1957. IFAC is a 
federation of member organizations each representing 
the technical societies interested in automatic control for 
the country represented. “For each country, one scien- 
tific or professional engineering organization or one 
council formed by two or more such organizations hav- 
_ ing an interest in the field of automatic control shall be 
_ eligible for membership of IFAC.” 


On this basis, the following 22 nations are represented 


listed. 

-Austria—Osterreichischer Arbeitsausschuss fiir Automatisierung, 
e Osterreichisches Produktivitats-Zentrum, Hohenstaufengasse 3, 
— Wien 1. 


‘7 


a 


* Manuscript received by the PGAC, October 22, 1959. 
+ General Electric Co., Schenectady N.Y. _ ap ; 
1, Oldenburger, “The American Automatic Control Council,” 


RE TRANS. ON AUTOMATIC CONTROL, vol. AC-4, pp. 5-7; May, 1959. 
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in IFAC by the member organizations and addresses 


Belgium—lInstitut Belge de Régulation et d’Automatisme, 3 
Ravenstein, Bruxelles. 

China—Chinese Association of Automation, c/o Institute of Auto- 
mation, Academia Sinica, Peking. 

Czechoslovakia—Laborator pro Automatisaci a Telemechaniku, 
Cekoslovenska Akademie Ved, Zitna 28, Praha. 
Denmark—Dansk Ingeniorforening, 29, Vester 

Copenhagen, 

Finland—Industrial Instrumentation and Control Engineering Pro- 
fessional Group, c/o EKONO, E. Espl. 14, Helsinki. 

France—Association Frangaise de Régulation et d’Automatisme 
(AFRA), 19, Rue Blanche, Paris. 

Hungary—Hungarian National Committee of IFAC, c/o Budapesti 
Muszaki Egyetem, Villamosgepek Uzemtana Tanszek, Stoczek-u. 
2 T.-EP. 111. 36, Budapest. 

India—The Institution of Engineers, 8 Gokhale Road, Calcutta. 

Israel—Ministry of Defense, Scientific Department, Hakirya, P.O. 
Box 7057, Tel-Aviv. 

Italy—Consiglio Nazionale delle Ricerche, Ufficio relazioni inter- 
nazionali, Piazzale delle Scienze 7, Roma. 

Japan—National Committee of Automatic Control, Science Council 
of Japan, Ueno Park, Tokyo. 

Netherlands—Koninklijk Instituut van Ingenieurs, Sectie voor 
Regeltechniek, Prinsessegracht 23, Den Haag. 

Norway—Norsk Forening for Automatisering, Kranprinsensgate 17, 
Oslo. 

Poland—Polski Komitet Automatyki, Naczelna Organizacja Tech- 
niczna, ul. Czackiego 3/5, Warszawa. 

Roumania—Comisia de Automatizare, Academia Republicii Popu- 
lare Romine, Calea Victoriei 125, Bucharest. 

Sweden—Svenska Centralcommitten for internationella ingenjors- 
kongresser, c/o Svenska Teknologforeningen, Brunkebergstorg 20, 
Stockholm. 

Switzerland—Schweizerische Gesellschaft fiir Automatik, Stern 
wartstr. 7, Zurich. 

Turkey—Teknik Universitesi, Rektorlugu, Istanbul. 

Union of Soviet Socialist Republics—Nazionalny Komitet SSSR po 
Avtomaticheskomu Upravleniu, Kalantschevskaia ulitza 15a, 
Moskva. 

United Kingdom—British Conference on Automation and Compu- 
tation, c/o Institution of Electrical Engineers, Savoy Place, 
London W.C. 2. 

United States of America—American Automatic Control Council, 
c/o Mr. W. E. Vannah, Secretary, 330 West 42nd Street, New 
York 36, N. Y. 


» fue 
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METHOD OF OPERATING 


The business of IFAC is conducted by its General 
Assembly, its Executive Council, and its Committees 
as described below. 


General Assembly 


The General Assembly is the supreme body of the 
federation and consists of delegations of the member 
organizations, each member organization being entitled 
to one vote. It meets infrequently, generally upon the 
occasion of an international congress of IFAC. An ex- 
traordinary meeting of the General Assembly was held in 
Chicago on September 16-18 to review plans for the 
First International Congress of IFAC in Moscow, in 
1960, to elect officers and to modify the Constitution 
and Bylaws. 
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Executive Council 


The management of the federation is vested in an 
Executive Council consisting of the President, the First 
and Second Vice-Presidents, the Treasurer, six Ordinary 
Members, and the immediate Past President. In addi- 
tion, serving the Executive Committee and appointed 
by it are the Secretary and the Editor. The present in- 
cumbents of these offices are listed below. The Executive 
Council has been meeting annually to conduct the 
necessary business of directing IFAC. 


President—A. M. Letoy, Russia. 

First Vice-President—KE. C. Gerecke, Switzerland. 

Second Vice-President—O. Benedikt, Hungary. 

Treasurer—M. Cuenod, Switzerland. 

Ordinary Members—J. F. Coales, United Kingdom; G. Evangelisti, 
Italy; P. J. Nowacki, Poland; J. Balchen, Norway; K. Kaneshigi, 
Japan; Z. Trnka, Czechoslovakia. 

Past President—H. Chestnut, United States. 

Secretary—G. Ruppel, Germany. 

Editor—V. Broida, France. 


Commuttees 


The Executive Council is assisted in its work by an 
Advisory Committee and is entitled to establish Tech- 
nical and Special Committees to deal with special sub- 
jects. The Advisory Committee is international in na- 
ture and consists of not more than one representative 
from each member organization. It has the duty of ad- 
vising the Executive Council regarding the technical 
work of IFAC. Also the Advisory Committee shall rec- 
ommend Technical and Special committees as required, 
specifying their areas of work. D. P. Eckman, United 
States, is the Chairman of the Advisory Committee and 
J. Loeb, France, is Vice-Chairman. 

The following Technical Committees have been estab- 
lished with Chairmen at the addresses noted. 


A pplications—J. Mozley, Johns Hopkins University, Baltimore, 
Maryland. 

Bibliography—W. Oppelt, University of Darmstadt, Darmstadt, 
Germany. 

Componenis—G. Boromissza, Villamosgepek Uezemtant Tanszek, 
Budapesti Muszaki Egyetem, Stozcek-u. 2, Budapest XI, Hungary. 

Education—A. Marino, Consiglio Nazionale delle Ricerche, Piazzale 
delle Scienze, 7, Rome, Italy. 

Terminology—E. C. Gerecke, ETH, Schweizerische Gesellschaft fiir 
Automatik, Sternwartstr. 7, Zurich 6, Switzerland. 

Theory—B. N. Petrov, Institute of Automatics and Telemechanics, 
Academy of Science, 9 leningrader Chaussee, Moscow D40, USSR. 


AACC representatives to these committees include: 


Bibliography—T. J. Higgins, University of Wisconsin, Madison, 
Wisconsin. 

Components—D. J. Pompeo, Shell Development Corporation, Emery- 
ville, California. 

Education—Otto J. Smith, University of California, Berkeley, Cali- 
fornia. 

Terminology—H. L.. Mason, National Bureau of Standards, Wash- 
‘ington 25, D. C. 

Theory—John G. Truxal, Brooklyn Polytechnic Institute, Brooklyn, 
New York. 


The functions of these committees in their respective 
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service in the field of automatic control on an inter- 
national basis. Included in these programs in the future 
will be the preparations for Congresses and other special 
meetings. In addition, these committees will prepare 
periodic reviews of the technical state of the art in auto- 
matic control and associated fields. At present, detailed 
programs and plans are being prepared by the various 
committees with reports of progress to be presented at 
the Moscow Congress in June, 1960. 

A special committee chaired by Professor C. J. D. M. 
Verhagen was formed to present for IFAC a program on 
Instruments and Transducers in Automatic Control at 
the recent ISA Conference and Exhibit at Chicago in 
September, 1959. These presentations by nine speakers 
were well received and served to show how helpful inter- 
national meetings can be from the viewpoint of speaker 
and audience alike. Although the Technical and Special 
Committees should be small for effective action, each 
member organization may recommend members to the 
chairmen of these committees. The members of each 
committee are selected by the committee chairmen with 
the concurrence of the Advisory Committee. In addi- 
tion, these committees may invite other experts to offer 
their services in particular areas, as for example, on 
specific subcommittees. 


PRESENT STATUS AND FUTURE OUTLOOK 


IFAC has grown in membership and strength during 
its two years of existence and looks forward to a future 
of continued growth and greater effectiveness. Five 
issues of the IFAC bulletin have been published and 
continued publication on a three or four issues per year 
basis is planned by Victor Broida of France, the Editor. 
From the bulletin, engineers and scientists are obtaining 
an increasing awareness of what is going on interna- 
tionally in the field of automatic control. Publication 
in the bulletin by more member organizations is antici- 
pated, and more letters on technical matters by indi- 
viduals will be included. 

International exchange visits of experts fron both in- 
dustry and university have been made between a 
number of member organizations in IFAC and further 
arrangements are being planned for student and pro- 
fessor exchanges in the future. These exchanges should 
speed the development of new ideas more rapidly than 
heretofore. 

Foremost in IFAC’s accomplishments is the work 
being done in preparation for the Moscow Congress re- 
ported elsewhere in this issue. The work of the National 
Committee of the USSR for Automatic Control has been 
thorough and well organized. The completeness of these 
preparations stresses the fact that everything possible is 
being done to make this Congress an outstanding one. 
In addition, the work of paper solicitation and review 


being done by the various member organizations is in” 
many cases a major undertaking. Arrangements for 
% 


‘ 


areas include the preparation of programs of activity of 


1960 


preparation of an English version of the preprints and 
publication of the proceedings have been made with 
Butterworths Publishing Limited by Professors John F. 
Coales of England and John R. Ragazzini of the United 
States as co-editors. 

IFAC has established a sound working organization 
and a cordial relationship among individuals in its many 
member organizations. IFAC has provided a means 
whereby the science of automatic control can grow and 
flourish more rapidly and effectively. With its newly 
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modified Constitution and Bylaws, its fine plans for the 
Moscow Congress in 1960, and its sound foundation of 
committees, IFAC can look forward with confidence for 
its new officers to carry it to a more successful future. 

Automatic Control can play an important role in 
helping to bring about greater productivity and a better 
world in which to live. IFAC can help engineers and 
scientists in America as well as throughout the world to 
fulfill their part in making automatic control more 
effective. 


Dynamic Programming and Adaptive Processes: 
Mathematical Foundation" 
R. BELLMAN} anp R. KALABAT 


Summary—In many engineering, economic, biological, and sta- 
tistical control processes, a decision-making device is called upon 
to perform under various conditions of uncertainty regarding under- 
lying physical processes. These conditions range from complete 
knowledge to total ignorance. As the process unfolds, additional 
information may become available to the controlling element, which 
then has the possibility of ‘‘learning” to improve its performance 
based upon experience; i.e., the controlling element may adapt 
itself to its environment. 

On a grand scale, situations of this type occur in the development 
of physical theories throgh the mutual interplay of experimentation 
and theory; on a smaller scale they occur in connection with the 
design of learning servomechanisms and adaptive filters. 

The central purpose of this paper is to lay a foundation for the 
mathematical treatment of broad classes of such adaptive processes. 
_ This is accomplished through use of the concepts of dynamic pro- 
gramming. 

Subsequent papers will be devoted to specific applications in 
different fields and various theoretical extensions. 


I. INTRODUCTION 
T HE PURPOSE of this paper is to lay a foundation 


for a mathematical theory of a significant class of 
decision processes which have not as yet been 
studied in any generality. These processes, which will 
be described in some detail below, we shall call adaptive. 
They arise in practically all parts of statistical study, 
practically engulf the field of operations research, and 
play a paramount role in the current theory of stochas- 
tic control processes of electronic and mechanical origin. 
All three of these domains merge in the consideration of 
Be problems of communication theory. 


7 * Manuscript received by the PGAC, September 9, 1958; re- 
1 “ised manuscript received, February 12, 1959. 
uot The RAND Corporation, Santa Monica, Calif. 


Independently, theories governing the treatment of 
processes of this nature are essential for the understand- 
ing and development of automata and of machines that 
“learn.” 

We propose to illustrate how the theory of dynamic 
programming [1] can be used to formulate in precise 
terms a number of the complex and vexing questions 
that arise in these studies. Furthermore, the functional 
equation approach of dynamic programming enables us 
to treat some of these problems by analytic means and 
to resolve others, where direct analysis is stymied, by 
computational techniques. 

In this paper, general questions are treated in an ab- 
stract fashion. In subsequent papers, we shall apply the 
formal structure erected here to specific applications. 


II. ADAPTIVE PROCESSES 


We wish to study multi-stage decision processes, and 
processes which can be construed to be of this nature, 
for which we do not possess complete information. This 
lack of information takes various forms of which the 
following are typical. 

We may not be in possession of the entire set of ad- 
missible decisions; we may not know the effects of these 
decisions; we may not be aware of the duration of the 
processes and we may not even know the over-all pur- 
pose of the process. In any number of processes occur- 
ring in the real world, these are some of the difficulties 
we face. 

The basic problem is that of making decisions on the 
basis of the information that we do possess. An essential 
part of the problem is that of using this accumulated 
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knowledge to gain further insight into the structure of 
the processes, using analytic, computational, and exper- 
imental techniques. 

From this intuitive description of the types of prob- 
lems that we wish to consider, it is clear that we are 
impinging upon some of the fundamental areas of scien- 
tific research. Obvious as the existence of these problems 


are, it is not at all clear how questions of this nature can 


be formulated in precise terms. 

Particular processes of this type have been treated in 
a number of sources, such as the works on sequential 
analysis [14]; the theory of games [13]; the theory of 
multi-stage games;! and papers on “learning processes” 


[2], [3], [S}- [8], [11]. 


III. THE UNFOLDING OF A PHYSICAL PROCESS 


In order to appreciate the type of process we wish to 
consider, the problems we shall treat, the terminology 
we shall employ, and the methods we shall use, it is es- 
sential that we discuss, albeit in abstract terms, the be- 
havior of the conventional deterministic physical 
system. 

Let a system S be described at any time ¢ by a state 
vector p. Let t1, te,-++, be a sequence of times, 
ty<te< +++, at which the system is subject to a change 
which manifests itself in the form of a transformation. 
At time f;, $1 is converted into 7;(p1), at time fe, 
b2=T (fi) is converted into T2(p2), and so on, with the 


result that the sequence of states of the system is given 


by the sequence {pe}, where 


peri = Tx(pr), 


The state of the system at the end of time ¢y is then 
given by 


pyar = Ty(Ty-i(- - + T2(Ti(pi)) - - - )), (2) 


where #1 is the initial state of S. 

If 7,(p) is independent of k, which is to say, if the 
same transformation is applied repeatedly, then the 
preceding result can be written symbolically in the form 


pyar = TN(pi). (3) 


The interpretation of the behavior of a physical sys- 
tem over time as the iteration of a transformation was 
introduced by Poincaré, and extensively studied by 
G. D. Birkhoff [4] and others. It furnishes the back- 
ground for the application of modern abstract operator 
theory to the study of physical systems, as, for example, 
in quantum mechanics [12]. The idea of using this 
fundamental representation in connection with the for- 
mulation of the ergodic theorem is due to B. O. Koop- 
man. 


IV. FEEDBACK CONTROL 


With all this in mind, we are now able to introduce 
the concept of feedback control. 


See R. Bellman [1], ch. 10. 
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Supposing that the behavior of the system as de- 
scribed by the foregoing equations is not satisfactory, 
we propose to modify it by changing the character of 
the transformation acting upon p. This change will be 
made dependent upon the state of the system at the 
particular time the transformation is applied. 

In order to indicate the fact that we now have a 
choice of transformations, we write 7(p, g) in place of 
p. The variable g indicates the choice that is made. Con- 
sequently, we shall call it the control variable, as opposed 
to p, the state variable. To simplify the notation and dis- 
cussion, we shall assume that the set of admissible 
transformations does not vary with time. 

If g, denotes the choice of the control variable at 
time tx, we have, in place of (1), the relation 


Pir = T (pe; qk)» k= iS 2; Pee (4) 


with py4+1 explicitly determined as in (2). 

The associated variational problem is that of choosing 
91,92, - * * , gv so as to make the behavior of the system 
conform as closely as possible to some preassigned 
pattern. We wish, however, to do more than leave the 
problem in this vague format. 


V. CAUSALITY 


Turning back, from the moment, to the deterministic, 
uncontrolled process discussed in Section ITI, let us note 
that the state of the system at time f,41 is a function of 
the initial state of the system and the number of trans- 
formations that have been applied. Consequently, we 
may write 


peri = felpr), (5) 


where py; is the initial state of the system. 

For the sake of convenience, let us merely write p in 
place of p;. Then, the function f;,(p) is easily seen to 
satisfy the basic functional equation 


fnen(p) = fm(fn(P)); OC G ec (6) 


This is the fundamental semigroup property of dy- 
namical systems. 


VI. OPTIMALITY 


With the foregoing as a guide, let us see if we can 
formulate the feedback control process in the same 
terms. 

To illustrate the applicability of the functional equa- 
tion technique, let us consider a finite process, of N 
stages, where it is desired to maximize a preassigned 
function, ¢, of the final state of the system, py. This is 
often called a terminal control process. 

The variational problem may now be posed in the 
following terms: 


Max 


[q1,92,-+-,¢n] 


(py). (7) 


This maximum, which we shall assume exists, is again a 
function of the initial state, , and the duration of the 
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process. Let us then introduce the function defined for 
all states p and N=1, 2, ---, by the relation 


ies) 


f(b) = Max 4(py), (8) 
g 
where g represents the set [gi, g2, - ++, qv). 

Let us now introduce some additional terminology. A 
set of admissible choices of the g:, [q1, gz, «+ + , qv], will 
be called a policy; a policy which maximizes $(py41) 
will be called an optimal policy. 

In order to obtain a functional equation correspond- 
ing to (6), we invoke the 

PRINCIPLE OF OPTIMALITY: An optimal policy 
has the property that whatever the initial state and initial 
decision are, the remaining decisions must constitute an 
optimal policy with regard to the state resulting from the 
first decision. 

The mathematical transliteration of this statement is 
the functional relation 


fw(p) = Max fy-s(T (4, q1)), (9) 


Nee. ay *. +". WILD 


f,\(p) = Max ¢(T(p, qv). (10) 


Further discussion, and various existence and unique- 
ness theorems for the functions {fi(p) } and the associ- 
ated policies will be found in [1]. 

In this way, the calculus of variations is seen to be a 
part of an extension of the classical theory of iteration, 
and of semigroup theory. 


VII. StocHAstTic ELEMENTS 


In order to treat questions arising in the physical 
world in precise fashion, it is always necessary to make 
certain idealizations. Foremost among these is the as- 
sumption of known cause and effect, and perhaps, even 
that of cause and effect in itself. 

To treat physical processes in a more realistic way, we 
must take into account unknown causes and unknown 
effects. We find ourselves. in the ironical position of 
making precise what we mean by ignorance. 

At the present time, there exist a number of ap- 
proaches to this fundamental conundrum, all based 
upon the concept of a random variable. Building upon 
this foundation is the theory of games. 

We shall discuss here only the direct application of the 
concept of stochastic processes, leaving the game aspects 
for a later date. 

The theory of probability in a most ingenious fashion 
skirts the forbidden region of the unknown by ascribing 
to an unknown quantity a distribution of values ac- 
cording to certain law. Having taken this bold step, it is 


- further agreed that we shall measure performance not in 


terms of a single outcome, but in terms of an average 
taken over this distribution of values. Needless to add, 
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this artifice has béen amazingly successful in the analysis 
of physical processes ; é.g., Statistical mechanics, quan- 
tum mechanics. | 

Following this line of thought, we begin to take ac- 
count of nae effects by supposing that the result 
of a decision ¢ is not to transform p into a fixed state 
T(p, q), but rather to transform p into a stochastic 
vector z whose distribution function is dG(z, p, q), de- 
pendent upon both the initial vector p and the decision 
q. Let us further suppose that the purpose of the process 
is to maximize the expected value of a preassigned 
function, ¢, of the final state of the system. 

Before setting up the functional equation analogous 
to (9), let us review the course of the process. At the 
initial time, an initial decision g; is made, with the result 
that there is a new state p,, which is observed. On the 
basis of this information, a new decision, gz, is made, 
and so on. 

It is important to emphasize the great difference be- 
tween a feedback control process of this type, in which 
the g; are chosen stage-by-stage, and process in which 
the g; are chosen all at once, at some initial time. 

In the deterministic case, the two processes are equiv- 
alent, and it is only a matter of convenience whether 
we use one or the other formulation.? In the stochastic 
case, the two processes are equivalent only in certain 
special situations. We shall be concerned here only with 
the stage-by-stage choice. 

The analog of (10) is then 


[AOORIEE f o@)dCG, p, 0); (11) 


and that of (9) is 
fv(p) = Max f fy-aaGl, 6,0), N= 2,3,-+-4 (12) 
qd z 


This type of process has been discussed in some de- 
tail [1]. 


VIII. S—Econp LEVEL PROCESSES 


Fortunately for the mathematician interested in these 
processes, the tale does not end here! For in a number of 
significant applications, it cannot be safely assumed 
that the unknown quantities possess known distribution 
functions. ; 

In many cases, we must face the fact that we are 
dealing with more complex situations in which far less 
is known about the unknown quantities. For a discus- 
sion of the importance of these processes in the general 
theory of design and control, see McMillan [9]; for a 
discussion of the dangers and difficulties inherent in any 
mathematical treatment, see Zadeh [15]. 


2 This corresponds to the choice we have of describing a curve asa 
locus of points or as an envelope of tangents. 

3 The descriptive version of this equation, when no control is 
exerted, is, of course, the Chapman-Kolmogoroff equation, the sto- 
chastic analog of (6). 
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A first attempt in salvaging much of the structure al- 
ready erected is to assume that the unknown quantities 
possess fixed, but unknown, distribution functions. Re- 
garding deterministic processes as those of zeroth level, 
and the stochastic processes described in Section VII as 
first-level processes, we shall refer to these new stochas- 
tic processes as second-level processes. 

Although it is clear that we now possess a systematic 
method for constructing a hierarchy of mathematical 
models, we shall restrain ourselves in the remainder of 
this paper to the discussion of second-level processes. 


IX. ADDITIONAL ASSUMPTIONS 


Some further assumptions are required if we wish to 
proceed from this point to an analytic treatment. These 
are 


1) We possess an a priori estimate for the distribution 
function governing the physical state of the sys- 
tem, which, until further knowledge is acquired, 
we regard as the actual distribution. 

2) We possess a set of rules which tells us how to 
modify this. a priori distribution so as to obtain an 
a posteriori distribution when additional informa- 
tion is obtained. 

3) We possess an a priori estimate for the distribution 
functions governing the outcomes of decisions, 
which, until further knowledge is acquired, we re- 
gard as the actual distribution, and, as above, we 


know how to modify this in the light of subsequent ~ 


information. 


In this paper, we restrict ourselves to the case of 
known physical states. 

In formal terms, our state vector is now compounded 
of a point in phase space p and an information pattern, 
dG(z, p, g). As a result of a decision q;, there result the 
transformations 


po pi (observed) 
dG(z, p*, q) > dH (2, p*, 9; po, G, q1, pi) (hypothesized). (13) 


On the basis of these assumptions, and considering a 
control process which continues in time as described in 
Section VII, we wish to pose the problem of determining 
optimal policies. For the first time, we are considering 
adaptive processes significantly different from those of 
the usual deterministic or stochastic control process. 


X. FUNCTIONAL EQUATIONS FOR SECOND- 
~ LEVEL PROCESSES 


As before, we introduce the function 


fn(b; G(z, p*, q)) =the expected value of (py, Gy) (14) 
obtained using an optimal policy for an N-stage 
process starting in state (p, G). 


Depending upon the objectives of the process, only one 
or the other of py and Gy may enter into ¢. Examples of 
both extremes abound. 
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Arguing as in the preceding sections, we see that the 
basic recurrence relation is 


fu(p; G(z, p*, 9)) 


= Max f fv-(w; A(z, p*, q; p,G,q1, w))dG(w, Pp, q1), (15) 


q1 
formiV SOP oye 
Ap; G(z, ee q)) 


, with 


ives i} $(w, H(s, 6*, 93», G, qx, »))dG(w, p, qx). (16) 


These equations are quite useful in the derivation of 
existence and uniqueness theorems concerning optimal 
policies, return functions, and in ascertaining certain 
structural properties of optimal policies [1], [2]. 

If, however, we treat processes which are too complex 
for a direct analytic approach, as is invariably the case 
for realistic models, we wish to be able to fall back upon 
a computational solution. The occurrence of functions 
of functions, e.g. the sequence {fy(p; G)}, effectively 
prevents this. 


XI. FURTHER STRUCTURAL ASSUMPTIONS 


In order to reduce the foregoing equations to more 
manageable form, let us assume that the structure of the 
actual distribution is known, but that the uncertainty 
arises with regard to the values of certain parameters. 

At any stage of the process, in place of an a priori 
estimate, G(z, p, gq), for the distribution function, we 
suppose that we have an a priori estimate for the distri- 
bution function governing the unknown parameters. 
Again, a basic assumption is that this distribution func- 
tion exists. 

The functional equations that we derive are exactly 
as above, with the difference in meaning of the distribu- 
tion functions that we have just described. 


XII. REDUCTION FROM FUNCTIONALS TO FUNCTIONS 


We are now ready to take the decisive step of reducing 
f(p, G) from a functional to a function. 

It may happen, and we will give an example, that the 
change in the distribution function, from G(z, p, q) to 
H(z, p*, q; p, G, qi, W) is one that can be represented by 
a point transformation. This will be the case if G and H 
are both members of a family of distribution functions 
K(z; a) characterized by a vector parameter a. Thus, if 


G(z, p, q) = K(z, p, g; a) 


H(z, p*, 9; ’, G, ui, w) = K(z, p, 9; 8), (17) 
the change from G to H may be represented by 
B= V(p, a, 1, w). (18) 


Then we may write 


fu(, G(z, P; q)) = fv(d; a), (19) 


1960 
and (15) becomes 


ty(p; a) = Max ec B)dK(w; a). (20) 


The dependence upon q; is by way of (18). 
XIII. An ILLUSTRATIVE PRocEss— 


DETERMINISTIC VERSION 


Let us now show how these ideas may be applied to 
the study of control processes. Consider a discrete 
scalar recurrence relation of the form 


(21) 


Here wu, is the state variable and », is the control vari- 
able. Suppose that the sequence {v,} is to be chosen to 
minimize the function 


Unt1 = AUn + Vn, Oyo 


N 
|u| +50 m2, (22) 
k=1 


subject to the constraints 


belay 3 = 0,-- <5) W~ 1, (23) 


Although the precise analytic form of the criterion 
function is of little import as far as the present discus- 
sion is concerned, we have used specific functions to 
make the presentation as concrete as possible. Further- 
more, the defining equation need not be linear. 

This is a simple example of a deterministic control 
process. Introduce the sequence of functions defined by 
the relation 


fy(c) = Min [ ux | +5 > uw? |, (24) 
{07} k=1 


where WN takes on the values 1, 2,---, and ¢ any real 
value. 
Then 


f,(c) = Min [| ac + 0 | + b(ac + w)?], (25) 


| vo] <r 
and for N>2, the principle of optimality yields the 
relation 


fy(c) = Min [b(ac + 00)? + f—i(ac + v0) |. 


|vo| <7 


(26) 


XIV. STOCHASTIC VERSION 


In place of the recurrence relation of (21), let us intro- 
duce a stochastic transformation 


Une) = Wn + Tn + Pnj uo = C. (27) 


Tere {rn} is a sequence of independent random varia- 

bles assuming only the values 1 and 0. Let 
tn = 1. with probability p 

= 0 with probability 1 — p. (28) 


The quantity p is known, and for simplicity taken to be 
independent of 7, although this is not necessary. 
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We now wish to minimize the expected value of the 
quantity appearing in (22). This is now a stochastic 
control process of the type described above in general 
terms. Call the minimum expected value fy(c). Then, 
following the procedures of Section VII, we have the 
relations 


min | f [| ac+ +10 


Jool <r T9 


fi(c) = 


+ b(ac + ro + 1146 (7) | 


Min 


|vo| <, 


| Pl ac + v6 + 1| + b(ac + v9 + 1)2| 


+(1—- P| | ac + v0 | + b(ac + oy | ; (29) 


and, for general N, 


in (c) = Min 


|vo] <, 


scant os Sein fara eee v0) | (30) 


| elb(ac + v9 + 1)? + fy-s(ae + v0 + 1)] 


XV. ADAPTIVE CONTROL VERSION 


Let us now consider the adaptive control version. We 
are given the information that the random variables r, 
possess distributions of the special type described above, 
but we do not know the precise value of p. 

We shall assume, however, that we do possess an 
a priori distribution for the value of p, dG(p), and that 
we possess a known rule for modifying this a priorz dis- 
tribution on the basis of the observations that are made 
as the process unfolds. 

If we observe that over the past m-+n stages, the 
random variables have taken on m values of 1 and 
values of 0, we take as our new a priori distribution the 
function 


dGnn(6) = pl = pyrac(e) / f PMA — p)°aG(), (31) 


a Bayes approach. 

Once we have fixed upon a choice of G(p), the a priori 
distribution function at any stage of the process is 
uniquely determined from the foregoing by the numbers 
m and n. This simple observation enables us to reduce 
the information pattern from that of the specification 
of a number or vector, in general, plus a function 
Gm,.n(p), to that of the specification of three numbers, ¢ 
and the two integers m and n. 

In this way, we reduce the problem from one requiring 
the use of functionals to one utilizing only functions. 
This is an essential step not only for computational pur- 
poses, but for analytic purposes as well. 

Let us then introduce the sequence of functions 
{ fu(c, m, n)} defined once again as the minimum ex- 
pected value of the quantity in (22), starting with the 

4 This is an assumption of the type called for in Section IX. Al- 


though reasonable, it is not the only one possible. There are analytical 
advantages in choosing G to be a beta distribution. 
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information pattern of m ones and » zeros, and state c. 
Then 


fi(c, m,n) = Min 


|vo| <r 


Ly fey well {Costa Naa Weeks 1 , (32) 


auc + vp 1)2-- | actu+1 | | 


where pm,» is the expected probability using the proba- 
bility distribution in (31), 2.e., 


[emma = prracio) 


Pn (33) 


Me 1 
f pmct = precip) 
0 
For N>2, we have the recurrence relation 


fy(c, m,n) = Min 
lvol <r 


auc + v) + 1)? 


+ fy-hac — 09-1, 9-1, 2)! 
Se ONE ee) Croats) 
+ fv-i(ac + v0, m,n + 1] . (34) 


In this fashion, we obtain a computational approach 
to processes with general criteria and an analytic ap- 
proach to processes with criteria of particular type. A 
thoroughgoing discussion of the analytic aspects of the 
solution of processes of this nature described by linear 
equations and quadratic criteria will be found in a forth- 
coming doctoral thesis by Marshall Freimer. Previous 
applications of these techniques have been made in [2 | 


and [3]. 
XVI. SUFFICIENT STATISTICS 


The fact that the past history of the process described 
in the preceding paragraphs can be compressed in the 
indicated fashion, so that functions rather than func- 
tionals occur, is a particular instance of the power of the 
theory of “sufficient statistics” [10]. 

Many further applications of this important concept 
will be found in the thesis by Freimer mentioned above. 

In a number of cases, this compression of data occurs 
asymptotically as the process continues; e.g., the central 
limit theorem. A number of quite interesting questions 
arise from this observation. 


XVII. Discussion 


In the foregoing pages, we have attempted to con- 
struct a mathematical foundation for the study of the 
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many fascinating aspects of the field of adaptive control. 
In further papers, we shall discuss a number of complex 
problems which arise from this approach. 

From the purely mathematical point of view, we are 
now able to contemplate a theory of continuous control 
processes of adaptive type, obtained as a limiting form 
of the theory of discrete control processes. A variety of 
significant convergence questions are encountered in 
this way. 

Furthermore, we can construct a theory of multi- 
stage games on the same foundations. 

Finally, the problem of computational solution is by 
no means routine, and there are a variety of interesting 
approaches based upon approximations in function 
space and approximations in policy space to be explored. 

From the conceptual point of view, we must face the 
fact that there are many further uncertainties to be 
examined in the state of the system, in the observation 
of the random effect, in the transmission of the control 
signal, in the duration of the process, and even in the 
criterion function itself. 
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The Properties and Methods for Computation of 
Exponentially-Mapped-Past 
Statistical Variables” 
JOSEPH OTTERMAN} 


Summary—The exponentially-mapped-past (emp) statistical 
variables represent an approach to the statistical analysis of a 
process when the interest is focused on the recent behavior of the 
process. An exponential weighting function, decreasing into the past, 
in the case of continuously observed processes, and a geometric 
ratio, in the case of discrete data, are utilized. This approach is the 
simplest from the point of view of ease of computation, and at the 
same time it possesses the advantage of some simple theoretical 
relationships, which are discussed. Analog computer circuits and 
digital computer flow diagrams which serve to compute the ex- 
ponentially-mapped-past statistical variables are presented. 


INTRODUCTION 


HIS PAPER is concerned with an approach to the 
oan analysis of a process when the interest 

is focused on the recent behavior of the process. 
This approach offers the advantage of simplicity in com- 
putation when the analysis is carried out by computers; 
and at the same time, it possesses the advantage of some 
simple theoretical relationships. It is an extension of the 
work reported earlier.! The material is primarily theo- 
retical, but analog computer circuits and digital com- 
puter flow diagrams which serve to compute the ex- 
ponentially-mapped-past statistical variables are pre- 
sented. The computational advantages, which have been 
realized by others,’ are pointed out. 

The exponentially-mapped-past statistical variables 
(which will be referred to as emp variables) are quanti- 
ties relating to a set of observations computed in such a 
way that the recent values of the observations contrib- 
ute more strongly than the values observed in the more 
distant past. The relative weighting is a geometric ratio 
in the case of discrete (naturally discrete or sampled) 
data and an exponential function in the case of con- 
tinuously observed functions. The emp variables are 
meant to apply to control or detection problems; there- 
fore, information about the process up to the present 
time only is utilized. The emp variables themselves are 
functions of time. 


* Manuscript received by the PGAC, September 8, 1959. This 
work was conducted by Project Michigan under Dept. of the Army 
contract (DA-36-039-SC-78801), administered by the U. S. Army 
Signal Corps. 

+ International Telephone and Telegraph Labs., Nutley, N. J., 
located at Bell Telephone Labs., Whippany, N. J.; formerly at the 
University of Michigan, Willow Run Labs., Ann Arbor, Mich. 

1 J. Otterman, “On Exponentially-Mapped-Past Statistical Vari- 
ables (Preliminary Report),” presented at the Joint Meeting of the 
Assoc. Computing Mach. and the Inst. of Math. Statistics, Cleveland, 
Ohio, April 3, 1959. (Unpublished. ) : 

2M. Blum, “On exponential digital filters,” J. Assoc. Computing 
Mach., vol. 6, pp. 283-304; April, 1959. 
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ContTINuOUS EMP VARIABLES 
EMP Average 


The emp average, fa(0), of a continuously observed 
variable’ f(¢) will be defined as 


7.(0) = a f f(dertdt, (1) 


where the past extends from 0 to — », and 1/aisa time 
constant for which the constancy of the regime can be 
effectively assumed. The coefficient @ in front of the 
integral sign comes to effect a normalized weighting 


function, since 
0 
if edt = 1/a. 
—o 


The concept of fa(t) is the average of f(r), determined 
experimentally with information about f(r) up to the 
time ¢ only. For an arbitrary moment t (value negative) 
in the past, the emp average will be given by 


(2) 


oe-2! {t “Ghee 


—3o 


LAD) = a f f(r)et 9dr 


—o 


Eq. (3) can be expressed in another way as 


Fal) = af flr)en-u(t = rar, (4) 
where u(é) is a unit step function [u(#)=1 for t>0, 
u(t) =0 for t<0]. This is a convolution integral of f(r) 
and an impulse response of a low-pass RC filter, or a 
“leaky” integrator, with leakage time constant 
RC=1/a. The computation of the emp average in ac- 
cordance with (3) or (4) is shown in Fig. 1. 


c=1 


R=2 


Fig. 1—Computation of fat) in accordance with (3). 


The above definitions apply in the case in which ob- 
servations extending indefinitely into the past can be 
assumed; 7.e., it is not necessary to have a good estimate 


3 The observed variable itself need not be continuous, but can 
assume discrete values only. 
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of the average during the initial stages of the analysis or 
control process. If such is not the case, we have (T desig- 
nates the beginning of the observation relative to the 
present and is negative) 


f.(0, T) = mar toe (S) 


and for an arbitrary time in the past 


hops a5 finer (6) 


This definition is adopted, since 


t ext aaa ett 
if etdr = ———_- (7) 
fy a 


It follows that if (3) is used instead of the more accurate 
(5), in the initial stages of observation when e?” is not 
negligible the absolute value of the emp average will be 
smaller than the “true” value, 7.e., the value resulting 
from the use of (5). The computation of the emp average 
in accordance with (5) is shown in Fig. 2. 


EMP Variance 


The emp variance s,2(0) is defined as follows: 


20) = af [f@) — Fa Pena (8) 


—oO 


and for time ¢ in the past the emp variance will be 


Sart) 


a {_[flr) — Jalr)Pere-Par 


ae—*t ip [f(r) — fa(r) |?ee7dr. (9) 


The scheme for computing s,?(¢) in accordance with this 
equation is shown in Fig. 3. This is based on the assump- 
tion that accurate computation in the initial stages is 
unimportant. If accurate computation at this stage is 
necessary, then the emp variance will be defined as 


6D = oe | Me) —Oerar, (10) 


as was done for the average in (5). For the remainder 
of the section, (8) only will be used as defining the emp 
variance. Introducing the value of fa(t) as defined by (3) 
into (8) we obtain 


Sa°(O) = af | — ae “f ik flejear | Caddie, = € |p) 


After expanding the right side of the above equation, 


$0) =a fe PWertdt — 202 f : flat ie f(aertdr 


eas fia tid floderar| 


(12) 
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Fig. 3—Computation of s,2(¢) in accordance with (8). 


Theorem: A simple relationship, parallel to one exist- 
ing in ordinary statistical parameters, will be now de- 
rived. The relationship involves the variable f.7(0), 
which is defined as 

0 
f2(0) =a f fd edt. (13) 


It will be shown that for s,2(0), as defined by (8), the 
following holds: 


= fa®(0) — fa2(0). (14) 


To prove this, the last expression in (12) is integrated 
by parts. The integrand is regarded as a product of 


| ff fear 
and, exc! 


fe 
-<| i sloderrtr | 
te =f « evar | fs oe 
fimo} 


ele 2 penne [sever 


Sa’(0) 


0 


(15) 


The above holds, since the lower limit, i>— ©», of the 


expression 
eat t 2 
— | if feerdr | 
a —o 
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vanishes, if f(r) has the property of any measured vari- 
able, viz., that there exists a number A such that 


a f ore — 20 f foa f Po 
refer f see} 

x f_poretas— 208 soa f pera 
AL] 

+f poa f" soerar 


fa?(0) — fa2(0), 


which constitutes the desired result. The scheme for 
computation of s.(¢) in accordance with (16) is shown in 
Fig. 4. 


Sa’(0) 


I 


(16) 


EMP Fourier Transforms 


Fourier transform F(jw) of f(¢) represents an analysis 
of the presence of different frequency components in 
the phenomenon under observation. 

The emp Fourier transform is defined as follows for 
an arbitrary moment ¢ in the past: 

t 
FAjie) =e f Flayera 2 e-ivedey, (17) 
The extension of the integration of + is permissible 


if f(t) is assumed to vanish for x>t, 1.e., for the future 
relative to the time of interest. Thus, 


Fi(jw) = af Haden — x)etdx, (18) 


where z(t) is a step function described before. 
Examination of the above equations shows that 


F,(jw) = aF(—a + jw); 


e., the emp. Fourier transform is not new but is the 
ordinary Fourier transform evaluated at a line, —a+ju, 
parallel to the jw axis, and multiplied by a. But the 
presentation here brings out the fact that aF(—a+jw) 
represents an analysis for recent frequency-wise be- 
havior of the observed variable, the “recentness” in- 
creasing with increasing a. It will be noted that 
F.(0) =fa(0); thus the emp average is equal to the 
dc (w=0) value of the emp Fourier transform. 

It may be noted that, if it is desired to concentrate on 
analysis in a certain region in the past, 7.e., to view the 
_ past through an aperture e@'—e! (plotted for a =2a 
in Fig. 5), then the frequency analysis is given by 
F,,(jw) — Fa,(jw). This statement actually will apply to 


— 


(19) 
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Fig. 4—Computation of 5,%(#) in accordance with (16). 


> 
t 


Fig. 5—Aperture in the form e“t—e (¢ is negative). 


other emp variables, such as the emp average or the 
emp variance. The problem outlined in this paragraph 
is rather different from the question of the determina- 
tion of the best recent analysis of the observed variable, 
for it is more related to interpolation than to extrapola- 
tion. Hence, it will not be discussed further in this paper. 


EMP Density of Occurrence 


Let f(t) be a random variable under observation. 
Assuming a stationary and ergodic process, let us con- 
cern ourselves with the probability, experimentally de- 
termined, that f(t) lies between two values; or more 
specifically, the probability that 


(a- 2) <n <(a+). 


This probability will be denoted as p(fi)Af. The im- 
portant question of whether a limit for p(f1) exists as 
Af-—0 need not concern us in this discussion. 

This probability can be defined on the basis of obser- 
vation of f(#) during a time interval T in the following 
manner: 


‘a(dt 


GAY pee see 


= (20) 


in which 6(¢)= 1 if 
*)<m<(e% 
( ey 2 i Th 2 ) 
and 6(t) =0 elsewhere. 


4 The use of the word “frequency” in the statistical sense will be 
avoided in this paper. 
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A parallel definition of emp density of occurrence 
p.(f1) in a slowly changing process is first made for an 
infinite interval of observation (in (20), 7-0), thus: 


0 

pal fr Af = af A(t) e%'dl, (21) 
in which 6(t) is defined as before. When summation of 
the density of occurrence is carried out over all the possi- 
ble range of f(#), without overlapping of the intervals, 
the summation will result in unity. Thus, p.(fi)Af has 
the basic property of a probability function, in that it 
ranges from zero to unity. Its physical significance is 
obvious from its mathematical definition: the experi- 
mentally determined density of occurrence is computed 
in such a way that the recent interval during which the 
observed function lies between (f; —Af/2) and (f:+Af/2) 
contributes more strongly than a similar interval in the 
more distant past. 

When an infinite interval of observation cannot be 
assumed, the emp density of occurrence will be com- 
puted according to the formula 


a 0 
plfaf=—“— f opera, (22) 
ip aie etl . 
and for an arbitrary time in the past, 
a t 
pal fi Af = —*_{ O(r)e*"dr. (23) 
ext — ect T 


EMP AuTOCORRELATION AND EMP 
POWER SPECTRUM 


The discussion here is almost identical to the concepts 
presented by Fano nine years ago.’ The difference lies 
primarily in the fact that the concept of the emp Fourier 
transform is used. The result is that Wiener’s theorem 
is reformulated exactly in the emp language, 12.e., the 
emp autocorrelation and the emp power spectrum form 
an exact emp Fourier pair. There is a difference in the 
coefficient of the emp power spectrum as compared with 
Fano’s corresponding function. The other, rather in- 
significant, difference is that active analog computer 
circuits rather than passive networks are presented. 


5 R. M. Fano, “Short-time autocorrelation functions and power 
spectra,” J. Acoust. Soc. Amer., vol. 22, pp. 546-550; September, 
1950. The writer did most of his work unaware of the work by Fano, 
but this paper has been modified in light of his work. Specifically, the 
writer originally considered the function Wo(r) rather than ¢,(r) as 
the emp autocorrelation function. Thus, the writer considered that 
the emp autocorrelation function should contain the e~** term. The 
notation in this paper has been revised to conform almost exactly 
with Fano’s notation. The exceptions are as follows: 


a) The subscript @ is used rather than the subscript ¢ to denote 
the emp variables. This was retained from the author’s original 
notation since it is anticipated that in the author’s later work 
comparison of emp variables with different time constants, 
1.€., a, and ae, will be involved. 

b) Fa(w) denotes the emp Fourier transform of f(t). The use of 
F(x) in Fano’s sense [his (14)| is continued. 


Fano’s text is followed closely in various paragraphs of this section. 
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The emp autocorrelation function of f(t) is defined as 
follows: 
t 
$.(r) = 2a f Kai@ — eet dx, (24) 
where 7 is a positive delay. The analog computer scheme 


for obtaining the emp autocorrelation function is shown 
1 Pao 0. 


Fig. 6—Computation of ¢.(r) in accordance with (24). 


Consider now the function 


Yale) = 2ae~ fflayfe— 7) 


-e a(t) y(t — x)u(t — x — r)dx, 


(25) 


where u(t) is a step function described before. This func- 
tion is symmetrical in 7, as can be seen by substituting 
y=x—r. Let the function ¢,(7) for negative 7 be inter- 
preted as a symmetrical image of the integral in (24) 
taken for positive 7. Then it is possible to write: 


Wa(r) a da(T) gal, 


Note that W.(r) represent a convolution integral of 
the function, F;(«) (here Fano’s t subscript is retained): 


(26) 


F(x) = (2a) *f(xjee" ult — x), (27) 
and of F;(—x). Therefore 
1 co 
vals) = 5 ftele)eat(wetrde, (28) 


where ga*(w) denotes the complex conjugate of ga(w), 
and where 


gle) = [ademas 


—o 


(2a) f fayeout — x)e—*dx 


9) 1/2 t 
= (=) a { fie rear 
a 00 


This last expression is equal to (2/a)!/? times Fy(w), the 
emp Fourier transform of f(x), thus 


(29) 


A AEE Gir Nh ha 
Walt) = a aie a(w) By*(w) eT da 


iy) 00 


2K: id ey, 


P.(co)e*"deo, (30) 
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where P.(w), which defines the emp power spectrum, 
is given by 


P.Aw) = Fa(w)Fa*(w) = | Fa(w) |?. (31) 
Eq. (30) can be rewritten in view of (26) 
adba(r)e—al7l 1 = ; 
: = x J Perera. (32) 
Taking the inverse Fourier transform of (32), 
a Sy ; 
Lie) =f a(t) et! e—irdr 
0 
= a f ba(r) e%e-1dr 
= emp Fourier transform of ¢a(r). (33) 


Wiener’s theorem is thus restated in emp term. 

Coming now to the question of experimental determi- 
nation of P,(w), it is noted that F.(w) can be written as 
follows: 


F.(w) act f f(aje* u(t — x)eie (ta) dy 


= cial f fadenrul — x) cosw(t — x)dx 


+5 f f(x)e* u(t — x) sin w(t — sds | . (34) 
It follows that 


P.(w) = F(w) |2 = al f f(a)e*“™ cos w(t — “ds | 


-- | [ f(aje*) sin w(t — nds | } (35) 


The first term of the right-hand side of (35) represents 
the square of the output of a filter with the impulse 
response 


E,|(t) = ae~*' cos wt u(t), (36) 


and with f(#) as the input. The transfer function of such 
filter is given by [a(a+ 6) |/[(a+j6)?+w?], where 6 is 
the frequency variable. The second term of (35) is 
similarly recognized as the square of the output of a 
filter with the impulse response 


E./'(t) = ae~* sin wt u(t). (37) 
The transfer function of such filter is given by 
aww / [(a +76) +0]. 

The computation of P.(w) by an analog computer 
setup is shown in Fig. 7. 
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Fig. 7—Computation of Pa(w) in accordance with (35). 


DISCRETE EMP VARIABLES 
EMP Average 


The emp average of a variable I() observed at dis- 
crete intervals A? is defined numerically by 


T,(0) = (1 — g)[Z(O) + T(1)g + 1(2)g? +--+ T(n)gr-- + | 


=(1-—g Di Inq", 


n=0 


N—> 2, 


(38) 


where the 0 in 7,(0) indicates that the summation ex- 
tends to the most recent observation, namely J(0). The 
factor 1—g is required to normalize the weighting func- 
tion, since 

N 1 — Nt 


Sere ea 


fae (39) 
and g¥*1 gradually vanishes as NV becomes large. 

The physical meaning of 7,(0) is apparent from its 
numerical definition. It is an average of the observed 
variable J, where the recent values are weighted more 
heavily, and the input of a long time ago contributes in- 
significantly to the average. How soon the past inputs 
stop contributing significantly to the average depends 
on q. Since g generally will be very close to 1, we can 
say that 1—q=e is a small quantity, and 


gla-a = (1 — de t/e, (40) 


Thus, an input which bears a number n=1/(1—q), 2.., 
an input that occurred approximately At/(1—q) before 
the present, contributes by a factor e less than the most 
recent input. The period At/(1—q) should be made 
smaller than the time of effective constancy of the 
regime. Within this restriction 1/(1—q) should be made 
large, 2.e., g should be as close to 1 as possible, in order 
to get an effective contribution to the average from a 
larger number of data. 

The average 7,(0) will be smaller in absolute value 
than the “true” average [as defined later in (41) | during 
the initial stages of the computation, during the stages 
when q+ is not negligible. When it is desired to have a 
reasonable estimate of the average during the initial 
stages, the average can be defined numerically as fol- 
lows: 
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x Le 
Tacos) = [egg HOO) + 1C)9 + LDQ" + TCV)" 
—q N 
= rar 2d I(n)q” (41) 


The flow diagram for computing the emp average in 
accordance with (39) is shown in Fig. 8. That corre- 
sponding to (41) is shown in Fig. 9. It will be seen that 
the second approach is somewhat more complicated to 
compute, and for the remainder of the report it will be 
assumed that the need for an accurate average during 
the initial stages does not arise. 

For a moment in the past, the emp average 7,(n) is 
numerically defined 


— 
Q 
en 
= 
= 
ll 


i gems (ee 
—— Di l(h)g', Nox. (42) 
Of Se 


This is a simple extension of 7,(0), with the difference 
that J,(m) represents an average of some time ago. 


EMP Variance 


The emp variance s,?(0) is defined in (43) and again 
in (44): 


Sq (0) ="{1'2tq) du [Z(m) — Ta(n) }?q"; 


5q°(0) si > | 1 


ana | 
10 vepiako q” 
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N 
MULTIPLY Dir a" BY 
n=1 


N 
n-1 
MULTIPLY Diuma 
n=1 


BYq - 1 - q TO FORM I q(0) 


Fig. 8—Computation of 7,(0). 


N N 
> I(n)g* becomes >> I(n)g"4, 
n=0 n=l 


once the new J(0) has been read. Initially, 0 is stored in the mem- 
ory cell containing 


N 
dX L(n)qr. 


n=0 


ADD 1(0) TO 
N 
ye I(n) q” 
n=l 
N 
MULTIPLY ma I(n)q" 
n=l 


1-q = 
BY : wal TO FORM Iq (0,N) 


Fig. 9—Computation of 7,(0, NV). g¥+1 becomes g¥, once the new I(0) 
has been read. Initially 1 is stored in the memory cell containing 
qv. 


(43) 


N 2 
ira eal 5. r(a)9| fa 


k=n 


= D1) = = QU) + Ht tg $+ + 1h" + + Wor 


n=0 


ie} 


= Di tl(n)q— (1 — DL + Ig + 1m + 2g? + 


n=0 


- ++ +b T(k)ghm +. +. J} eqn 


= {10 q¢—-(— [7+ 12qt---+ gr +--- JP 


+ {1(t)¢ — 1 — g)[Z(2)q + 1(3)g2 + + (ng + ++ +d hq 


+ {I(n — 3)qg — (1 — g)[Z(m = 2)g + Tn — 1g? + T(n)gi +--+ hegre 
+ {L(m — 2)q — A — QT — 1g + Tg? t+ --- }}%gr 
+ {I(m—qg—-A—g[Zi)qt -- + ]}2gr 
Fy eu a) N20 


(44) 
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TABLE I 


eq 9) 1(0)geg® 
Sar Oly=\¢)11)o*F¢ 


= — 29(1 = gn = 3)¢8 qn 
2g 4 yg) el 2) gg? 
= 241 — g)l(n — 1)gq" 


+ ¢I(n)q” 


*{T(n)q?r + ee —1)gr1...27 


(1) ntl] 0 
LT (n)g? + +» » 27()qrIq! iy 


rays) 
&~ 


(i aI (n)q® e a(n - 1)g8 a2 21(n — 2)q'|qr 
(eg) ng? 21 (n = 1)g*lg” 
(1 — g)*{Z(n)@?la” 


The terms with J(m) that contain I(k), k<n, are pre- 
sented in Table I. The terms are culled from (44) above 
in such a way that each row in Table I corresponds to 
an appropriate row in (44). In the first column of Table 
I (which contains only one expression) is the squared 
term of the first term of each row of (44); in the second 
column of the table are the cross-terms of the first term 
and the terms in brackets of (44); and in the third col- 
umn are the terms of multiplication inside the brackets. 

Table I can be rewritten using the formula for the 
sum of geometric series, as follows: 


n—1 
I(n)grt? — >) 2gr*(1 — g)I(k) — 2g"*(1 — g)I(0) 
k=1 
ae (1 ats g)2I(m) (qrt} a gtt? =. mart + ge") 
n—1 
Beebe eget. gett 2) get") 


k=1 
n—1 

= I(n)grt? — Dd) 2gr*(1 — g)I(k) 
k=1 


1 — gq” 
+ (1 — g)Z(n) qr — 


m1 1 leas 
— 2gr*H(1 — g)T(0) + 21 — 9)? 2 T(h)gtt 5 — 


k=1 


= gt [I(n)q + I(n) — I(n)q — T(n)qr + Ing? 


— 2(1 — HO) = 21 = a) & ree". (45) 
Theorem: Let I? be defined as follows: 
720) = (1-9 VPOe. (46) 


n=0 


The following relationship will now be demonstrated: 


g50'(0) = T2(0) — Ts (0). (47) 
Consider the expression 7,2(0) —T,(0): 
7,40) — 7,(0) 
= d-9Srmr-a-9{ Lime ]s os 
n=0 n=0 


a a -— 


and, after dividing by (1—gq), expand as follows: 


72(0) — 7, (0) 
ro =PO)+PH¢+ Pe2+--- 
+ P(n)gr+---—(—g[Z0) + 1g 
+ £(2)9? + > »< I(m)q™ +. « (2, (49) 


The terms with J(”) which contain only I(k), such that 
k<n, are 


I(nm)qr — 21 — g)[Z(0) + (1g + ++ ++ Tn — 1)gr-Jq" 
— (1 — g)I(n)q? = qr[I(n) — I(n)qr + I(n)gr*! 


ott g)I (Oi 2(1) 9) ErH0'|. 


k=1 


(50) 


A comparison of terms of (50) with terms given in (45) 
establishes the validity of (47). 


CONCLUSION 


The exponentially-mapped-past statistical variables 
offer considerable advantages in either digital or analog 
computation. The analog computer circuits are very 
simple and easy to realize. Neither perfect integrators 
nor delay devices are required, elements which are 
necessary for instance if statistical analysis is based on 
an equally-weighted finite aperture extending from the 
present to time Jin the past. The digital computer pro- 
grams involve storing only one or several emp variables, 
rather than storing all the data for the interval of 
interest. 

The theoretical relationships which have been pre- 
sented are simple and show parallelism with ordinary 
statistical variables. 

While this report does not present a complete theory 
of exponentially-mapped-past statistical variables, it 
attempts to present a unified approach to the emp sta- 
tistics. 
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Generalized Weighting Function and Restricted 
Stability of a Linear Pulse-Modulated 
Error Feedback System‘ 


WILLIAM A. JANOS+ 


Summary—A generalized weighting function is obtained for a 
linear feedback system with a pulse-modulated error signal. This is 
expressed in the form of a matrix operator acting on an input vector, 
the components of which are the first R—1 derivatives of the input, 
where R is the order of the unmodulated closed-loop system. In 
addition, since the system operator takes the form of a finite dimen- 
sional matrix, it has been possible to make and realize more stringent 
conditions on the transient stability; namely, a preassigned bounded 
output after a preassigned time. 

Although the solution has been obtained generally for non- 
uniform pulsing, the latter stability investigation has been made 
only for the uniform case. 


HE present article is concerned with the inversion 

of the integral equation relating output to inputed 

in a linear feedback system with a pulse-modulated 
error signal. Both uniform and nonuniform pulsing are 
considered. The formal solution is expressible in closed 
form as a matrix operator of rank equal to the order of 
the closed-loop system, acting on an appropriately de- 
fined input vector. 

Previous work on this problem has been carried out 
primarily by Farmanfarma.! In such an investigation 
the analysis was performed in the transform (s) domain, 
to yield a comprehensive treatment of the problem in 
terms of the P-transform. Stability studies of an asymp- 
totic nature were also included in the latter work. 

Although the Farmanfarma study is most commenda- 
ble and penetrating, it seems perhaps even formally un- 
wieldy. The present treatment essentially obtains a 
generalized weighting function for the system in the 
form of a matrix operator acting on an input vector, the 
components of which are the first R= 1 derivatives of the 
input, where R is the order of the unmodulated closed- 
loop system. Thus, the present approach may have cer- 
tain pictorial merit. Although it is formally more con- 
cise, the calculational work necessary for an arbitrarily 
ordered system can become rather involved too. In addi- 
tion, since the system takes the form of a finite dimen- 
sional matrix, it has been possible to make and realize 
more stringent conditions on the transient stability; 
namely, a preassigned bounded output after a pre- 


* Manuscript received by the PGAC, March 3, 1959; revised 
manuscript received August 17, 1959. 

+ Convair (Astronautics) Division of General Dynamics Corp., 
San Diego, Calif. 

1G. Farmanfarma, “General Analysis of Finite Pulsed Linear 
Systems,” Ph.D. dissertation, University of California, Berkeley; 
November, 1957. 


assigned time. Although the solution has been obtained 
generally for nonuniform pulsing, the latter stability 
investigation has been made only for the uniform case. 

The following system is given (see Fig. 1) where G, 


Fig, 1—System block diagram. 


and G2 label the components with corresponding impulse 
responses Gi(r) and G2(7), respectively. The pulse 
modulation is represented by the crossed switch of dura- 
tion h, over the nth interval. The integral equation 
governing this system has the form, 


f ete ule) fi . ip ates ry fle’yan'h de 


=f) (1) 
where 
ur) =1, Tat Tit%,, allintegers n, 
= 0, otherwise; (2) 
fi is assumed turned on at T=0. Let 
aD serthe Vanlete Gantt f Gi( 
ace transform 
Di(p) p 0 1 cae 
NA) the transfo { G2(r) 3 
rm o : 
Dalb) i i 


Further, let us make the following definitions, for all 
functions f concerned: 


fr(r) = f(r), St Sy He hy. 
= 0, otherwise. 

BAD) = WT nat hy Sir (4) 
= 0, otherwise. 


Note that the transforms of f,, f, would be the P-trans- 
forms of f as expressed previously.! 
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Hence, over the interval (T,, Tn+hn) we have the 
following operational form of (1). 


Ny 
tone + T) SS Sins + tT’) = 


N2(p) 


On fe 7! 
Pek } 


OS etn fg C5) 
where the integral operators, 

Nie), Nilo) Nal) 

Di(p) Di(p) Ds(p) 


refer to the variables 7’, r’’. 
After transposing and integrating by parts, (5) takes 
the form 


fon(Tn + 7) = {Ni(p)N2(p) + Di(p) Do(p)}> 
‘ { D2(p)Ni(p)fin( Tn si T’) 


R m—1 
+ Drm Ye Fon(Tr)} (6) 
m=1 r=0 
where 
to (Tn) = (Fooled) 
PE a TT, 
and 
R 


Ni(p)N2(p) + Di(p) Do(p) = 


m= 


Ymp™; (7) 
0 
1.€., Ym is the coefficient of p” in the polynomial in p, 
N,N2+D,D2. Notice that the last term on the right of 
(6) represents the influence of the initial conditions at 
7=T,+ ,which are the terminal conditions of fon: How- 
ever, during the preceding interval, (7,1+,-1,T,), we 
have 


fon(T a1 - Vong. 3 + 7’) 


1 Q py m—l1—r (r) \ 
oS ve gxthile toe 
Dp = B {> p Son Ta-1 1) 
0) =< ql! Ss hss — sf AS i. Nini (8) 
where 
“4 
Di(p) = 22 Bmp™. (9) 


m=0 


Since there is no driving function or feedback over 
this interval, the response is the decay characteristic of 
G, with the initial conditions being, again, the terminal 
conditions of the response over the previous interval, 
fon—1. Thus, we can establish recursion relations over 
consecutive intervals. 

Let us first give interpretation to the operational 
symbolism used. Consider the following, where the van 
der Pol-Bremmer transform pair notation is used.? 


2 B. van der Pol and H. Bremmer, “Operational Calculus Based 
on the Two-Sided Laplace Integral,” Cambridge University Press, 
Cambridge, Eng., 2nd ed.; 1955. 
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D2o(p)Ni(p)Fin(Tn + 7’) =~ §(T, +2"). (12) 


Ni(p)N2(p) + Di(p) D2(p) 


On using (10) in (8) and evaluating at 7’=T,—T,-1 
—hn-1, it is useful to define the following matrices; here 
{A }.; denotes the i, jth element of the matrix A. 


{ Vi V (qa+m—1—r) 
{Wilt) bar = DD Bnwis(t), om >1+r 
j=1 m=1-+-r’ 
(RixeV ae, es Ones), ie nea, 
ee (q’-+m—1-r') 
{Wa(r)} ae = DD) Ymwvox(r), m>1-+1', 
k=1 m=1 
(VSR) Or = One ee Ra 
i 
F on—1(7) = MF Pear An at ote Lees = T), Fe sie 
V—1) 
fecal Tuan) (Ve 


(1) 


Pi = {foo(0), foo (0), - - 
Faye dgnt. Po 5 ae ee 


W=1) 
Din 


foo (0)},(R X 1) G6) 


(Tr + 7)}, VX 1). 


The JT superscript denotes the transpose, and “(xy)” 
denotes x rows and y columns. Further, it is clear that 


(17) 


fon (EV fue ta 
namely, the initial conditions of the first V—1 deriva- 
tives of fon(7,+7) are the terminal conditions of 
fon(Ln—-1+hn—1 +7). Thus, on using matrices (13) through 
(17) and (18), the first V—1 derivatives of (6), evaluated 
at the terminal points of each interval, may be expressed 
by the matrix equation 


Fon(hn) = W o(hn) W1(Tn = (gia ee aay)? ogee) 
+ Fin(hn). 
Eq. (19) is a recurrence relation between Fo,(h,) and 


Fon(Mn—1). Hence, we may recur back to @o, (16), and 
obtain 


(19) 


n 


Fon(An) = i {9% oon) (AT n ae! lin’ 1)}°W o(ho) 0 


n=l 
+ DS IL [ate UAT a — Barna) } Fenn) (20) 
n’=0 n'’4+1 
where 
AT, = HES 4 (21) 


n’'—1 
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and [[ denotes the ordered product; i.e., for n>n’', 


matrices Gr, 


LL GusiGeceaiece 


[| is defined to be unity in the following sense: 
n+1 


(LGC. (22) 


n+1 
The first term on the right of (20) denotes the response 
due to initial conditions at t=0, the second denotes 
that due to the driving function f;. 

Thus, the response at the terminal point of each 
pulsing interval, T,,+h,, has been obtained. The output 
behavior during the interval is given by a substitution 
of (20) into 


Fon(7) = ae (Welt) Wi(AT, — Tin—1) F on—1(Mn) + F;, (7) 


Ota (23) 


The response between pulsing intervals is given by (8), 


with (10) and (20) substituted. 


UNIFORM PULSING 


The general solution (20) takes on a less complicated 
form for the uniformly pulsed case. Here 


Tle WEE 3 hy eh (24) 
Then we have 
Fon(h) = {Wo(h)Ws(T — h)}"Wo(h) oo 
+3 (waywir — Ay” FB). 25) 


n’=0 


or, to express the weighting function for the response at 
the end of each pulse interval, 


Fan(t) = Dy Ws” [Fuh + 07 -(0) 60) 
where 
Ws = We(h)Wi(T — h) (26) 
and 
6no = 1, n = 0 
= 0, otherwise. 


Hence, W3” is the generalized weighting function at the 
pulse terminal points for the uniformly pulsed case. Eq. 
(26) could have been obtained also by substituting (24) 
into (19) and using the z-transform representation.® A 


3E. I. Jury, “Analysis and synthesis of sampled data control 
systems,” Trans, AIEE, vol. 73, pt. 1, pp. 332-346; September, 1954. 
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substitution of (24) and (26) into (23) and (8) then de- 
termines the continuous response during and between 
pulsing intervals, respectively. 

It is seen that in the solution (26) the matrix 
Wo(h)W1(T —h) must be raised to integral powers. Con- 
sidering (13) and (14), we see that this matrix must be 
square, VXV. Assuming that it is of full rank V 
of distinct eigenvalues (of zero multiplicity), then 
{ 7.(h) W1(T —h) }m, m> V, is reducible by the Cayley- 
Hamilton theorem, to a polynomial of degree not 
greater than V—1. This can be stated in another way. 
Let 


Wo(h)Wi(T — h) = SAS, 
where A is the matrix of eigenvalues of We(h)W1(T —/h) 
and S is the generator of the similarity transformation 


which diagonalizes it. 
Then 


{wWo(h)wi(T — h)}™ = SAS. 
Therefore, in order to raise such a matrix to a power, tl. 
eigenvalues and matrix S are necessary.* 


STABILITY CONSIDERATIONS FOR UNIFORM PULSING 


Let us require the following behavior on the transient 
response. Given (J, ey) such that 


|| Fonl| = max | Fon(h) | <ev, foralla>WN (27) 


where 


max | Fon(h) | 
denotes the largest member of the R-tuple: 


Ch fon. PS | fon CB) | ee fan Gye era 


Let us define the following norm for matrices: 


Al) = max E | {Aha = sup 


where x is a column matrix; 7.e., the norm is the maximal 
row vector. It is evident that the norm in the case of a 
column matrix reduces to that of (27). The consequence 
of the following derivation, from functional analysis,} is 
to be found very useful. Given matrices A, B and col- 
umn vector x. Then from (27), (28) 


||BAal] <||Ball-lal] < [Bll All-llall, 


i al = = a (28) 


or equivalently 


|B.4|| < |]3I| |] ll. (29) 


4B, Friedman, “Principles and Techniques bd ee Mathe- 
matics, Ife Wiley and Sons, Inc., New York, N. Y., 1956. 

. Kolmogorov and S. V, Fomin, “Functional Analysis: Metric 
aay) Normed Spaces,” Graylock Press, Rochester, N. Y., vol. 1; 1959, 
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It is assumed that we are free to modify the feedback 
component Gz in order to attempt to fulfill (27). Since 
the response to an initial condition vector @o1 (16), is 
explicitly a function of We(k), (26), which contains the 
parameters of G2 implicitly, [(7) and (17) ], the problem 
is first to obtain bounds on W satisfying (27), and then 
to derive suitable parameters for Gz. Thus, from (27) 


Fon = (W2W1)"Wodo. (30) 
Then, according to definition (28) and result (29) 
|| Foal] = |] Cvs)¥ Wado] < |lwel|”#*-||4l]¥-[]G0l] (31) 
Thus, if 
1/(N+1) 
sal] < (a7) (32) 
© ~ \Ifsvill* lel] 
it follows that 
|| Fox|| < ea. (33) 


But now we must also have boundedness for all n> JN. 
Consider the original recurrence relation with zero 
“driving” function. 


Fo, = Wh) WiA(T — h)Fo,. (34) 
Thus 
|Foxaall < || v2)||||-w1(T — H)||l| Foul]. (35) 
Hence, if 
| ey (N41) 
wl < min) req) 69 
then (32) will be fulfilled and also 
| ||Fonsil| < ||¥o,]| for alla > N, (37) 
all of which imply 
|| Fo,|| < ew for all p> WN. 


It also follows that if only the inequality in (36) is ful- 
filled, the vector sequence { Fo,} converges. to zero 
faster than some geometric sequence, 


{om}, |e] <1. 


EXAMPLE A 


In the example we shall consider the response to a 
unit step of a second order system in which g; is of 
second order and gs is a gain constant. Although calcu- 
lational details will be carried out to an extent, they are 
to be secondary to the derivation of the general form of 
the solution. 

Let us first consider the general form of 22 matrix 


raised to a power 
ze) 
Wee 


WW; = ee (38) 


Wai 
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The eigenvalues of 2 are determined by 
ina aN Wie 


det. (3; — XZ) = | =0. (39) 


Wa1 W359 d 


Thus 


Wi1 + We (wii == We)” 
Al = eras any a a — W11W22 + Wi2We1, 


(wir + Wie)? 
A2=A1 — 2 fos — Wi1W22 + Wy2W21. (40) 


Now, to determine the generator S of the similarity 
transformation, 


i 0 
Se Sa = ie ) (41) 
Oe Xe 
and 
Ce = 7S STAM Ss (42) 
So the following must be solved: 
w= AAS (43) 
or 
(*) & a) 
Q — 
S12 Ar Siz 
and 
o(<") * @ * (44) 
S20 deo Soe 
which yields as a possible solution 
—WepS —WerS 
He 3 2) a5 Sot ES 210 22 : (45) 
Wii — AG Wir — A2 
Therefore, we may set Si:=S22=1 and obtain 
— Wai Wa 
Siu = ’ Soa = (46) 
% Wii — At . Wi1 — Ao 
The inverse of S, S?:SS7=S4S=J 
aie a os 
O21 O22 
1 (wir — Ar) (wir — Deo) 
CU ty hcg aaah een ne ag ge a? Ene ere eee eg aa OM I) 
Wie Xe Ni 
(w11 — Az) 
(pl) Se Wont SS Oy) (47) 
Thus 
Devs iS Ans 


ate + ordo™S21, 71141512 + ey (48) 
oo1AS11 + 22201, 7211S 12 + T22d25 22 


We shall keep this result for later reference. 


22 IRE TRANSACTIONS ON AUTOMATIC CONTROL 


As already mentioned, G; is a second order system and 
’ a 
Gz is a constant gain factor. Hence 


N,= Ap+B, (49) 
D, = (p + a1)(p + a2) = p? + (a1 + a2)p + aie 
2 
= 2) Bmp” (50) 
m=0 
J — ! { —ayt —aat} pas i (54) 
Decik 8S) Soa a8 4) : aks ei 
1 
No» eS ils D2= K, G.=—») (52) 
K 
NiN»2 + D,Dz2 
= K( + ( | \pebenay ae 
= pe (a4 = a) + =)? a1a2 =) 
2 
= 12) ¥nd™ (53) 
m=0 
1 
NiN2+ DiD2 


-(004/ Yoro 2) —(omr By 


| e(artasta/K+V )t — g—(artarta/K—v" yt 


== 93()) = Kole ott — C2). 


(54) 
Thus, from (13) and (17) 


Biri(t) + Borr'(t), Borr(t) ) (55) 
Birr’ (t) +Ber1'"(t),  Birr'(t) + Borr'(d)7’ : 
y102(t) + yove"(4), Y2v2(t) 


2 SEs) . 5 
cy) Co + yave(t), —-yite(t) + a oe 


Wilt) = ( 


Let the input f; be a unit step 


1 

= HO) (57) 

fi(ni + 7) = HG — nT) — A(t — nT +h). (58) 
Thus, from (12) 


DN, 


F,(nt = ——_——_—_-f,, (nT 59 
(nt + 1) apLoON ee at) (59) 
and from (49) and (54), 
DN, ( d ) 
——_——— =(|KA—+8B 60 
WiN2 + DiD: Pr ase kG) 
nit+r d 
~ F(nt + T) = if (xa — + B) vo(T) dr. 
nT? dr 
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F(nT +7) = KAKo\ (e-" — 1)eCint 4 (Car — 1)e-CanT} 


1 
= BRKs \= Cas Le ieee 
G 


1 


1 
tale = neem, (61 


or 


Fi (nT + T) =S fooe C1"? + fore CO? 


I 


Fi (nT + 7) = five CO"? + frre C2"? . (62) 


Zero initial conditions are to be assumed for this case. 
Then (26) becomes 


Fon(h) a » W 3" F (kh), 


n=0 
where 


Hence the output and first derivative of the output take 
the form 


™m 


foals) = DS { {we} 1 Fi(nT + h) 


n=0 


+ {97 9"—"} oF! (nT + h)} 


™ 


fon (h) = » { {we} os Fi(nT + h) 


+ {awn} ook y' (nT a h) } (64) 


Disregarding calculational details, usage of (48), and 
(62) yield the following form of the output vector 
Fon(h) wis@iciee 


™m™ 


fon (it) = > {rim { Fa%e-Can? + 12% C27 


0 


+ dam—nf F50¢-Cinr + S22% C27} } 


™m 


fon'(h) = Do {rmmnt Sule C1"? $F y9/e Cam? 


0 


+ rom 2f Foy’eCumP Soo"e Caml (65) 
where Ai, Az are the eigenvalues of W3 = (We(k) Wi(t—h)), 
(Ai) —(As3), and the ¢;;°, £;,’ are obtained from (48) and 
(62) as products of elements of the S, S-! matrices (41) 
and (47) and the coefficient f;; of (62). Hence, the vector 
Fon(h) is available in closed form since the geometric 
progressions indicated are easily summable. 
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2 d;” + A 
fom’(h) = ay = ’ 
om (I) x ON Pa cM wile (66) 
‘ \y™ se Ay le Cum Tr | ‘ Ay” ees Ay te C2mT No” 2M Ao te Cit do” <a eg Le Com 
fom(h) Fu iene r £12 po aowr So an) roa Ve) F £29 eakeeeys 
Fom h) — = 
. . i” = Arte Cum T : \i™ eee Ay te CwmT do” a Ng te CimT do” — hg le-CamP 
Fom’ (i) fu r $12 + £21 + £20" 


1 — Ay te Cir 


EXAMPLE B 


Here, the system described in Example A will be 
stabilized. It will be assumed that the forward loop com- 
ponent G; is inalterable; hence, we must obtain a proper 
bounding value for G;=1/K in order to have the tran- 
sient response behave according to prescription. Given 
the initial condition vector, 


Foo 
go = | Dale 
Foo 
it is required that the transient response to ¢» decay to 


ey or less on and after an output time NT. Thus, ac- 
cording to (35) 


1 — Ar te C2? 


1 — Ag te C1? 1 <a No te Camr 


and from (50) 


Cie OCs aes, ee paeeh L, (69) 
and on referring to (55) and where 7—h=7, 
Wi(T — h) 
di Ce } Le-0.7 2(¢0.35 == a) 
(Cie) ace 28 a (e0.8) 
0.45, 0.62 
% (ea el 70) 


It is seen that the norm of W,(7—hA), as defined by 
(28), is 


|| 3% 2(A)|| |wi(T — h)|| = 0.93 + 0.45 = 1.38. (71) 
< min (( eae jes <=). Further, from (67), 
Ie ale) 7 Ree =a) a eae ae 
The following numerical values for the parameters con- 
Hence, 
cerned are to be used... 
foo = 0, eyed 001, |sv2(H)| ek 
1 
foo = 1, N= 110; < min (( oy ) ‘ 1 ) 
ene. T = 10 Iwi(r — 1)]|*|| 0 I(T — || 
? ? 1 1 
lad BED; = min (= —— (0.00138)*"). (73) 
1Se adt3s 
ie ORS: 
a. = 1, (67) Now (0.00138)"!1=0.55. Thus 
Thus, from (51) and (67) I|972(%)|| < 0.39. (74) 
vy(t) = 2{ e-5¢ — e~*}, (68) Referring to (54) and (67), 
( A 2 =] 
a1 + a2 + =) , ile 
K B ! re 
K, = |2K ~—__*f = (esa: + =) = ( -2k +1) ; (75) 
4 K 
2 
ies oD Yas i =) 
A ( B 3 2 1 2 
ODS ie omar 4 ~ (ae + 2-54 2+ oe (76) 
/ =) 
(« + a2+ =) B { 5 1 ign 
PRC oe eer ee 4 ACCS Se a. 
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where A substitution of values (67) referred to (75), (76), 
vo(t) = Ko(e-Cut — ¢-02!), (78) and (77) yields the following data for K. 

K M 
It is also seen from (53) that 
Gat tr 1/10 0.13 
i1= K(a+ a+) => 4 ay 1/9.5 0.26 
Sos te 1/8.5 0.39 
yo = K. (79) 1/8 0.42 


Thus, a suitable bound for || 9%72(2)|| is obtained. By re- 
ferring to (56), 


I|2(2)|| <M | (1 z= Cry2) Kee~@1" 
Se yao 41) Ke~°2* | + max 
; { | y2Ko(e—O1" — e— Cah le 


| (Cry2 — Ciy1) Kee" (80) 
+ (11C2 — y2C2?) Koe~©2" | } 
and (74) is certainly satisfied if 
M < 0.15. (81) 


Therefore, it appears that the transient response is defi- 
nitely bounded in absolute value less than 0.001 for 
positive K less than or equal to 20. 


CONCLUSION 


It has been shown that a generalized weighting func- 
tion can be developed for a linear pulse-modulated error 
feedback system which permits a formal solution in 
closed form as a matrix operator acting on an appropri- 
ately defined input vector. Two examples have demon- 
strated the mathematical procedures involved. 


Optimization Based on a Square-Error Criterion with 
an Arbitrary Weighting Function’ 
G. J. MURPHY} anv N. T. BOLD 


Summary—Several important criteria for the performance of 
communication systems and control systems are reviewed, and a 
new criterion (the mean-weighted-square-error criterion) is then intro- 
troduced. This is shown to be a special form of a very general 
criterion proposed earlier, but to have special significance in that it 
is a generalization of the familiar mean-square-error criterion. 

The minimization of the mean-weighted-square error is treated 
in detail, and a solution for the optimum physically realizable fre- 
quency function of the system is given. 


INTRODUCTION 


N recent years there have been many attempts to 

establish practical performance criteria for com- 

munication systems, but in most cases the measure 
of system performance has been limited to a specific 
function of the system error. 

Usually, the first step in designing a system for opti- 
mum performance is to attribute “importance” to the 
error as a function of the magnitude of the error, and 
the next step is to minimize the appropriate function of 


* Manuscript received by the PGAC, March 2, 1959; revised 
manuscript received, August 17, 1959. 
ft Electrical Engrg. Dept., Northwestern University, Evanston, III. 


the error. To provide a background for the discussion 
of a square-error criterion with an arbitrary weighting 
function, several well-known criteria are reviewed below, 
and the concept of “importance of error” is discussed 
briefly. 

The work done in the United States on the statistical 
anathesis (the word “anathesis” is defined herein to 
mean “analysis and synthesis”) of systems has its origin 
in the work [1] done by Norbert Wiener in 1942 for the 
United States government. Wiener’s objective was to 
develop a technique for the analytical design of an opti- 
mum linear system, and he established a performance 
criterion based on the (nondeterministic) input signals 
received by the system under normal operating condi- 
tions. This criterion, referred to as the mean-square- 
error criterion, is characterized by 


E; = e(8) (1) 


where 


e(1) 


ca(t) — cal?) (2) 
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is the system error, ca(t) is the desired response of the 
system, and c¢,(¢) is its actual response. That is, 
1 ft 
Ey = lim — 
T->0 21 =f 


[ca(t) — ca(t)]?de. (3) 
For deterministic or transient input signals, with 
et) =0, #<0, (4) 


Hall [2] has proposed the use of 


i= | ” 621) di (5) 


as a measure of system performance. 

For application in situations in which (4) is satisfied 
Graham and Lathrop [3] have proposed the use of the 
“integral of the time-multiplied absolute error” (ITAE) 
criterion, 


A 


E; = | “t| e(t) | at, (6) 


and Nims [4] has suggested the use of the “integral of 
the time-multiplied error” (ITE) criterion, 


E,(t) = f rede (7) 


Other integrals [3] that have been used to measure 
system performance include 


Ee= if “| es) | a, (8) 
Es = i “A| eo| ds, (9) 
ae f ” te2(t)dl, (10) 
Pete P * 2e(i)dt (11) 


and others similar to £3 but differing from it in the 
powers of ¢ and e(t) employed. 

A more general criterion, suggested by Schultz and 
Rideout [5], is 

E,(r) a ip W(t — r)Fle(t, 7) |d1, (12) 
0 

_ where W(é) is an arbitrary weighting function and the 
error is made to depend on a parameter 7 by letting 


ale b=), (13) 


where r(t) is the input to the system. Schultz and Ride- 
out have investigated the use of the special form 


y= J mf 203) | a (14) 
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no 
or 


and Spooner and Rideout [6] have studied the use of 


4 


[e(t, 7) |2dd 


ay: 1 
Ex,(7) = lim —s 
To 2 —T 


(15) 


which is also a special form of (12). 

The integrand in each of the integrals in (3) and 
(5)—(11) may be regarded as the importance of the error. 
For example, the importance of error implied by the use 
of a square-error criterion is shown in Fig. 1, where it is 
evident that errors greater than unity are emphasized 
while errors less than unity are de-emphasized. 

As a second example, the importance of error associ- 
ated with the IAE criterion is shown in Fig. 2. In this 
case, the importance of the error is directly proportional 
to the magnitude of the error. 

Another important function, suggested by Truxal 
[7], is shown in Fig. 3. In this case proportional emphasis 
is placed on error magnitude for errors less than é9, and 
all errors of magnitude greater than é are regarded as 
being equally undesirable. 


THE MEAN-WEIGHTED-SQUARE-ERROR CRITERION 


In all the criteria listed above, with the exception of 
E,(r), some particular weighting of the error is specified. 
The advantage of /9(r) over the other criteria is, there- 
fore, very great; it is considerably more general and, in 
fact, it encompasses the other criteria. However, be- 
cause of its generality, E9(r) is rather difficult to employ 
to advantage. 

A slightly less general criterion, which is readily ob- 
tained as a special form of E,(r) by letting F[e(t, r) | be 
e*(t) and setting 7 equal to zero, is 


Ey = E,(0) (16) 


= f woewa, (14 


0 
which, from now on, will be referred to as the “integral 
of the weighted square error” (IWSE) criterion. 

The corresponding criterion for optimum anathesis of 
linear systems with random inputs, instead of determi- 
nistic inputs, is 

A 1 4 | 
Ex3 = lim — W (t)e2(t)dt, (18) 
T © oh =r 

which will be referred to, from now on, as the “mean- 
weighted-square-error” criterion. Furthermore, these 
two criteria encompass all the others listed above except 
E,(t) and, of course, E,(r), if W(t) is allowed to be a 
function of e(¢). However, for the purposes of this paper, 
W(t) is understood to be a function of ¢ alone. 

For convenience of notation, let £13 be denoted by £. 
Then the mean-weighted-square-error criterion becomes 


1 r 
E = lim — W (t)e?(t)dt 


TH © 7 


(19) 


no 
(op) 


e(t) 


e(t) 


Fig. 1—Error importance function for the error associated 
with the square-error criterion, e7(¢). 


le(*)| 


e(t) 


Fig. 2—Error importance function for the error associated 
with the absolute error criterion |e(¢)|. 
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Fig. 3—Saturating error importance function. 


and the associated problem is that of determining the 
optimum system transfer function Gopi(s), for which E 
has its minimum value for a given W(t). 

If the transfer function of the system is denoted by 
G(s), then, by application of the real convolution 
theorem, 


call) = [gle rl = adr, (20) 


where g(t) is the inverse Laplace transform of G(s). Since 


IRE TRANSACTIONS ON AUTOMATIC CONTROL 


January 


e(t) = calt) — call), (21) 


substitution of the right-hand member of (20) into (19) 
yields 


1 
E = lim — “Ww [ca2(t) — rea) g(r)r(t — 7)dr 
Too 2T 2G) 
+f eore=1 [sere —verdr|ar 22) 
li : W (t)di 
nee Lr i. 
(a) 1 She 
aS tis g(r) lim = Bb onli 
T 
T ) lim —— 
+f of x 7) tim > Dsk 
-W(i)r(t — r)r(t — y)dtdydr. (23)! 
Now let 
A 1 e. 
Gime weal Ts 7) = lim. W (t)ca(t ate T)ca(t = y)dt, (24) 
Ps Dd an 
i 1 
Dace T, y) = fim = [e W(t)cali + r)r(t + y)dt, (25) 
Too 2T == yt 
and 
Peer Ene ty =f" Writ rr(t+ dl. (26) 
T> 0 
If 
W(t) = 1, —x<ti<cw (27) 
then 
@rneaea\T, Y) ae PeacalY ae T (28) 
Cy a Gre 7) = Pearl¥ = T)s (29) 
and 
Pwrr(T, 7) = br a T), (30) 


where ¢egcg, Pear, and ¢,, are the well known auto- and 
cross-correlation functions of ca(t) and r(é). Since for 
arbitrary W(t), Pwcaegr Pwegr) ANd Purr are functions of the 
two real variables 7 and y, instead of the one real variable 
y—rt. As for the W(é) defined in (27), the functions de- 
fined in (24)—(26) will be referred to as “correlation 
functions in two-space.” 


1 The inversion of the order in which the indicated operations are 
performed can be justified if the system is stable and W(t), r(#), and 
the functions defined in (24)-(26) are bounded. 


, err 
a 
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Substitution of the left-hand members of (24)—(26) 
into (23) yields 


Foils, (0-0).— 2 f Peen(O, ar 


—oo 


+f 0) [ sdbunl—7, —v)avdr (31) 
= E=zc 7) = J s0dbecals 1 x) dx 
wi f £(X) bwearl¥, Y a x)dx 
+f g(a) f £(9) bwrr(t — my — yaya | A o2) 


where it is evident that the mean-weighted-square error 
is completely determined when the unit-impulse re- 
sponse of the system and the two-space correlation func- 
tions ducar(T, Y), Pwegea(T, Y), aNd hwrr(7, Y) are specified. 


MINIMIZATION OF THE MEAN-WEIGHTED- 
SQUARE ERROR 


The mean-weighted-square error can be minimized by 
application of variational calculus as follows. Let &(x) 
be a continuous function of x such that 


E(o) = 0 (33a) 


and 


E(x) = 0, (33b) 


and let it be assumed that g(x) is continuous for x>0. 
If g(x) +(x) is substituted for g(x) and g(y) + B&(y) is 
substituted for g(y) in the right-hand member of (32), 
then the left-hand member of (32) becomes E+AE. If 
(32) is subtracted from the equation thus obtained, the 
result is the variation in E which is thus found to be 


a Ke 


ie: | -2 J ietecat y = a)dx 
8 f EG) bec 1 x)dx 
+8 fi J on ~ 2,4 — y)dydx 
ee 
+e f 


(34) 


[dbl = 7 - s)dyd | 


T=y=0- 
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The necessary condition for an extremum of AE is 


JAE 


— == (). (35) 
OB |g=0 
From (34) it is found that 


OAE 
dp 


B=0 


== {- [ © [Swear(7; Dinan «) = PwearlY, T om x) |dx 


+ { g(x) | Eo)bunlr = x1 = ydyde 


+f f gONbunlr — 2,7 = y)dyae (36) 


T™=y—0 


Changing the order of integration? in the first double 
integral on the right-hand side of (36) yields 


OAE 


op 


B=0 


_ ‘S f &@ [bwear(T, Pee «) aie wor; Ube x) |dx 


= f #(x) f PON ae Ceerases 


+ fe) [ s0%enlr—27=)ayart| 
= | fe 1 fs0)lbon(a ee ae) 

=f Purr(T en tS eae y) |dy a Pwear(T; a} aa x) 

i te Disoar Cys oa 2) ax | (38) 


According to (35), the right-hand side of (38) must 
equal zero for arbitrary &(«) satisfying (33a) and (33b). 
It follows that 


{fo [Dwrr(¥ Peake y a ¢ y) +E Pare (T. Saito dens y) |dy 


0, 


ar Pwoar Ts 7) ae x) ae PwearlY, Urs a) 


T=y=0 


D> 0239) 


2 A justification of this procedure is given in the Appendix. 


28 IRE TRANSACTIONS ON AUTOMATIC CONTROL 


Taking the double Fourier transform of both sides of 
(39) gives 


Lf ci [om feo) lbenls = 5,7 — 9) 


={- Pier? Sree yt Vie an y) |dydydr 


ie) eo 
= —0 


+ Oocar(¥, 7 — «)|ayart Le 0 


= 0, (40) 


r=y=0 


if the integral on the left-hand side of (40) is absolutely 
convergent. Changing the order of integration’ then 
gives 


\ J (9) [Yrorr( jor, joo’) Galt ea iee Joe Vwrr( Go’, jw) e 10! 2e—iw u|dy 


= CF ye ar( jer’, Jw) ae “fa je j0!)} 


t=y=0 


“>0, (41) 


where Ww, is a double Fourier transform of ¢.;;, and so 
on. 
Eq. (41) can be rewritten in the form 


[Gj Yoorr jor jo’) it Wwear( Jo’, jw) Je-i#= 
+ [G( jee) Prove ja", jee) — Wrcar( jo, joo") ]e-H"*} 


TtT=y=0 = 0, 


G > Go (42) 


Taking the double Fourier transform of both sides of 
(42) gives 


1 rr) C-} . 
‘Ag? f if [G(jo! Porr(jo, jo’) — Wucar(jo’, jo) |e**dwdu! 


+ = ii y f [66s ani a) 


= Wwcar( jo, jo’) |e-%’*dwdw’ = 0, «> 0. (43) 


It can be shown that the first integral on the left-hand 
side of (43) cannot be the negative of the second integral. 
Therefore, if a solution exists it is obtained when both 
integrals are zero. This will be the case if the integrands 
are identically zero for x >0. Thus, it is found that the 
mean-weighted-square error is minimized if the system 
frequency function G(jw) satisfies 

G(jo" Wrorr( jo, jo’) = Wwear( Jo’, jo) (44) 
and 
G(jo)Pwrr(Jo’, jo) ae Pwear(jo, jo’). (45) 


3 The justification for this operation is Fubini’s theorem. See 
Appendix. 
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The solution to these equations is 
Wwear(Jo, Jo’) 


(46) 
Vrorr( jo’, jw) 


Gopt( Jw) aa 


which serves to define the optimum frequency function 
for minimization of the mean-weighted-square error. 

In a typical situation of considerable interest, the in- 
put to the system is the sum of a signal s(#) and noise 
n(t) which are uncorrelated‘ in the two-dimensional 
space (7, y), and the desired output is 


t 
petit f THE (47) 
In this event, 
Vwear( Jo, jo = Gal(jw)Wwss( Jo, jo") (48) 
and therefore (46) becomes 
; Ga(jo)Wuss(Jw, Jo’) 
Gopt( jw) = au 5 (49) 


Vwss (Jo, jo’) +: Voonn(Jo, jo’) 


PHYSICAL REALIZABILITY 


The existence of the solution indicated in (46) does 
not imply physical realizability of the optimum fre- 
quency function. Accordingly, the problem of designing 
the optimum system on the basis of the weighted-mean- 
square-error criterion is not in general completely solved 
by application of (46). 

The determination of the optimum physically realizable 
frequency function is based on the fact that a transfer 
function G(s) is physically realizable only if 


a) g(t) = 0, be 0 
and 
b) g(¢) is bounded’ for ¢ > 0. 


Since the existence of poles of G(s) in the right-hand half 
of the s plane requires that g(t) differ from zero for t <0, 
it is necessary that no poles of G(s) lie to the right of the 
imaginary axis in the s plane. Therefore, if the optimum 
frequency function given by (46) or (47) is a rational 
function of w with no poles below the real axis of the w 
plane, the optimum frequency function is physically 
realizable. 

If Gops(jw) is a rational function® of w with one or more 
poles below the real axis in the w plane, the following ap- 


4That is; 


dwsn(T, 7) = duns(7, 7) =0. 


5 From a practical point of view, this condition can usually be 
replaced by 


lim g(t) = 0. 


6 A rational G(jw) multiplied by the exponential factor e~* may 
be treated similarly. 
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proach may be used. As in the optimum synthesis pro- 
cedure based on the rms error criterion, let the ordinary 
power spectral density of the system input r(t) be 


®,,(s) = ,,+(s)®,,--(s), (50) 


where ®,,+(s) has poles and zeros only in the left-hand 
half of the s plane and ©,,-(s) has poles and zeros only 
in the right-hand half of the s plane. The input r(¢) is 
first converted to white noise by passing it through a 
unit with the physically realizable transfer function 


Gi(s) = (51) 


&,,+(s) 


The white noise must be passed through a physically 
unrealizable unit with the transfer function 


Gopt() 
Gi(s) 


G2(s) — (52) 


to obtain the optimum response. This optimum response 
is the sum of two statistically independent components 
[7]: 1) a completely predictable component due to all 
pulses of the white noise that have previously occurred 
and the pulse (if one exists) that is simultaneously oc- 
curring, and 2) a completely unpredictable component 
due to pulses of the white noise that will occur in the 
future. 

On the average, the best physically realizable measure 


’ of the second of these components is zero. Therefore, the 


optimum physically realizable frequency function 
Goor(jw) is obtained by making 


(@) = ‘s ee 
20r -~ g-1 [@,,*+(jw) Gopt(jw) ] ) 


By combining (53) and (46) it is then found that the 
optimum physically realizable frequency function 


Gor(jw) is 
jo 
f eft 
—jo 


dwdt. 


$>0,° (53) 


1 eaten 
‘eee I = is gat Z 
(je) 5 FGA Oni 


prihseveolicrdit (54) 
Wore (ja, jo’) 


Similarly, for the special case where the system input 
is the sum of a signal s(t) and noise u(t) which are un- 
correlated in the two-dimensional space, and where the 
desired response is given by (47), it is found that the 
optimum physically realizable frequency function is 


pearson 9p sachs 
Gor 10) = : | Gi" : f é “” 
: ®,,*(jw) Jo 2aj J 


i Ga(jo) Bret (Jo) Wwes( Jo, jo’) 
Vwss(Jo, jw’) 5 Vann (Jo, jo’) 


duwdt, (55) 
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LIMITATIONS 


Although (35) is a statement of a necessary condition 
for minimization of the mean-weighted-square error, it 
is not sufficient. In order that the extremum resulting 
from satisfaction of (35) be a minimum, it is necessary 
in addition that 


AE 
sates (56) 
In general, 
AE 
OB? |pm0 
=2f 6) [ t0)benlr— 5,7 - vdyde] 7) 
. : 1 T ee 
=2f «f © lm 
-W(t)r(—«)r(—y)didydx. (58) 


Changing the order of integration’ and averaging yields 


CAE “ 2 T oe) we) 
i eee) J eon-» ie 
-E(y)r(t — y)dydaxdt (59) 
. 2 “4 2 
= aa ae eds he (t)dt, (60) 
where 
m(t) = i eerie Hane (61) 


Thus, it is seen that the extremum obtained by satisfy- 
ing (35) is a minimum if 


W(t) > 0, ae ee ee ye (62) 


which is acceptable as a sufficient condition that does 
not impose a severe limitation on the choice of W(é). 


CONCLUSION 


A fairly general new criterion for the performance of 
communication and control systems has been presented 
and investigated in detail. A technique for minimizing 
the mean-weighted-square error (MWSE) has been pre- 
sented as has been a solution for the optimum physically 
realizable frequency function. It has been found that the 
procedure for the use of this criterion is very similar to 
that for the use of the familiar mean-square-error cri- 


7 The justification for this operation is similar to that used in ob- 
taining (37) from (36). 
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terion. However, in using the new MWSE criterion, 
the anathesist has some freedom in the manner in which 
he may assign importance or undesirability to errors at 
different instants in time. 


APPENDIX 
If E(y) is chosen so that 


[8@) f e006 = 7 = ydayde (63) 
is absolutely convergent, then by Fubini’s theorem,® the 
repeated integrals taken in any order exist and are equal. 
It has previously been stated that &(y) is continuous and 
arbitrary except for the physical realizability restric- 
tions, 


Ey) = 0, ~~ y 50 (64) 


and 
(co) = 0. (65) 


8 E. W. Hobson, “The Theory of Functions of a Real Variable,” 
Harren Press, Washington, D. C., vol. I, p. 630; 1950. 
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In addition, it is now understood that only &(y) for 
which (63) is absolutely convergent will be permitted. 
It follows that then the order of integration may be 
changed as required. 


BIBLIOGRAPHY 


[1] N. Wiener, “Extrapolation, Interpolation, and Smoothing of Sta- 
tionary Time Series,” The Technology Press, Mass. Inst. Tech., 
Cambridge, Mass.; 1959. 

[2] A. C. Hall, “The Analysis and Synthesis of Linear Servomecha- 
nisms,” The Technology Press, Mass. Inst. Tech., Cambridge, 
Mass.; 1943. 

[3] D. Graham and R. C. Lathrop, “The synthesis of ‘optimum’ 
transient response: criteria and standard forms,” Trans. AITEE, 
vol. 72, pt. 2, p. 273; November, 1953. 

[4] P. T. Nims, “Some design criteria for automatic controls,” Trans. 
AIEE, vol. 70, pt. 1, p 606; 1951. 

[5] W. C. Schultz and V. C. Rideout, “A general criterion for servo 
performance,” Proc. Natl. Electronics Conf., vol. 13, p. 549; 1957. 

[6] M. G. Spooner and V. C. Rideout, “Correlation studies of linear 
and nonlinear systems,” Proc. Natl. Electronics Conf., vol. 12, p. 
321511956. 

[7] J. G. Truxal, “Automatic Feedback Control System Synthesis,” 
McGraw-Hill Book Co., Inc., New York, N. Y., p. 467; 1955. 

[8] N. T. Bold, “Optimization Based on an Error Criterion with an 
Arbitrary Weighting Function,” Ph.D. dissertation, Northwest- 
ern University, Evanston, IIl.; August, 1958. 

[9] H. W. Bode and C. E. Shannon, “A simplified derivation of linear 
least square smoothing and prediction theory,” Proc. IRE, vol. 
38, pp. 417-425; April, 1950. 


Multiple-Rate Sampled-Data Systems” 


LESTER A. GIMPELSONT 


Summary—tThe characterization of multiple-rate sampled-data 
systems by the ordinary z-transform of single-rate systems is 
shown. Single-rate sampling, or impulse modulation, of continuous 
signals is performed by an impulse modulator, M; the sampled, or 
“starred,” function is described by the z-transform. In an analogous 
manner, a submultiple-rate modulator is introduced; its presence in 
a branch allows the passage of every nth pulse, or a train of pulses 
at a submultiple rate; the nomenclature of single-rate systems is 
continued through the performance of submultiple-rate ‘‘starring” 
of discrete signals and discrete filters. Table I permits starred ex- 
pressions to be rewritten as functions of the z-transform in closed 
form. Techniques are shown for the reduction of discrete, and 
mixed continuous and discrete systems via flow graphs, so that, 
after the modulators are removed from the feedback loops, the 


* Manuscript received by the PGAC, March 12, 1959; revised 
manuscript received, July 28, 1959. This paper is based on a thesis 
submitted in partial fulfillment of the requirements of the S.M. de- 
gree in the Dept. of Elec. Engrg. at the Massachusetts Institute of 
Technology, Cambridge, in January 1959, This report was prepared 
at the M.I.T. Instrumentation Laboratory under the auspices of 
DSR Project 52-140, sponsored by the Weapons Guidance Laboratory 
of the Wright Air Development Center, through USAF Contract 
AF 33(616)-3892, Air Force Task Number 50496. 

{ Dept. of Elec. Engrg., Massachusetts Institute of Technology, 
Cambridge, Mass. 


analysis may proceed by standard methods. Representation of single- 
and submultiple-rate modulation in the s- and z-planes is used to 
demonstrate that submultiple-rate modulation of discrete signals is 
analogous to the impulse modulation of continuous signals. 


INTRODUCTION 


N many applications computing elements form the 
if nucleus of the control, calculation and display oper- 

ations. Since the computer is communicating with 
physical components, definite requirements are placed 
upon it which govern its operational speed. Even after 
the most judicious programming, more efficient use of 
the computation capacity of a digital computer may be 
possible by computing less critical loops less often. 

In the wholly discrete control system of Fig. 1(a), 
each element of the block diagram represents a portion 
of the computer program. If the contribution of the 
upper loop is small compared with that of the lower 
loop, why not compute the upper loop less often? To 
accomplish this, the output of H;(z), shown in Fig. 1(b), 
either is maintained at its previous value for several 
interations before recomputation, as in Fig. 1(c), or is 


} 
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t t t 
(b) (c) (d) 


Fig. 1—(a) Block diagram of discrete control system. (b) Output of 
H(z) if computed at every iteration. (c) H;(z) computed less 
often; previous output held. (d) H(z) computed less often; no 
output between computations. 


just computed less often without being “held” between 
computations, as in Fig. 1(d). Now the contribution of 
the upper loop will be at a lower rate than the basic rate 
of the system. Assuming that each block requires the 
same computation time, the determination of the output 
of H3(z) at every fifth iteration reduces the average 
iteration time by almost 20 per cent. 

This is but one of a large class of digital systems em- 
ploying several pulse rates either for economy or be- 
cause of their basic nature (e.g., radar systems).! This 
paper outlines techniques for the analysis of these sys- 
tems; in addition, their relation to single-rate systems 
is shown. 

A multiple-rate sampled-data system will be charac- 
terized by the presence of synchronous pulses whose 
rates are harmonically related. The synchronous and 
harmonic qualifications require that pulses occur at only 
the basic clock instants and at rates which are sub- 
multiples of the basic rate; as will be shown, there are 
frequently approximations which permit these restric- 
tions to be relaxed. 


SAMPLING OPERATIONS 


Sampling is accomplished by an impulse modulator, 
M [Fig. 2(a)], which transmits pulses (approximated 
by impulses) whose heights are equal to the instantane- 
ous magnitude of the function being sampled. The out- 
put of the modulator, f*(#) [“f(é) starred”], is obtained 
through the multiplication of the input signal by a train 
of unity area impulses separated by T seconds. When 


1G. M. Kranc, “Input-output analysis of multirate eed 
-I, Pp. 
21-28; November, 1957. Also, “Multi-Rate Sampled Systems,” 
Columbia University, New York, N. Y., Rept. No. T-14/B; May, 
1958. ¢ 

G. M. Kranc, “Additional techniques for sampled-data feedback 
problems,” 1957 WESCON ConveEnTION REcorD, pt. 4, pp. 157-165. 
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Fig. 2—Sampling devices. (a) Impulse modulator, M, operating at 
basic system rate. (b) Submultiple-rate modulator, Mijn, sam- 
pling impulses. 


considered in the frequency domain, the impulse 
starring of F(w) produces sidebands? at intervals of w,, 
the sampling rate (Fig. 3). In the s-plane F*(s) is char- 
acterized by an infinite string of poles for each pole 
of F(s). For example, when the exponential function 
shown in Fig. 4(a), 


1 
Fey = 
( S+ a 
is sampled, the single pole at s= —a, is replaced by an 


infinite number of poles [Fig. 4(b) | located at s=njw—a: 


n=00 1 


il 
F*(s)=— > 


be ede Hh) a wah 


These observations follow from the two equivalent fre- 
quency domain expressions for f*(), 


Fes) = fla Tyee”, 
or, 
1 n=00 
FH) = — XS Fs + na), 


where T is the sampling period.’ The difficulty, which 


2W. K. Linvill, “Sampled-data control systems studied through 
comparison of supine with amplitude modulation,” Trans. AITEE, 
vol. 70, pt. 2, pp. 1779-1788; 1951. 

3B. Widrow, class notes for “Pulsed-Data Systems,” Course No. 
6.54, Mass. Inst. Tech., Cambridge. 

J. G. Truxal, “Automatic Feedback Control System Synthesis,” 
McGraw-Hill Book Co., Inc., New York, N. Y., pp. 500-558; 1955. 

4 If f(t) is nonzero at £=0*, the right side of the second equation 
for F*(s) should be appended by the term 3f(0*); in most systems, 
however, this addition is unnecessary. 
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Fig. 4—Frequency domain view of sampling. (a) Single pole of F(jw). 
(b) Pole plot of F*(jw) showing multiplicity of poles. (c) F*(z) in 
z-plane. 


results from infinite numbers of poles, is overcome by 
the substitution, 


gus ef, 


which allows F*(s) to be written as 


n=co 


F(z) = DO f(nT)z. 


For positive time only, : 
Perils fre yiSsrje +, 


which is a description of the sampled function at only 
the sampled instants. (The z-transform, as defined 
above, differs from the notation used by some authors. 
The use of a negative exponent permits z to be consid- 
ered as a unit of delay; expansions are then written in 
powers of z, which is very convenient in the following 
material. Expressions using the alternate definition, 
z=etsT, are the same as those obtained here with every 
z replaced by 271.) 
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The z-transform maps the left half of the s-plane onto 
the area outside® of the unit circle in the z-plane; the 
unit circle corresponds to the jw-axis; in general, circles 
concentric about the origin of the z-plane correspond to 
lines of constant 0, where o is the parameter of the hori- 
zontal axis of the s-plane. Thus, system characteristics 
may be analyzed in the z-plane with techniques similar 
to those employed in the s-plane. 

Representation of a signal occurring at a rate which 
is (1/7)th of the basic rate is possible using powers of z. 
For example, 


= 2[(b2)° + (62)! + (2)? +--+ J, 


1 — bz 
becomes 
br—ign 


page Oe Oe Ca ibe 


when the original function is resampled at every nth 
instant, or at a submultiple rate. 

Resampling is performed by a submultiple-rate modu- 
lator [Fig. 2(b)|], which multiplies the input signal by 
a train of finite width, unity height pulses, separated by 
nT seconds.®:? The output of the submultiple-rate modu- 
lator is termed F*/ to indicate the reduced sampling 
rate. The expressions above can be written as 


k vA *1/n br—lign 
=| ~ 1— (bz) 


Thus, there are two types of starring operations: im- 
pulse and submultiple-rate. Tables of z-transforms® are 
useful in the evaluation of impulsively starred functions; 
determination of submultiply starred functions® is sim- 
plified by the use of Table I for which starred functions 
are expanded into partial fractions of the form 


aml 


where k is the number of units of delay; d governs the 
multiplicity of the poles; 1/d is the location of the poles 
(possibly complex); and 1/n indicates the submultiple- 
rate sampling. (The development of this table, instruc- 


5 The convention employed here defines z = e7*7, 

6 The use of a multiple-rate modulator and the “family” of modu- 
lation processes are based upon suggestions by Prof. B. Widrow, 
Mass. Inst. Tech., Cambridge. 

7 More rapid samplers are inserted in single-rate networks by sev- 
eral authors in order to obtain information concerning the operation 
of these systems between the sampling instants: 

W. K. Linvill and R. M. Sittler, “Extension of conventional 
techniques to the design of sampled-data systems,” 1953 WESCON 
CONVENTION RECORD, pt. 1; 1953. 

G. V. Lago and J. G. Truxal, “The design of sampled-data feed- 
He: systems,” Trans. AIEE, vol. 70, pt. 2, pp. 247-253; November, 
1954 


8 R. H. Baker, “The pulse transfer function and its application to 
sampling servo systems,” Proc. IEE (London), vol. 99, pt. 4, pp. 302- 
317; December, 1952. 

E, I. Jury, “Analysis and synthesis of sampled-data control sys- 
tems,” Trans. AIEE, vol. 73, pt. 1, pp. 332-346; September, 1954. 

E. I. Jury, “Additions to the modified z-transform method,” 1957 
WESCON ConveEnTION RECORD, pt. 4, pp. 136-156. 
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TABLE I 
SUBMULTIPLE-RATE EQUIVALENTS® 
ee ban-kgdn 
[1 — (bz)n]2 
ssi [n —(k — 1) ]b@-Ynmkg (a-1)n [k—- (¢@— 1) |ben-kgan 
[1 — (bs)*]# 
eS pedeigteain + ( Sd ye ) Bee Dete(tvm + Do bdnrkgdn 
ne n—(k-1) k-+n—(d—1) k—(d-1) k—(d-1) 
[1 — (bz)"]4 
> >: bh (4-8)n—-kz (d-8)n ( 2. >» =a De >) p (d-2)n=Ky (d—-2)n ( *3 Ss i) ss >) (4-1) n—ke (d-1)m . So banked 
BRE, n—(k—1) 2n—(k—1) n—(k=1) k4n—(d—1) k—(d—1) k—(d—1) 
[1 — (bz)"]2 


tions for its use and examples of equivalents with real 
and with complex poles are included in the Appendix.) 


NATURE OF MULTIPLE-RATE SYSTEMS 


The use of this new modulator and Table I will be 
illustrated through a simple example. The system in 
Fig. 5 consists of two delay loops; in the flow graph solid 
lines indicate pulses at the basic rate, while the dotted 
branch experiences pulses at the submultiple rate only, 
governed by the closure rate of the switch. The trans- 
mission of the system is 1/(1—z), an integrator, for 
closure of the switch at each basic clock instant. If the 
switch were to close at every third basic clock instant, 
the pulse rate at node a would be one-third of the basic 
rate. By considering the transmission of signal through 
the network, the response of the system to a pulse co- 
incident with a switch closure, and a closure rate of one- 
third, is readily obtained as’ 


To = 1+ 32+ 32? + 22° + g2¢ t+ eae t+ tee! 
a. aa2" — on 
Arrival of the input pulse at one iteration after a switch 
closure leads to 


T, = 1+ 42 + 42? + 428 + fst + Z2° + 2° 


ree ee oe 
Also, 
To = 12a $2? + 428 + 3st + G8 + tee? 
Ape ge Rae 


In closed form, the impulse responses are, 


1 + 92+ 42" 


a 1,—. 42° 
1+ 22+ 22° 

of oy gi nas a 
1 — 42 
1+2+ 32° 

te 
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Fig. 5—Multiple-rate system. 


These responses show three important characteristics of 
multiple-rate sampled systems: they have time-varying 
impulse responses, which are a consequence of the sys- 
tem’s cyclic changes in structure; the system’s degree 
is increased (above responses are third-order, while the 
single-rate system had a first-order response); the poles 
of each response are always the same, but the numera- 
tors are not unique, indicating that while the stability 
and general character of the system is not a function of 
the timing of an input, the exact nature of the response 
is time-varying. 
As Table I shows, 


gk Alln ( Ge gee ) 
ie ee, a [1 — (b2)"]2 
the denominators of all submultiply starred functions 
undergo the same transformations, only dependent upon 
the values of 7 and d. Prior to starring, poles of fractions 
of the form z/(1+6z)? are located at 1/b, and their mul- 
tiplicity is equal to the power, d, to which the denomina- 
tor is raised. After starring, the multiplicity is still 
governed exclusively by d, however; the locations be- 
come a function of m. There are now 7 poles, equally 
spaced about a circle of radius 1/b; see Fig. 6. (Since 
points which are equidistant from the origin of the 
z-plane, lie on lines of constant o in the s-plane, sub- 
multiple-rate sampling cannot affect a system’s settling 
characteristics; rise time and overshoot are, however, a 
function of the lower sampling rate.) If d were increased 
to 2, each pole would be double order; if d=3, triple 
order, etc. Clearly this is so since the modulator cannot 
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Im(Z) 


(d) 
Fig. 6—Zeros of [1 —(42)”]. 


change the fundamental form of the signal it passes, and 
therefore the multiplicity cannot be governed by n. 

Impulse modulation causes each pole of the continu- 
ous function to be repeated an infinite number of times 
in the s-plane [Fig. 4(b) ]. The infinite number of initial 
conditions which would seem to be required is explained 
by the infinite number of ways of phasing the impulse 
modulator with respect to the continuous signal. The 
application of the z-transform coalesces this multitude 
of singularities into a single pole in the z-plane for each 
infinite string of poles in the s-plane [Fig. 4(c) |. When 
a sampled function is again sampled, as by a sub- 
multiple-rate modulator, the resulting function will have 
a higher-order representation in the z-plane (Fig. 6) and 
a multiply infinite number of poles in the s-plane. The 
effect of submultiple-rate sampling is therefore very 
similar to the single-rate sampling of a continuous sig- 
nal, since the occurrence of new poles in the g-plane has 
the same significance as the development of infinite 
numbers of poles in the s-plane. Due to the qualification 
here that the sampling rates be submultiples of the basic 
rate, the number of possible phasings of the modulators 
with respect to external discrete signals is limited, and 
therefore the number of additional poles is finite in the 
z-plane. These poles distribute evenly about circles 
whose radii are those of the original poles. (If sampling 
rates are not commensurate, there will be a continuous 
distribution of poles about these circles.) 

To illustrate these points, consider an integrator 
which has the transforms® 1/s and 1/(1—3). If sub- 
multiply sampled at (1/7)th of the basic rate, the trans- 
form becomes 1/(1—2”). This is still a first-order func- 
tion; however, the sampling is less often. A system in- 
corporating this element together with its modulator has 
n impulse responses, or alternately, there are ” possible 
phasings of the submultiple-rate modulator. If there are 
storage elements in a single-rate system, the degree of 
the system will be equal to their number (e.g., 1/s and 
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1/(1—s), or 1/s? and zT/(1—2z)?). The addition of a sub- 
multiple-rate, 1/n, increases the function’s order to n 
times the original order (e.g., 1/s and 1/(1—3"), or 1/s? 
and g"7/(1—z")?). Inclusion of several modulators re- 
sults in increases proportional to the product of the 1; 
designators. 

The broad areas of similarity in the mathematics and 
the interpretation of these operations leads to the de- 
velopment of a hierarchy of modulations: namely, im- 
pulse modulation is to continuous signals as sub- 
multiple-rate modulation is to discrete signals. 


REDUCTION TECHNIQUES 


Returning to the system of Fig. 5, the switch will be 
replaced by a submultiple-rate modulator, passing 
pulses at one-third of the basic rate. The flow graph of 
this completely digital system [Fig. 7(a) ] is reduced by 
the selection of node a as the “discrete node” since pulses 
occur there at the lower rate only. The flow-graph of 
Fig. 7(b) is obtained by noting 1) the transmission with 
the branch containing Mj1,;3; removed, 1/(1—42); 2) the 
transmission from the input to node a, [$2/(1—4z) |M1s3; 
3) the transmission out of node a, 1/(1—4z); and 4) the 
self-loop about that node, [42/(1—4s) | M13. 

The presence of modulators before and after a con- 
tinuous element permits impulse starring of the element 
and the elimination of the second modulator. In the 
same manner, a discrete filter placed between two sub- 
multiple-rate modulators is submultiply starred to- 
gether with the removal of either modulator. The func- 
tion in the self-loop of Fig. 7(b) meets this criterion, and 
starring is therefore permissible, thus removing the 
modulator from the feedback loop. The transmission 
is now 


22 1 1 
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1) 132 a 52 
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Evaluation of the starred function leads to 


1 $3 1 — $23 1 
—+| : |atus( ———. 
Lis ssee) SLL iog2 1— 428 1 — 2, 


Note that the ordering of the terms in this expression 
carries the same significance as in single-rate sampled- 
data systems. The first term indicates the signal which 
feeds through the network irrespective of the modu- 
lators; the input signal is filtered by the term in square 
brackets and its output is sampled at every third itera- 
tion; finally, the one-third rate signal is filtered by the 
term in parentheses. (At this point it is evident that the 
system is absolutely stable, as every element in the for- 
ward branches is stable, and the modulators cannot, of 
themselves, introduce an instability into a forward 
path.) 
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(c) 
Fig. 7—Reduction of Fig. 4. 


Since any representative’ impulse response will pro- 
vide the system’s pole configuration (proof of this point 
will follow), it is only necessary to determine one con- 
venient response, usually 7». When the input is a single 
pulse at ¢=0, the system is unaffected by the placement 
of a modulator in the input branch [Fig. 7(c) ]. This ad- 
dition allows the branch joining the input with node a 
to be starred, and the response is obtained as 

Li8 
Ih == ela ze 


1 — 32 


14+ 92+32 


t= 
i= 


IE | Cole 
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If the input, R(z), were at the lower rate, the output is 
simply 


C(z) = R(z)- Tol), 


at the lower rate. Through a change of viewpoint, net- 
works which include a higher sampling rate are analyzed 
by considering the fast rate as a pseudo-basic rate, and 
the actual basic rate as a pseudo-submultiple rate. 
Reduced flow graphs, as in Fig. 7(b), are complete 
descriptions of the original system’s terminal character- 


9 Degenerate situations sometimes arise. If, for example, Mijn is 
located at the input, 7» will be the only finite transmission. The 
other transmissions have the same denominators, but are considered 
to have zero numerators. 
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istics at the sampling instants. The cyclic, time-varying 
nature of the impulse responses allows investigation of 
any 7, since an input at time f; can be considered as 
contributing to initial conditions for the response 7; at 
some later time, ¢;. Therefore, the poles of these re- 
sponses cannot depend upon the timing of the input. 
Alternately, by removing the modulators from the feed- 
back loops, this method of reduction will yield flow- 
graphs composed only of forward path transmissions, 
T,’s. The total transmission, 7;, is the sum of the 1 
In this summation, the poles of the 7,’s contribute to 
the poles of 7;; none of the numerators of the 7,’s, 
which are the only functions dependent upon the input, 
can affect the denominator of 7;. These two demonstra- 
tions of the nontime-varying character of the poles of 
the impulse responses permit use of any T; (usually To, 
as in the example above) for analysis. In addition, this 
entire analysis technique employs only the ordinary 
z-transform; as a result, all of the methods of analysis 
(root locus, Nyquist diagrams, etc.) previously de- 
veloped for sampled systems, are readily applicable. 

As a further example, the wholly digital system of 
Fig. 1 will be reconsidered; the flow graph, Fig. 8(a), 
shows two integrators in tandem, the first one providing 
the necessary delay for both feedback loops.?° If, for the 
moment, the modulator is shorted, the impulse response 
is 


abz 
2°(1 — ad) + z(abk + ad — 2) +1 


Elimination of the upper loop replaces T, with 


abz 
2? + 2(abk — 2) 21° 


which has both of its poles on the unit circle in the 
z-plane, indicating that any input would lead to an in- 
definitely persistent output. This system is therefore one 
which is naturally oscillatory, receiving corrections 
when the switch in the upper loop closes. 

Reduction about node a leads to Fig. 8(b), where the 
modulator has been removed from the self-loop by sub- 
multiple-rate starring. 7) is then obtained as 
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roljeilsttish 9(abR 2) pedi let (dies ad) eet pead sls 


13 


abz 


if the upper loop sampling rate is one-half of the basic 
sampling rate; K,; and Kz are constants. In all cases 
where 1/n#1, the presence of two poles on the unit 
circle is shown in 7). The expression in the first paren- 
thesis is exactly the same as 7;, obtained earlier for the 
transmission with the upper loop deleted. The second 
parenthesis, which is similar in form to the response 74, 


10 Each loop in a sampled system requires a delay to prevent the 
return of a signal to an element during the same iteration. If each 
loop were not so designed, the approximation of pulses by impulses 
would not be valid. 


36 IRE TRANSACTIONS ON AUTOMATIC CONTROL 


January 


(a) 


abz 
(1-z)?+abzk 


-a bz 
(1-z)*+abkz 


“U-Z)adZ_] ayes f_U-ziedz_] Xe 
(iz)? tabkz| "| U-z®+abkz 


(b) 


(1-Z)adZ : 
a ft 
(1-Z)* +a bkz 


abz 
(1-2)? +abzk 


(1-Z)adz M -abz 
(1-2)? +a bkz 


(1-2)? +abkz 
ples bare le ares ane 
ip Paezadz at. 2°27 (kK) +1 
[j-zi2+eorz | Z4(1-ad) +Z°( Ky) +1 J oot 
Liat J 
(c) 


Fig. 8—Reduction of Fig. 1. 


shows the effect of the correction of the upper loop at 
every nth sampling instant. Intuitive interpretations of 
this type are frequently possible. 


CONTINUOUS AND DISCRETE ELEMENTS; 
APPROXIMATIONS 


When continuous components occur together with 
discrete elements operating at different rates, the most 
rapid sampler determines the basic pulse rate; all others 
are designated as submultiple-rate pulses. After the con- 
tinuous elements have been individually starred at the 
basic rate, the resulting network is wholly discrete (ex- 
cept possibly at the terminal branches), and the analysis 
is similar to the previous two examples. This approach 
is illustrated in Fig. 9, which shows a control system 
whose output pulse rate is four times as rapid as the in- 
put rate; in addition, the feedback element is pulsed 
half as often as the output. Placing impulse modulators 
at the inputs and outputs of the continuous elements, 
H,(s) and H,(s), is redundant since these modulators 
operate more frequently than the Mj,,,’s; obviously a 
continuous element between two impulse or submulti- 
ple-rate modulators can be impulsively starred, as in 
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Fig. 9—Mixed system with three rates. 


Fig. 9(b), in which the redundancies have been removed. 
Node a is selected since pulses at this point are at the 
lowest rate in the circuit, and reduction to a single 
branch is obtained. The term 


{ Hi(s)*+ M12- H2(s)*} Mays, 
is evaluated by impulse starring of continuous functions, 
| Hi(2) + M1j2- H2(2)} Maps, 
and then by submultiple-rate starring. 
{ [Hi(z) ]*¥2- Ho(z)} #14, 


Since the element H(z) in the self-loop is preceded by 
Mi, and followed by My,2, insertion of another M4). be- 
tween Mi, and H(z) is redundant, and starring at 
1/n=} is valid! 

The inclusion of modulators whose rates are not inte- 
gral multiples of one another is representable by these 
techniques, but reductions with the complete elemina- 
tion of modulators is not always possible. In these cases 
approximations can frequently be made. For example, 
when a continuous filter is preceded by My; and fol- 
lowed by Mi;s, the latter modulator might be replaced 
by M,,9; of course, these substitutions are only valid 
when there is a continuous filter separating the modula- 
tors, and, even then, care should be exercised in these 
approximations. 


_ ™ In general, when a filter is followed by a modulator whose rate 
is a multiple of the input signal rate, the filter may be starred at the 
more rapid rate. 
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CONCLUSION 


Systems which lend themselves to representation by 
a multiple-rate sampled-data system model can be 
analyzed with facility by the methods developed above. 
These systems are characterized by the ordinary z- 
transform, using an extended “starring” technique. A 
table for the evaluation of submultiply starred functions 
(with either real or complex poles), as closed form ex- 
pressions in zg, allows the removal of modulators from 
the feedback loops. The result is a network composed 
of forward transmissions which may be analyzed by 
conventional methods. 

The stability and settling characteristics of an ele- 
ment are not affected by starring since the new poles 
produced by this operation always remain on a circle of 
the same radius as that of the original unstarred func- 
tion. It is seen that the evaluation of one representative 
response of the time-varying system will provide suf- 
ficient information for most analyses. The demonstra- 
tion that impulse modulation of continuous signals, re- 
sulting in the multiplication of poles in the s-plane, is 
analogous to the multiple-rate modulation of discrete 
signals resulting in a multiplication of poles in the 
z-plane, leads to a unified description of these modula- 
tion or sampling operations. 


APPENDIX 


The development of Table I’ was motivated by the 
need for a ready means of obtaining reduced rate and 
closed form equivalents of functions of zg. For example, 
when 


1 


1—iz 


=144s+ Ge)? + G+ 


is sampled at one-third of the basic rate, pulses are ob- 
tained at the Oth, 3rd, 6th, 9th, - - - , sampling instants; 


their heights are 1, (})*, (4)®, (3)°, -- +, respectively. 
Therefore, 
1 *1/3 : De 4 G ys 4 3 1 
a SARA To ~ 1— Ga) 


Attempts to obtain a single expression for equivalents, 
using a generic d, led to prohibitively lengthy terms. 
However, rather economical formulas were obtained for 
specific values of d up to d=4, which should be suf- 
ficient for almost all applications. The formulas were 
arrived at by long division of the general fraction 
z'/(1—bz)4, selection of terms, and arrangement into 
closed form. (The generic d has been left in Table I to 
show its form for future extensions.) 

To use Table I, the function is expanded into partial 
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fractions and substitutions are made rigorously. The 
summation notation should be interpreted as follows: 
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For example, 
Dee 20, ee ae 
4 2 


Evaluation of a function proceeds as follows: 


gt w1/4 
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>, bot + ( >» - 3 x) + >> 6812 
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Peat 
zt + 12b428 + 358312 
[1 — (bz)4]8 


Note the simple check that the sum of the coefficients of 
the numerators equals n@—); in the example above, 
>> (coef.) =16=4¢-», 

Another example illustrates the determination of a 
submultiple-rate equivalent of a function with complex 
roots: 


a (1 — 42)2z 
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The b’s are ($-+4/) and ($—4)): 
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Automatic Control of Three-Dimensional 


Vector Quantities—Part oe 
A. S. LANGE 


Summary—In Part 1 of this paper a vector algebra was de- 
veloped using a three-element column matrix to represent the vec- 
tor, and a three-by-three matrix to represent a vector transformation 
operator. Problems in spherical trigonometry were analyzed with 
the use of a position vector, and the design of automatic computers 
to solve such problems was considered. In Part 2, the angular 
velocity vector is introduced for the purpose of analyzing and de- 
signing geometric stabilization systems. 


V. GEOMETRIC STABILIZATION 
(Ci. stabilization is the process of iso- 


lating some controlled member, such as a gun or 

a radar antenna, from the motion of the platform 
on which it is mounted, such as a ship, an airplane, or a 
tank. The geometric stabilization system consists of a 
device to measure the motion of the base, and servo- 
mechanisms which move the controlled member with 
respect to the platform. The stabilization system func- 
tions in such a way that the motion produced by the 
servos is opposite in sense to that of the platform, so 
that the net motion of the controlled member with re- 
spect to some given set of reference coordinates is zero. 
In general, the platform has three degrees of (angular) 
freedom, so that it is necessary to have two or more 
degrees of freedom in mounting the controlled member 
on the platform. These several degrees of freedom are 
achieved by a series of gimbals, so that we speak of a 
three-gimbal mount or a three-axis mount, since each 
gimbal is designed so as to rotate with respect to the 
member on which it is mounted. Several of the more 
commonly used two- and three-axis antenna mounts are 
illustrated in Cady, et al.4 The six ways that the three 
axes of a three-gimbal mount can be arranged are shown 
in Fig. 6. Associated with each gimbal is a servomechan- 
ism, so that each gimbal can apply one restraint to the 
controlled member. When the controlled member is a 
gun or a telescope, there is a line associated with it (viz., 
the barrel of the gun, or the optical axis of the telescope) 
called the controlled line. In many instances, it is only 
necessary to isolate the controlled line from the motion 
of the base, rather than the controlled member itself. 
For this purpose, a two-axis mount is sufficient, since it 
is usually not required that a gun barrel not rotate about 
its own axis. However, if the controlled member itself 


* Manuscript received by the PGAC, August 10, 1959. 
{ Systems Div., Bendix Aviation Corp., Ann Arbor, Mich. 
13 A, S. Lange, “Automatic control of three-dimensional vector 


quantities—Part 1,” IRE Trans. on AUTOMATIC CONTROL, vol. » 


AC-4, pp. 21-30; May, 1959. 

14 W. M. Cady, M. B. Karelitz, and L. A. Turner, “Radar Scan- 
ners and Radomes,” M.I.T. Rad. Lab. Ser., McGraw-Hill Book Co., 
Inc., New York, N. Y., vol. 26 p. 106, Fig. 4.1; 1948. 


must be irrotational with respect to the reference frame, 
then three degrees of freedom are required. 

The requirement that the controlled member be iso- 
lated from the motion of the base on which it is mounted 
may be expressed by the vector equation 


(29) 


where Wi.=angular rate of S, with respect to S;; Sq is 
a set of coordinates fixed to the controlled member, and 
S; is a set of reference coordinates. It may be recalled 
that in the previous section, it was necessary to use posi- 
tion vectors in order to describe angular quantities, be- 
cause angles do not obey the rules of the usual vector 
algebra. For our purposes, however, angular rates follow 
these rules, and hence may be treated as vector quanti- 
ties, whose direction is given by the right-hand rule, and 
whose length is given by the magnitude of the rotation.® 

Eq. (29), then, states that the angular velocity of S, 
(a coordinate set which defines the controlled member) 
with respect to S; (a coordinate set which defines the 
reference frame) must be zero. Eq. (29) may be written 
as 


Wie = Wa t+ Wa = 0 
or 
(30) 


where Wa is the angular velocity of Sz (a coordinate set 
which defines the platform on which the controlled 
member is mounted) with respect to S;, and Waa is the 
angular velocity of Sa with respect to Sg. Eq. (30) shows 
that if the controlled member, defined by Su, is to be 
stationary with respect to the reference frame S;, then 
the angular velocity of the controlled member with re- 
spect to the platform Waa, must be equal to the angular 
velocity of the platform Sz with respect to S;, Wia. 

One of the design problems associated with base mo- 
tion isolation systems is the arrangement or order of the 
gimbals. In some cases, the designer has no choice in 
selecting the configuration; however, some aspects of 
the design which may influence the choice of gimbal 
sequence will be discussed in Section VIII. 

The designer of an automatic control system is also 
concerned with defining such servo characteristics as 
bandwidth, maximum velocity, maximum acceleration, 
gear ratio, torque, and horsepower. In the following dis- 
cussion, the determination of the angular velocity re- 


1 H. Goldstein, “Classical Mechanics,” Addison-Wesley Co., 
Reading, Mass., pp. 124-134; 1950. 
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Fig. 6—Types of three-axis mounts. Key to symbols: Gzi=gimbal angle; z indicates system i=1, 2, 3+ ++ 
out) 7=1, 2, 3. 1: between base and outer gimbal (0). 2: between outer gimbal (0) and middle gimbal ( 


and inner gimbal (J). 


quired is considered in detail, and the methods used for 
determining the other servo parameters are considered 
only briefly. 


VI. STABILIZING THE Two-Axis Mount! 


Consider the two axis mount discussed in Section ITI 
of Part 1. It has been stated that two-axis mounts are 
used to isolate a controlled line from platform motion, 
and we consider here some details concerning the opera- 
tion of such a mount. Fig. 2(a) of Part 1 illustrates the 
space-flow diagram for this mount, and Fig. 2(b) is a 
three-dimensional sketch of the angles which describe 
the problem; the definition of these angles and the co- 
ordinate sets associated with them given in Part 1 is 
repeated here.!” 

Let the target be located in space by the angles EF and 
B and range, where B (bearing angle) is measured in the 
horizontal plane from the bow of the platform, and EL 
(elevation angle) is measured about an axis in the hori- 
zontal plane. Let the instantaneous attitude of the plat- 
form be given by Ezo (pitch angle) and Zo (roll angle). 
Eio is measured about an axis in the horizontal plane, 
and Zo is measured about the longitudinal axis of the 
platform. Further, let the telescope or antenna be posi- 
tioned with respect to the platform by two angles, Bd’ 
(antenna train angle) and Ed’ (antenna elevation angle), 


16 Cady, et al., op. cit., Fig. 4.1(b). 

17 The auxiliary sets of coordinates located between S; and Sg are 
of no interest here, since we assume that the motion of Sz with respect 
to S; can be specified without regard to how these sets are defined. 
This topic is discussed in more detail in Appendix II. 
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Fig. 2 (Part 1)—(a) Space-flow diagram. (b) Coordinate 
converter geometry. 
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where Bd’ is about an axis normal to the deck of the 
platform and Ed’ is about an axis in the deck of the 
platform. Fig. 2(a) is the space-flow diagram which de- 
scribes this problem, and Fig. 2(b) is a three-dimensional 
sketch of the problem. Fig. 2(a) contains virtually all 
the information presented in Fig. 2(b), and at the same 
time, serves to indicate the flow of information in a way 
which is useful to the design of the coordinate converter. 
The S’s in Fig. 2(a) represent sets of Cartesian coordi- 
nates which may be defined as follows: 


S;=a reference set of coordinates, with the (xy), 
plane horizontal, and the z; axis directed ver- 
tically down. 

S,=an auxiliary set of coordinates which results from 
a rotation B about 2;; 3;=z, and the plane (xy); 
is coincident with the plane (xy),;. The (xz), 
plane, like the (xz); plane, is vertical, and con- 
tains Rzs and hence, the target. 

S;=target coordinates, resulting from a rotation 
through the angle E about y, so that y,=y23 x4 
is directed toward the target, so that Rx is along 
xz, and z; forms an orthogonal set. 

S,=an auxiliary set which results from the rotation 
through the angle Lio about y;=y,; x, is the 
longitudinal axis of the platform. The (xz), plane 
is vertical and y, is horizontal. 

Sa= platform coordinates, which result from the rota- 
tion Zo about x, so that (xy)a defines the deck of 
the platform; if the platform is thought of as a 
ship, then the mast is directed along the negative 
Za axis. 

S;=a set of coordinates fixed to the outer gimbal, or 
train member. S, results from the rotation Bd’ 
about 2¢ so that zgz=2,, and the (xy), plane is co- 
incident with the (xy)¢ plane. 

S,=a set of coordinates fixed to the inner gimbal or 
elevation member, resulting from the rotation 
Ed’ about the , axis so that y,=7ya. The vector 
Rrz is directed along the x, axis, and is therefore 
contained in the (xz), plane. 


It is customary to consider S;as the Newtonian frame 
with respect to which the gyroscopic instruments meas- 
ure motion. A Newtonian frame is usually defined,!® in 
essence, as a frame in which Newton’s second law is 
applicable. The physical significance of this somewhat 
circular definition is that the selection of a Newtonian 
frame depends on the particular problem at hand. For 
example, in the present problem of isolating some con- 
trolled member from the motion of the base on which it 
is mounted, it is usually sufficient to consider that the 
reference frame is moving with the uniform transla- 
tional motion of the platform, but not participating 
with the platform’s angular motion. However, if the 
problem is to design an inertial navigation system 


18 J. L. Synge and B. A. Griffith, “Principles of Mechanics,” Mc- 
Graw-Hill Book Co., Inc., New York, N. Y., p. 32; 1959. 
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which must indicate position very accurately for many 
hours, it is necessary to use, as a Newtonian frame, a 
set of coordinates centered at the earth, but not rotating 
with it. Presumably, in order to design an inertial sys- 
tem for extraterrestrial navigation, it would be neces- 
sary to use a Newtonian frame defined by the fixed 
stars. The point here is that there are inaccuracies intro- 
duced by the choice of the Newtonian frame which may 
be made arbitrarily small by the proper choice of a 
Newtonian frame. When stabilizing a gun, for example, 
the errors introduced by neglecting the rotation of the 
earth or the translational motion of the platform are 
essentially undetectable. 

It was previously stated that the requirement that 
the controlled member be isolated from the motion of 
the base may be written vectorially as 


Wea =o Wia. 


But, from Fig. 2(a), as well as the definitions listed 
above, it can be seen that Waa= Wa+ Wea. Further, if 
this vector equation is to be useful for design purposes, 
it must be expressed in terms of components which 
represent measurable quantities. Therefore, (29) may be 
written as 


Wa + Wea = — Wis, (31) 
where the superscript a indicates that the vectors are to 
be expressed in their S, components. The three vectors 
in (31) will each be examined in detail. Consider first the 
term Was: 


—=0 ae 
Was = TasWar 


cos Ed’ 0 -—sin Ed’ 0 
ay Be era 0 0 (32) 
sin Ed’ 0 cos Ed’ || DBd’ 


where D (_)!° is the differential operator d/dt ( ). Ta is 
the transformation matrix which operates on a vector 
expressed in S, components in order to express that 
vector in its S, components. Since, as can be seen by 
Fig. 2(a), Bd’ is the angle between S; and Sz, and since 
this angle represents a rotation about the %, zg axis, it 
follows that |Wa|, the magnitude of Wa, is equal to 
DBd’, and that the direction of Wa is along the % axis. 

Likewise, the second term in (31) may be expressed as 


0 
DEd' | ; 
0 


We = (33) 


1.e., the angular velocity of S, with respect to S; is equal 
to the time rate of change of the angle between them. 
As can be seen from Fig. 2(a), this angle is Ed’, so 
| Wea| =DEd’, and since the rotation Ed’ occurs about 
the y, axis the direction of Ws, is along that axis. 


19 Following “Standard Fire Control Symbols,” Dept. of -the 
Navy, Washington, D. C. BuOrd Publication No. OP-1700; 1950. 
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Finally, the last term in (31) may be written as fol- 
lows: 


cos Ed’ 0 —sin Ed’ 
Wia = Die DigW sa = 0 1 0 
sin Ed’ 0 cos Ed’ 
cos Bd’ sin Bd’ 01|| w, 
-|—sin Bd’ cos Bd’ 0 || w, (34) 
0 0 1 {| w, 


where 7}, is the transformation matrix required to trans- 
form a vector expressed in Sz; components into its Sy 
components, and w,, wy, and w, are the x, y, and z com- 
ponents, respectively, of Wiain Sz. (The several ways in 
which these components may be expressed are discussed 
in Appendix IT.) 

It is worthwhile to observe that (34) is more easily 
evaluated if the operations are performed from right to 
left. That is, first evaluate 


wz cos Bd’ + w, sin Bd’ 


Wo TW — | Ho. sin Ba’ + w, cos Bd’ 
Wz 
and then 
Wia= Taye 


(w, cos Bd'+w, sin Bd’) cos Ed’ —w, sin Ed’ 
= —w, sin Bd’+w, cos Bd’ 
(w, cos Bd’+w, sin Bd’) sin Ed’+-w, cos Ed’ 


. (35) 


In contrast, if (34) is expanded from left to right, we 
have, first, 


Dap = Tel -ta 
cos Ed’ cos Bd’ cos Ed’ sin Bd’ —sin Ed’ 
a —sin Bd’ cos Bd’ 0 
sin Ed’ cos Bd’ sin Ed’ sin Bd’ cos Bd’ 
and secondly, 
= 2 —— a 
W a= TarW ia 
w, cos Ed’ cos Bd’+w, cos Ed’ sin Bd’—w, sin Ed’ 
= —w, sin Bd’+w, cos Bd’ . (36) 


w, sin Ed’ cos Bd’+w, sin Ed’ sin Bd’+w, cos Ed’ 


Of course, (35) and (36) are identical, but the steps fol- 
lowed in obtaining (35) are simpler because the result 
of each operation is a 3-element column matrix, whereas 
the first step in obtaining results is a 3X3 square matrix. 
Such matrices become more and more cumbersome with 
each successive operation, so that the right-to-left se- 
quence is even more important when a number of trans- 


formations are to be performed. 
4 
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When the indicated operations are performed, (32)— 
(34) can be substituted into (31), giving 


— DBd' sin Ed’ 
_— DEd’ 
DBd' cos Ed’ 


(w, cos Bd’ + w, sin Bd’) cos Ed’ — w, sin Ed! 


a —w, sin Bd’ + w, cos Bd’ (37) 
(w, cos Bd’ + w,sin Ba’) sin Ed’ + w, cos Ed’ 
so that 
DEd’ = w, sin Bd’ — wy, cos Bd’, (38) 
DBd' = — [w. + (wz cos Bd’ + w, sin Bd’) tan Ed’]. (39) 


Eqs. (38) and (39) express the outputs of the elevation 
and train servos, respectively, which are required to iso- 
late the controlled member from the motion of the base. 
It may be observed that Wia is not exactly zero; i.e., 


(w, cos Bd’ + wy, sin Bd’) sec Ed’ 
(40) 


However, if x. (the x axis of the coordinate set S,) is the 
barrel of a gun, or the tracking line of a radar antenna, 
or a telescope, then the controlled line itself is stationary 
with respect to S;, although the gun or the antenna (the 
controlled member) is rotating, with respect to S;, about 
the controlled line. In many instances, this rotation is 
not important. Note that the vector equation W;.=0 
expresses three conditions which must be satisfied. How- 
ever, the two-gimbal mount supplies only two con- 
straints, corresponding to the outputs of the train and 
elevation servos. If it is desired to satisfy (31) identi- 
cally, it is necessary to introduce a third constraint, in 
the form of another gimbal. A three-gimbal mount 
which satisfies (31) is discussed in the following section. 
It will suffice here to observe that to stabilize a line, two 
gimbals are required, and to stabilize a plane, three 
gimbals are required. (It is presumed in this statement 
that there is no redundancy in the gimbal system. Re- 
dundancy has been considered in Section IV of Part 1 
and is discussed further below). 

It is necessary, next, to examine in detail how a 
mechanism capable of solving (38) and (39) can be in- 
strumented. That is, (38) and (39) represent the values 
of DEd’ and DBd’ required to stabilize the tracking line, 
but these are the outputs of the elevation and train 
servos. What are their inputs, and how are the servo 
loops closed? Two methods using gyroscopic instru- 
ments are common: 1) measure the motion of the plat- 
form with respect to inertial space, or 2) measure the 
motion of the tracking line with respect to inertial space. 
The electrical signals from the gyroscopic instruments 
form the inputs to the servos, and the motion of the 
tracking line, in one case with respect to the platform, 
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and, in the other, with respect to inertial space, act to 
close the servo loops. One of these methods is called 
“director stabilization,” and the other is called “on- 
mount stabilization.” In the director system, a vertical 
gyro (or, in general, a two-gimbal gyro) is used to meas- 
ure the motion of the base with respect to inertial space. 
The signals from the gyro are modified by means of a 
coordinate converter (see Part 1, Section III) to obtain 
the servo inputs. In fact, as was pointed out, if the 
servos which drive the BD’ and Ed’ resolvers in Fig. 
3(b) of Part I are also used to position the controlled 
member, the controlled line is isolated from the motion 
of the base; it can be seen that (38) and (39) can be ob- 
tained by differentiating (21) and (17) of Part I, re- 
spectively. The name director stabilization is derived 
from the fact that the vertical gyro and the coordinate 
converter are usually located in the fire control director. 

The name on-mount stabilization, on the other hand, 
is derived from the fact that the gyros are mounted 
directly on the controlled member. Single gimbal, or 
rate gyros have been widely used for this purpose,?° and 
this is the configuration which will be assumed in the 
following discussion. Fig. 7 is a servo block diagram of 
one axis of the base motion isolation system utilizing on- 
mount stabilization. Note that the servo loop is closed 
by means of the motion of the controlled member with 
respect to inertial space, as measured by the gyroscope. 


INTEGRATING RATE GYRO Wa 


Fig. 7—Servo block diagram for stabilization loop (single axis). 
W.= Wia commanded, Wie =angular rate of Sa with respect to S;, 
W,a=angular rate of Sg with respect to S;, Waa=angular rate of 
Sa with respect to Sq. 


The rate gyros (a detailed discussion of the operation 
of the rate gyro will be given in Part 3) mounted on the 
controlled member measure the y and zg components of 
W ia. (In this case, the x component is not measured for 
base motion isolation purposes because there is no way 
to constrain it to zero.) The outputs of the y and z gyros 
can be determined from the expression for W%,, where 


Wis = Wia + Was + Woo. 
Then, from (32)—(34) 


(Wia)y = DEd' — wz sin Bd’ + w, cos Ba’ (41) 
(Wia)s = (DBd' + w,) cos Ed’ 
+ (wz cos Bd’ + w, sin Bd’) sin Ed’ (42) 


20 C. S. Draper, et al., “The Floating Integrating Gyro, and its 
Application to Geometrical Stabilization Problems on Moving 
Bases,” Inst. Aeronaut. Sciences, New York, N. Y., Rept. No. 503; 
January, 1955. 
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where (Wia)4 is the y component of Wi, measured by 
the y-axis gyro, and (Wi)? is the g component of Wi, 
measured by the z-axis gyro. 

Consider first the elevation servo; from the block dia- 
gram shown in Fig. 7, the output of the elevation servo 
is 


(KiKo)e 


fe a3 ae 3 
DEd ‘ [Wee — (Wia)y] (43) 


where K, is the gyro sensitivity, Ke is the servo gain, and 
and the subscript e designates the elevation servo. 

Now, if W.. (the commanded elevation rate) is zero, 
then (43) becomes 


Dida (wz sin Bd’ — w, cos Bd’) 


ip 


(44) 


(Kika). 


and for (K,K2)e;>D, (44) may be approximated by 


DEd' ~ (w, sin Bd’ — wy, cos Bd’). (45) 
In a similar fashion, for the train servo, 
ppg Soke yee (46) 
D 
where the subscript ¢ designates the train servo. 
Roraye;— 0) 
DBd' = — a ee [w, + (w, cos Bd 
D sec Ed’ 
(Kik2) 
+ w, sin Bd’) tan Ed’]. (47) 
As before, if (K,K2); cos Ed’>>D, (47) becomes 
— DBd' = |w, + (wz cos Bd’ + w, sin Bd’) tan Ed’]. (48) 


Since (45) ~ (38) and (48) ~ (36), it can be seen that 
the train and elevation servos do, in fact, stabilize the 
controlled line when they are suitably connected to the 
y and zg gyro outputs, respectively. 

It may be observed from (47) that the time constant 
of the train servo varies inversely with cosine of the ele- 
vation angle Ed’. This is because the input axis of the 
z-axis gyro is misaligned with the train servo output axis 
by the angle Ed’. As a consequence, the loop gain is 
low for large values of Ed’ (going to zero as Ed’—7/2) 
and the performance of the train servo varies with Ed’. 
If Ed’ varies over a significant range of values, it is cus- 
tomary to insert an automatic gain control (AGC) de- 
vice to eliminate this effect. The automatic gain control 
may be accomplished by multiplying the output of the 
z-axis gyro by sec Ed’. This may be done with a resolver 
(see Appendix I) but usually it is sufficient to use a non- 
linear potentiometer, wound so as to modify the applied 
voltage by the secant of the potentiometer shaft angle. 
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If the shaft is driven by Ed’, (46) becomes 


(K,Ko)t 


DBd' = [Wet — (Wia)z] sec Ed’ (49) 


from which it follows that 


1 
DBd' = — 


a cee to a 
(KiK2): aE 


[w, + (w, cos Bd’ + wy sin Bd’)|; (50) 


1.e., the variation of the train servo time constant with 
Ed’ has been eliminated. It must be emphasized that the 
purpose of this operation is to maintain a constant loop 
gain for wide variations in Ed’. This AGC function is 
not so clear in three axis mounts, so it is important that 
this function be clearly understood in this more simple 
example. The function sec Ed’ is difficult to mechanize 
as Ed’—7/2, since sec 7/2—> ». However, it can be seen 
that for Ed’—7/2, this two-axis mount exhibits a trig- 
onometric singularity which means that operation for 
Fd’ near 7/2 must be avoided. Therefore, the AGC 
operation will be performed satisfactorily for all values 
of Ed’ for which the system may be expected to operate. 


VII. A THrREE-AxtIs STABLE PLATFORM 


It was shown in the previous section that a two-axis 
gimbal arrangement, although geometrically stabilizing 
a controlled line, allowed rotation of the controlled 
member about the controlled line, where, it may be re- 
called, the controlled line is an axis fixed to the con- 
trolled member. It is of interest, therefore, to consider a 
three-axis gimbal arrangement, which has three condi- 
tions of restraint, one of which can be used to restrain 
(Wia)2* to zero, thus eliminating the rotation of the con- 
trolled member about the controlled line. 

As an example, we will take for our model the cross- 
level gimbal arrangement# shown as System I of Fig. 6. 
The “cross-level” arrangement is similar to the train- 
elevation system described in the previous section, ex- 
cept that an extra gimbal has been added between the 
train and the elevation gimbals. The function of this 
extra gimbal is to maintain the elevation axis level, or 
horizontal. The name cross-level refers to the extra axis. 
which is oriented 90° with the elevation axis; hence it 
crosses with the level axis. Notice that the elevation axis 
_may be thought of as the controlled line of a two-axis 
mount; hence, it can be demonstrated by arguments 
similar to those used in Section VI that the controlled 
line (the elevation axis) is isolated from the motion of 
the base. It follows, therefore, that if a controlled mem- 
ber is mounted on this axis, the controlled member can be 
completely isolated from the motion of the base, since 
it can rotate about the controlled line in such a way as 
to remain fixed in space. 


1 Cady, et al., op. cit., Fig. 4. 
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The space-flow diagram which represents this gimbal 


configuration is shown in Fig. 8. The S’s in Fig. 8 repre- 
sent sets of Cartesian coordinates which are defined as 


follows. 


S:=a reference set of coordinates, with the (xy); 
plane horizontal, and the z; axis directed ver- 
tically down. 

Sa=platform coordinates, with (xy)¢ defining the 
deck of the platform. 

»=a set of coordinates fixed to the outer gimbal or 
train member. S, results from the rotation Bd 
about zg so that z= and the (xy), plane is co- 
incident with the (xy) plane. 

S,=a set of coordinates fixed to the middle or cross- 
level gimbal. S, results from the rotation Zd 
about the x, axis so that x,=x-. The (xy), plane 
is horizontal. 

S,=a set of coordinates fixed to the inner or elevation 
gimbal. S, defines the controlled member, and 
results from the rotation Ed about the y, axis so 
that ye=Ya. Since the (xy). plane is horizontal 
the (xz)_ plane is vertical. 


ANY 
CONVENIENT 


Fig. 8—Space-flow diagram cross-level system. (See Appendix II.) 


Bd and Ed are the train and elevation angles, re- 
spectively, analogous to Bd’ and Ed’. (In general, 
primed quantities are used in OP-1700'* to represent 
“unstabilized” quantities, as distinguished from “sta- 
bilized” quantities which are so named because the 
angle Ed lies in the vertical plane (xz),, which remains 
fixed with respect to the reference space when directed 
at a stationary target.) 

The requirement that the controlled member be iso- 
lated from the motion of the base is given by the follow- 


ing vector equation: 
Win = 0. (Sta) 


But, from Fig. 8, Wia=WiatWatWee+Wea so that 


(51a) may be written 
Was t+ Wo + Waa = — Wia (51b) 


where the superscript a indicates that the vectors are to 
be expressed in S, coordinates, and 


Wa = angular rate of Sa with respect to S;, 
Wap = angular rate of S, with respect to Sa, 


Ww = angular rate of S, with respect to Ss, 
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angular rate of S, with respect to S¢. 
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The four vectors in (51b) will each be examined as fol- 
lows: 


cosHd Q -—sin Ed 
Wee AW aetna Oats hake 
sin Hd 0 cos Ed 
1 0 0 0 


Omcos Za..cin aa. 0 
0 —sinZd cosZd|| DBd 


(52) 


Ta. is the transformation matrix which operates on a 
vector expressed in S, in order to express it in S,, and 
T.» is the matrix transforming components from S, to 
S,. Since Bd is the angle between Sz and S;, | Wa|, the 
magnitude of Wa, is equal to DBd, and since the com- 
mon axis between Sz and 5; is the g axis, it is directed 
along this axis. 


cos Hd 0 -—sin Ed DZd 
Wie=TuWie=| 0 1 0 0 (53) 
sin Hd 0 cos Ed 0 


Zd is the angle between S, and S,, with x,»=x,. There- 
fore, | Wee| = DZd, and W,, is directed along the x direc- 
tion in both S, and S.. 


0 
DEd |. 
0 


Wace =a (54) 


Fig. 8 shows that the angle between S, and S, is Ed, and 
the axis of rotation is their common y axis. Therefore, 
the rate of change of Ed is given by DEd, directed along 
Ya: 


=o ——d 
W ia = Tacl cb! baW ia 


cosHd O -—sin Ed 1 0 0 
= 0 1 0 0 cosZd_ sin Zd 
sin Ed 0 cos Ed || 0 —sinZd cosZd 


cos Bd sin Bd 0|]| wz 
—sin Bd cosBd 0|| wy 
0 0 1 || w, 


(55) 


where 7}q4 transforms the Sg components of Wig into S, 
components. Note that in Fig. 8, no particular set of 
angles between S; and Sq is specified. Wi represents a 
measurable physical quantity, and the representation 
shown in Fig. 8 is intended to emphasize that the vector 
Wa is completely independent of the various ways it 
can be described mathematically. It is customary to 
speak of the three angles relating the attitude of an 
object with respect to some reference coordinates as 
“Rulerian angles.”” Expressions for Wia in terms of 
Eulerian angles are discussed in more detail in Appen- 
cx.) Te 


2 Goldstein, op. cit., pp. 107-109. 
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Eq. (51) can be rewritten by combining (52)-(55), 
giving 
DZd cos Ed — DBd cos Zd sin Ed 
DEd + DBd sin Zd 
DZd sin Ed + DBd cos Zd cos Ed 
X cos Ed — w, cos Zd sin Ed 
+ Y sin Zd sin Ed 


=— Y cos Zd + W, sin Zd 
X sin Hd + w, cos Zd cos Ed 
— Ysin Zd cos Ed 


(51ic) 


where 
X = (wz cos Bd + w, sin Bd), 
Y = — (w,sin Bd — w, cos Bd). 


The following scalar equations may be obtained from 
(ate): 


DZd cos Ed — DBd cos Zd sin Ed 
= — |—X cos Ed + Y sin Zd sin Ed 
— w, cos Zd sin Ed| (56) 
DEd + BDd sin Zd = — [VY cosZd+w-sin Zd], (57) 
DZd sin Ed + DBd cos Zd cos Ed 
= — [X sin Ed —Y sin Zd cos Ed + w, cos Zd cos Ed| (58) 


From these three equations we seek expressions for 
DZd, DBd, and DEd, which are the outputs of the 
servos which drive the three gimbals. These expressions 
may be found by the following steps. 

1) The sum of cos Ed (56) and sin Ed (58) gives 


DZd = — X = — (wz, cos Bd + w,sin Bd). (59) 


2) The sum of—sin Ed (56) and cos Ed (58) gives 


DBd = — w,+ Y tan Zd 
= — [w. + (wz sin Bd — wy cos Bd) tan Zd]. (60) 


3) Finally, the substitution of (60) into (57) gives 


DEd = sec Zd(w,z sin Bd — w, cos Bd), (61) 


Eqs. (59)—(61) express the outputs of the three gimbal 
servos, for the middle gimbal, outer gimbal and inner 
gimbal, respectively, which are required to make 
Wia=0. The work of deriving these expressions may be 
simplified somewhat by recalling that the function of the 
cross-level servo drive, whose output is DZd, is to main- 
tain the elevation axis y,, horizontal. Mathematically, 
this may be expressed as (W;.)$=(Wie)?=0. That is, 
if the cross-level gimbal is not rotating about its x axis 
with respect to the reference frame, then its y axis is not 
rotating with respect to the reference frame. Therefore, 
if we express W}, as 
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——b ———6 = =e 
Wie = TuWia + Was + TW 


wz cos Bd + w,sin Bd 0 DZd 
=| —w,sin Bd + w,cosBd|+|} 0 |+!] 0 (62) 
Wz DBd 0 
then the condition that (W;.)3=0 gives 
DZd = — (w, cos Bd + w, sin Bd) (63) 
and (62) reduces to 
0 
Ws = |—w,sin Bd + w, cos Bd (64) 
DBd+ w, 


The vector equation specifying that the controlled 
member S, be irrotational with respect to the reference 
frame is given by 


Wis = FET aw xe Wen = 0 


which may be expressed in terms of the now simplified 
W.- to give 


Y sin Zd sin Ed—(BDd+w-) cos Zd sin Ed 
Y cos Zd+ (DBd+w.)+sin Zd+ DEd = 


0 
0}. (65) 
— Y sin Zd cos Ed+-(DBd-+-w,) cos Zd cos Ed 0 


Either the x or the z component of (65) can be used to 
give the expression for DBd as 


DBd = — [w, + (wz sin Bd — w, cos Bd) tan Zd]. (66) 


Eq. (66) may be substituted into the y component of 
(65), to give 


DEd = sec Zd(w, sin Bd — w, cos Bd). (67) 


Eq. (63), (66), and (67) are of course identical to (59)— 
(61), respectively. However, it may be noted that the 
second method of derivation is more simple and requires 
fewer algebraic operations; also, this second method 
demonstrates the “levelling” operation of the middle 
gimbal servo. 

The discussion thus far has been concerned with the 
determination of the gimbal servo rates required to iso- 
late the controlled member S, from the motion of the 
base Sz. The next step is to determine how the servos 
must be arranged to accomplish this base motion iso- 
lation. From (51c), the outputs of x, y, and z gyros 
mounted on S, may be inferred as 


(Wia)s = DZd cos Ed — DBd cos Zd sin Ed + X cos Ed 
+ Ysin Zd sin Ed — w, cos Zd sin Ed (68) 
DEd + DBdsin Zd + Y cosZd+ w,sin Zd (69) 


DZd sin Ed + DBd cos Zd cos Ed + X sin Ed 
— V sin Zd cos Ed + w, cos Zd cos Ed. (70) 


(W ia)y 
(W ia) 


I 


The servo output quantity is equal to K,K2/D times the 
difference between the commanded rate and the actual 
rate measured by the gyro, where K; is the gyro sensi- 
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tivity, and Ky is the servo gain, as is shown in Fig. 7. 
Therefore, if the subscripts J, M, and O are used to 
designate the inner, middle, and outer gimbal servos, 
respectively, we have 


(Ki Kk») 


Did= [Wee — (Wia)a (71) 
DEd = SEs [We — (Wia)yl (72) 
DBd = ae [Wee — (Wia)e] (73) 


where, as before, D=d/dt ( ). 
Now, for W.., Wye, and W,,, the x, y, and zg commands, 
respectively, all equal to zero, we have 


DZd = [(DBd + w.) cos Zd tan Ed 
D sec Ed 
(KiKe) 
— X — Vsin Zd tan Ed] (74a) 
1 
Did = 
eae eh 
(KiK2)r 
[—(DBd + wz) sin Zi — Y cos Zd| (75a) 
1 
DBd = ———_—_—_——_ [—-w, + Y tan Zd 
D sec Zd sec Ed 
(KiK2)o 


— (DZd + X) sec Zd tan Ed]. (76a) 


It can be observed that the servo loop time constant for 
the cross-level and train servos, (74) and (76), respec- 
tively, varies with the angle Ed, for the former, and with 
the angles Ed and Zd, for the latter. This variation in 
the loop gain can be compensated for in the same man- 
ner described in Section VIII, by the use of secant po- 
tentiometers for the automatic gain control operation. 
However, before discussing the AGC operator require- 
ments, it is of interest to first demonstrate that, in 
theory at least, the gimbal drives function properly even 
without the AGC feature. 
Assuine that 


(K1K2)u > D sec Ed 
(KiK2); = D 
(KiKe)o = D sec Ed sec Zd; 


I 


1.e., assume that 
sec Ed 
——— D*Zd 
(K1K2)r 
is much smaller that DZd and 
sec Ed sec Zd 


(K1K2)o 


D° Bd 
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is much smaller than DBd. If the system is stable, so 
that DBd and DZd reach some equilibrium value, then 
this is a valid assumption, and our purpose in making 
this assumption is to examine the equilibrium behavior 
of the system. (The determination of stability is a more 
difficult problem; the consideration of stability in a 
simple example is considered in Appendix I). These as- 
sumptions are somewhat analogous to the application 
of the final value theorem of Laplace transformation 
theory. Using these assumptions, (74a), (75a) and 
(76a) become, respectively, 


DZd = [|(DBd + w,) cos Zd tan Ed 
— X — Y sin Zd tan Ed], --(74b) 
DEd =~ |—(DBd + w,) sin Zd — Y cos Zd], (75b) 
DBd = [—w, + Y tan Zd 
+ (DZd + X) sec Zd tan Ed]. (76b) 


The sum of Equation (74b) and cos Zd tan Ed (76b) 
gives 

ae ee an, Con Bd secon, Sin Ba (90) 
The substitution of (77) into either (74b) or (76b) gives 
DBij-— "0, +0) tan Zd 


[—w, + (wz sin Bd — w, cos Bd) tan Zd| 


I 


(78) 
and finally, the substitution of (78) into (75b) gives 


DEd = sec Zd(w; sin Bd — wy cos Bd). (79) 


Since (77)—(79) are identical to (59)—(61), respectively, 
it can be concluded that the fact that the gyro axes and 
the servo axes are not always aligned does not effect the 
steady-state performance of the system. However, the 
degeneration of the dynamic performance implicit (74a) 
and (76a) is not generally tolerable, so that it is neces- 
sary to compensate for the misalignment of the gyro and 
servo axes by multiplying the output of the gyro by the 
secant of the angles of misalignment. For example, (71) 
becomes 


(KiK2)u 


DZd = [Wee — (Wia)2| sec Ed (80a) 


and for Wz.=0, we have 


DZd = [(DBd + w,) cos Zd tan Ed 
—— +1 
(KiK2)u 


— X — Ysin Zdtan Ed]. (80b) 


#3 M.F. Gardner and J. L. Barnes, “Transients in Linear Sys- 
tems,” John Wiley and Sons, Inc., New York, N. Y., vol. 1, pp. 
265-267; 1942. 
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Since the input axis of the y gryo is always aligned with 
the inner gimbal servo output axis, no AGC is required 
for the elevation servo. However, the train servo axis is 
misaligned by both the Ed and Zd angles, so that the 
AGC operator used to correct the z axis gyro has the 
form sec Ed sec Zd. Therefore (73) becomes 


KiK 
pBe state [Wee — (Wia)z| sec Ed sec Zd. (81a) 
For W,.=0 (53a) becomes 
DBd = [—w, + Y tan Zd 
pte: ees BS -h 1 
(KiKe)o 
— (DZd + X) sec Zd tan Ed]. (81b) 


It can be seen that the term in the brackets of (80b) is the 
same as the similar term in (74a); likewise, the brack- 
eted term in (81b) is identical to the bracketed term in 
(76a). Therefore, it may be concluded that the AGC 
operator effects only the servo-loop time constant, and 
does not change the steady-state performance of the 
system. 

In many applications, Ed and Zd are relatively small 
angles (7.e., less than 7/2). In such cases, it is con- 
venient to use potentiometers to perform the AGC func- 
tion, since potentiometers are less expensive, lighter and 
smaller than resolvers. In fact, since high accuracy is not 
required for this operation, linear potentiometers may 
be used.% However, if Ed or Zd perform larger excur- 
sions, the use of resolvers is required to adequately per- 
form these AGC operations. 

The location of the resolvers in the stabilization sys- 
tem may be determined by the following considerations. 
The three gyros measure the x, y and zg components of 
Wj,. The components of W,. in Sc may be obtained by 
applying (Wia); and (Wia)f to a resolver driven by Ed, 


since (Wia)2= (Wa); Then 
Win = T, sa W ie 
cos Ed 0 sin Ed (Wia)s 
= 0) 1 0 (Wia)y 
—sin Ed 0 cos Ed || (Wia)e 


(Wiale cos Ed + (Wiz), sin Ed 
{ 


= (Wia)y 
= Was) sin Ed + (Wia)e cos Ed 
(Win)e 

as (Wia)y 
(Wia)e 


*4 J. A, Greenwood, Jr., J. V. Holdam, Jr., and D. MacRae, Jr., 
“Electronic Instruments,” M.I.T. Rad. Lab. Ser., McGraw-Hill 
Book Co., Inc., New York, N. Y., vol. 21, pp. 102-103; 1948. 
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The x component of Wj, can be evaluated from (68) 

and (70) to give 
(Wa D7d-- xX. (82) 


If (Wia)z is applied to the cross-roll servo, the servo 
output may be written as 


(K1Ke)u 


DZd = (Wa Didi xX) 


which reduces to the familiar form 


Wee X 


DZd = mWee x) (83) 


(KiKo)u = ; 
and for W.,,=0, DZd~—X, which is the same result 
found earlier in (59), (63), and (77). Inspection of (83) 
shows that the Zd servo gain is independent of the value 
of Ed. Consequently we may conclude that the process 
of generating the Zd servo signal by passing the X- and 
Z-axis gyro signals through an Ed resolver is the correct 
operation. Since DZd ~ —X, it follows that (Wi.)$=0. 
The term (Wia)$ can be evaluated from (67) and (69) 
to give 


sin? Zd 


cos Zd 


(Wia)y = (DBd + wz) sin Zd — Y (34) 


Likewise, the z component of Wj, can be evaluated from 
(68) and (70) to give 


(Wia)e = (DBd + w,) cos Zd — YsinZd. (85) 
Therefore, 
0 
: sin? Zd 
Wis = | (DBd+w.) sinZd— V (86) 
cos Zd 
(DBd+w.) cosZd—Y sin Zd 


The vector W;, can be resolved into its S, components 
as follows: 


1 0 0 0 
Wi=T.Wi =| 0 cosZd —sinZd|| (Wia)y 
0 sinZd cosZdl! (Wia)z 
0 
= 0 (87) 
(DBd + w.) — tan ZdY 
so that 
Kik 
DBd = eS: [—(DBd + w,) + tan ZdY| 
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or finally, 


DBd =~ — w,+ Y tan Zd (88) 


which is the desired result. 

Fig. 9 shows the schematic arrangement of the gyros, 
resolvers, and servos required for these operations. The 
continuous lines show the signals and the dashed lines 
show mechanical connections. Note that the y-axis gyro 
signal goes to the y-axis servo and to the Zd resolver. 
The output of the y-axis servo is shown as driving both 
the three gyros and the Ed resolver, to which is applied 
the outputs of x- and z-axis gyros. One of the signals 
from this resolver is applied to the x-axis servo, which 
generates Zd. Zd drives the inner gimbal (via the dashed 
line to the y-axis servo, on the diagram) and the Zd re- 
solver. The other signal from the Fd resolver, together 
with the y-axis gyro signal, is applied to the Zd resolver. 
One of the outputs of this resolver is zero; the other is 
applied to the Bd servo, which drives the middle gimbal. 
Platform motion is shown as driving the outer gimbal. 


ELECTRICAL 
SIGNALS 


eee ee eR OHAL 
ROTATIONS 
KINEMATIC 
OPERATIONS 


Fig. 9—Schematic diagram of three-axis stabilization system. 


VIII. DESIGN CONSIDERATIONS FOR 
MULTIPLE-GIMBAL SYSTEMS 


The preceding sections have considered in some detail 
the analysis of two base-motion isolation systems. These 
analyses serve two purposes: one, to illustrate an ana- 
lytic procedure which is useful in determining the dy- 
namic relations which affect the performance of a com- 
plex automatic control system, and, two, to obtain 
quantitative information necessary for the rational 
design of such control systems. The purpose of this sec- 
tion is to briefly consider certain additional aspects of 
those multiple-gimbal design problems which are of in- 
terest to the automatic controls engineer. 
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In Section IV of Part 1, which was concerned with 
gimbal lock, the effect of trigonometric singularities was 
discussed in some detail. It was asserted there that these 
singularities or poles are associated with the problem of 
tracking a target near the extension in space of the outer 
gimbal axis. For example, consider that the controlled 
member is a radar antenna tracking a target; it was 
demonstrated in Section IV, that for small changes in 
either the input or the disturbance quantities, large 
changes in the output quantity are commanded, requir- 
ing unachievable angular rates and angular accelera- 
tions. As a consequence, the servo drives cannot keep 
the tracking error small enough to continue tracking, 
and the target is lost. Thus, it is said that there is a 
“hole” in the space coverage of the controlled member. 
Since the trigonometric singularities move with the 
craft, it is clear that the holes in the coverage of the 
antenna are determined by the orientation of the gimbal 
axes, the motion of the craft, and the servo rates.» The 
motion of the craft is generally beyond the control of 
engineer charged with the design of the stabilization 
system; in fact, it represents the operational environ- 
ment of the system. Further, the designer is limited in 
the torques and speeds he can expect the servos to de- 
velop, so it is necessary to consider how these holes in 
the antenna coverage may be minimized or eliminated. 
One method of alleviating this difficulty is to align the 
outer axis so that it is directed through some spot where 
the antenna is already blind. Fig. 10 shows a sketch of a 
shipboard installation,?® in which the poles associated 
with a radar antenna mount are made to coincide with 
the fantail and the superstructure of the ship. Another 
method of minimizing the holes in the antenna coverage 
is to add another axis, utilizing a four-axis mount in- 
stead of a three-axis mount. This redundancy can be 
employed to eliminate holes and reduce the servo rates 
required to stabilize the antenna, but creates new prob- 
lems (e.g., the design of a computer to determine which 
of two or more servos has control authority at any in- 
stant); systems of this type will not be considered here. 

Another consideration in the selection of a gimbal 
configuration is the relative value of the three com- 
ponents of platform motion. For example, let w,, w,, and 
w, be the x, y, and z components of the angular velocity 
of the platform with respect to the reference frame. If 
W,>>W,, W:, then it is advisable to align the outer gimbal 
axis as closely as possible with w,. Then the inner gim- 
bals are, to a large extent, isolated from the platform 
motion, reducing the velocity and acceleration require- 
ments imposed on their servo motors. This concept is 
demonstrated with a numerical example following the 
discussion of Fig. 6. 

Fig. 6 shows the six ways the three axes of a three-axis 
gimbal system can be arranged. It may be observed that 

% Cady, et al., op. cit., refer to more extensive studies of this prob- 
lem performed by H. M. James, as for example, “Train Rates in Two- 
Axis Director,” Rad. Lab., M.I.T., Cambridge, Mass., Rept. No. 8; 


September 18, 1943. 
26 R. S. Sanford, unpublished note. 
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Fig. 10—Tilted-axis two-gimbal system. 


a characteristic of these six configurations is that when 
all three gimbal angles are zero, the three axes are mu- 
tually orthogonal. It is possible to have a three-axis con- 
figuration, in which the inner and outer axes are parallel, 
when the middle gimbal angle is zero. Such a configura- 
tion may be termed a redundant system (v7z., the four- 
axis mount alluded to above) and is not considered here. 

The servo rates required to stabilize each axis of the 
six gimbal configurations shown in Fig. 6 are listed in 
Table I. These rates are expressed in terms of wz, Wy, Wz 
and the geometry which applies to each arrangement. 
It should be noted in each case by setting the inner 
gimbal angle equal to zero, the system reduces to a two- 
axis mount. For example, the three-axis mount in sys- 
tem II reduces to the two-axis mount discussed in Sec- 
tion VI, if G23=0, and the elevation and train rates 
given in Table I correspond to (38) and (39), respec- 
tively. Because there are six configurations, most of 
which have no assigned notation, it is convenient to 
denote all the gimbal angles by the letter G, followed by 
two digits: the first runs from 1 through 6 and designates 
the system; the second runs from 1 through 3 and desig- 
nates the location of the axis, with 1 being located be- 
tween the platform and the outer gimbal, 2 between the 
outer and middle gimbals, and 3 between the middle 
gimbal and the inner gimbal, or controlled member. 

The maximum servo rates associated with each axis 
of each type mount are also indicated in Table I, to- 
gether with a recommended location of the gyros re- 
quired for stabilization. The locations are determined 
so as to use single gimbal gyros for stabilization, and 
required AGC operators associated with these gyro loca- 
tions, are also tabulated. 

As an example of how Table I may be used, consider 
the case where 


40° per second 


4° per second 


4° per second 


That is, the angular rate of the platform, in platform 
coordinates, is 40° per second about the x axis (roll) and 
4° per second each about the y and z axes (pitch and 
yaw.) 
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Then, from Table I, the rates for each servo in all six 
cases may be determined as indicated in Table II. Note 
that the rates required for configurations III and IV are 
significantly smaller than those required for the other 
four types. This means that smaller servos and narrower 
servo bandwidths may be used in either the system III 
or IV, as compared with the other four configurations. 
Finally, it may be observed that if a conventional verti- 
cal gyro?’ is used in system IV (measuring pitch and roll) 
the stabilization orders from the director can be taken 
directly from the vertical gyro axes, without the use of 
a coordinate converter; this greatly simplifies the design 
of the system, and improves its accuracy and reliability. 
Consequently, the system IV gimbal configuration is 
recommended for this installation. The physical sig- 
nificance of this choice is clear from the sketch of system 
IV shown in Fig. 6. Since the outer gimbal rotates about 
the x axis of the vehicle, it “unwinds” the most violent 
component of the vehicle’s motion, thus allowing the 
two remaining gimbals to move quite slowly. The let- 
ters, O, M, and J in the lower left-hand corner of each 
box designates the order of the gimbals from the base 
out. (See Fig. 6.) The quantities a, y, and 6 are defined 
on Table I. In all cases the middle gimbal angle is taken 
as 45°, considered representative of normal operation. 


TABLE II 
SERVO RATES, DEGREES PER SECOND 
Suan Middle. | ote, Gim: 
(see Train Roll Elevation| Gimbal fee Man - 
Fig. 6) Angle cake 
I 45 41 58 Gigi 45° 
O M 16 
——— a= 6° 
II 45 58 41 G22 =45° 
O I M 
Ill 6 46 8 G32 —452 
M O if 
|  y=45° 
IV 8 46 6 G42 =45° 
i O M 
V 41 58 45 Gs2—45. 
M if O 
ae le | ee ee a 
VI 58 41 45 G62 =45° 
i M O 


It may be observed that if the controlled member is to 
be isolated from the motion of the base, its angular ac- 
celeration must be zero as well as its angular velocity. 
That is, to geometrically stabilize the controlled mem- 
ber, it is required that 


W ia =o 0, De. = 0. 
Knowing the angular acceleration for each servo, as well 
as its angular velocity, and given the dynamic accuracy 


27 .. Becker, “Gyro pickoff indications at arbitrary plane atti- 
tudes,” J. Aeronaut. Science, vol. 18, pp. 718-724; November, 1951. 

M. J. Abzug, “Application of matrix operators to the kinematics 
of eee motion,” J. Aeronaut. Sciences, vol. 23, pp. 679-684; July, 
1956. 
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requirements of the system, an estimate of the servo 
bandwidths may be determined. Further, from the 
gimbal inertias, required servo torques may be esti- 
mated. Then the designer can determine the required 
gear ratios, and make estimates of horsepower require- 
ments. Two remarks should be made concerning these 
calculations. First, the angular acceleration has two 
components—one due to the relative angular accelera- 
tions, and the other, to products of angular rates (the 
Coriolis accelerations). Some of the problems involved 
in evaluating angular acceleration are discussed in Ap- 
pendix II. Second, in determing the torque require- 
ments care should be used in applying Newton’s second 
law. For example, if the angular acceleration of the con- 
trolled member with respect to inertial space is zero, 
then the torque required to accelerate it is zero. This is 
true if a one-speed servomotor®® is used; however, in 
most cases, the controlled member is driven through a 
gear train, and the acceleration of the motor armature 
with respect to inertial space may be considerable. 

For example, assume that the controlled member is 
connected to the armature of a servo motor by a simple 
spur gear train. Then the axes of the motor, controlled 
member, etc., are all parallel, and from Newton’s second 
law, the following scalar equations may be written: 


J, DWiget eas (89) 
JuDW im = Tn — NT 3: (90) 
Jn,DW ip a i (91) 


where 


J,=moment of inertia of controlled member, 
Jm=moment of inertia of servo motor, 
J,=moment of inertia of platform, 
W;,=angular velocity of controlled member with re- 
spect to Newtonian frame, 
Win =angular velocity of servo motor with respect to 
Newtonian frame, 
Wiyp=angular velocity of platform with respect to 
Newtonian frame, 
T,:=torque applied to controlled member by gear 
train, 
T4= disturbing torque applied to controlled member, 
I‘, = torque applied to gear train by motor, 
N =gear ratio, 
DO) =di/dt (@. 


Further, the kinematic relations are 


Wim = Wip 3 Wom (92) 
Wie = W ip ot W oc (93) 
Wey Wosk (94) 


Eq. (92) states that the angular rate of the motor with 
respect to the Newtonian frame is equal to the angular 
rate of the platform with respect to the Newtonian 


28 F. M. Bailey, “Performance of drive members in feedback con- 
trol systens,” IRE Trans. oN AUTOMATIC ConTROL, vol. AC-1, p. 
74; May, 1956. 
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frame, plus the angular rate of the motor with respect 
to the platform. Eq. (93) states that the angular rate of 
the controlled member with respect to the Newtonian 
frame is equal to the angular rate of the platform with 
respect to the Newtonian frame, plus the angular rate 
of the controlled member with respect to the platform. 
Eq. (94) states that the angular rate of the controlled 
member with respect to the platform is N times the 
angular rate of the motor with respect to the platform; 
N is the gear ratio, a number much smaller than one as 
defined here. 

The servomotor functions so as to make W,, and 
DW. approach zero; that is 


Wi. 0 DW: — 0 


so that 


Win — Woe hee Ta: (95) 


Therefore, combining (90), (92), (94), and (95) gives 


4 


NT, 


T2DWe. (96) 


AT 


i 


For N<1, (96) 


ImDW ip 
phe eee (97) 
N 


and it may be observed that while J,;,DW:, may not be 
absolutely large, the factor 1/N may be considerable, so 
that the required motor torque 7; is usually not 
negligible. 


CONCLUSION OF PART 2 


The vector algebra described thus far has been shown 
suitable for the manipulation of three-dimensional posi- 
tion vectors (Part 1) which may be used to design co- 
ordinate converters, and for the treatment of three- 
dimensional angular velocities (Part 2), which may be 

_used to design stabilization systems. In the third and 
concluding part of this paper, we will extend our vector 
algebra to problems in kinetics; that is, we will examine 
the forces and moments which produce the motions 
analyzed in Parts 1 and 2. Such studies are of interest in 
the design of stabilization systems, but the primary 
emphasis in Part 3 is on gyroscopic phenomena. Gyro- 
scopic devices are of special interest since they play such 
a basic role in all of the devices described so far, and, in 
Part 3, the equations which described the performance 
of these devices will be developed in detail. 


APPENDIX I 
INSTRUMENTING THE COORDINATE CONVERTER 


In Part 1 the principles which govern the design of co- 
ordinate converters were examined. It was shown that 
the design of such devices can be accomplished by the 
use of plane geometry and implicit computation tech- 
‘niques without the necessity of explicitly solving the 


d 


J 
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pertinent equations of spherical trigonometry. The pur- 
pose of Appendix I is to describe the electrical com- 
ponents which may be used in instrumenting such com- 
puters, and to indicate some of the detailed problems 
which are involved in the design of these devices. In ad- 
dition, a proof of the implicit computation technique 
presented in Part 1 is given. 

Figs. 11 and 12 show two components which may be 
used to instrument coordinate converters.2® Fig. 11 
represents an electromagnetic resolver, capable of re- 
solving a plane vector from one set of Cartesian co- 
ordinates to another. Let @ represent the position of the 
shaft, analogous to the rotation of Cartesian coordinates 
about a common axis; further, let the vector in one set of 
coordinates be represented by the ac voltages Ex, and 
Ey,. Ex; and Ey, are the analogs of the x and y com- 
ponents of the vector in S;. The electromagnetic re- 
solver produces the two output voltages Ex. and Eye 
such that 


Ex» => 
Ey» = 


Ex; cos 6 + Ey, sin 6 
— Ex,sin 0 +-- Lyi cos 6 


so that Hx, and Ey» are the ac voltage analogs of the x 
and y components of the vector in So. 


E. =e cos 0+ Ey sin 


9 
Ey 
6 
Fy = Ex sin O+E, cos 0 
E 
ay 


Fig. 11—Electromagnetic resolver. 
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Fig. 12—Potentiometer resolver. 


Usually, when a number of such devices are cascaded 
together to form a coordinate converter, buffer ampli- 
fiers are inserted between stages to prevent adverse 
loading and coupling effects. A number of manufac- 
turers, however, have developed components which do 
not require such isolation; in fact, most of these com- 
panies provide complete resolver chains built to the 
customers’ requirements. These resolver chains include 
carefully designed mechanical structures and gear trains 
integrated into highly compact packages, operating with 
great precision over a wide temperature range. 


29 Cady, et al., op. cit., pp. 123-126, and Greenwood, et al., op. cit.; 
the latter reference is an excellent source work for problems of this 
nature and includes a good discussion on implicit computation tech- 
niques as well. 
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The device shown schematically in Fig. 12 makes use 
of two nonlinear potentiometers, whose resistances are 
so arranged that the output voltages vary sinusoidally 
with the rotation of the shaft. In this device, it should 
be noted that voltage-summing devices are required to 
properly form the two components of the output vector. 
Note, also, that the sine taps of the two potentiometers 
must be mechanically phase 180° apart. 

Potentiometers can be used with either ac or dc volt- 
ages; consequently they are frequently used in dc analog 
computers. One of the difficulties encountered in the use 
of sine-cosine potentiometers is the resolution of voltage 
from wire to wire and the action of the wiper, both of 
which result in noisy operation. In computer use, large 
diameter mandrils minimize the effects of resolution 
errors, and, in addition, the use of deposited materials 
in place of wire winding greatly minimizes the problem 
of noise. 

Fig. 13 shows a partial space-flow diagram for a co- 
ordinate converter, and Fig. 14 illustrates in block dia- 
gram form the servo driven resolver which performs the 
operation indicated by Fig. 13. The vector equation rep- 
resenting the ideal operation of the converter is given by 


1 


Ue = 


Xd 
= TarT ra} ya | ete. 
2a 
To analyze the servo performance, let Ed’ =(Ed’) 
+AEd’, where (Ed’)) is some equilibrium value, and 
AEd’ is a perturbation about that value, to which the 
usual trigonometric approximations can be applied, 
namely, 


sin AEd’ ~ AEd’, 
cos Aid’ = 1... 
Then it follows that 


(%a)o = 1 = (a)o cos (Ed’)y — (22)o — sin (Ed’)o 
(Za)o = O = (%)o sin (Ed’)o + (z)9 cos (Ed’) 
so that 
(x»)o = cos (Ed’) 
(zs)0 = — sin (Ed’)p. 


From Fig. 14, the servo relations can be written as 


(Ed’)) + AEd’ = Y(D)zq (98) 
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Fig. 13—Partial space-flow diagram for coordinate converter. 


Fig. 14—Partial servo block diagram for coordinate converter. 


where 


cos (Ed’)o sin |[(Ed’)) + AEd’| 
— sin (Ed’)o cos [(Ed’)) + AEd’| 
Y(D) = servo transfer function. 
For AFd’ small, (98) becomes 
(Ed’)o 
V(D)+1_ 


Za 


AEd! ~ (99) 


A typical form for Y(D) is given by 
Ky 

DCD ete ees 

D(TD + 1) 


so that (99) may be written as 


D(TD + 1)(Ed’)o/K, 
ang _ DED + DEA) Ky 


(100) 


D? + 


DD tel 
K, Rares 


If (Ed’)) changes at a constant rate W, the equilibrium 
value of AEd’ is given by 


(AEd’) “f 
fo K, 
Therefore, as long as W/K, is sufficiently small, so that 
the small angle approximations are valid, the resolver 
loop is stable and the usual linear design techniques for 
the selection of K, and T are applicable. 

When the small angle approximation is not used, it 
can be seen that the system behaves like a pendulum, 
with a sinusoidally varying restoring force. Ku,* among 
others, considers this problem in detail, and the inter- 
ested reader is referred to this work for a more extensive 
analysis. 

In Section III of Part 1 a coordinate converter was 
considered which was characterized by the inputs E and 
B, the disturbances Eio and Zo and the outputs Bd’ and 
Ed’. Eq. (11) stated that to solve for Bd’ and Ed’ in 
terms of FE, B, Eio, and Zo, it is necessary to solve the 
following matrix equation: 


30 'Y, H. Ku, “Analysis and Control of Non-Linear Syst ” Th 
Ronald Press Co., New York, N. Y., pp. 69 ff; 1958. ystems, e 
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Since 7a:7T14 =I, d2 =a3=0, so that (101) becomes 
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Lo asamas Now, let 

Tor =| 0 be bz (101) A, Az As 
Or Cam Ce Tar =| Bi Bz Bs 


Tat = 


and 74; is given by 


where 


The solution of (102) is simplified by noting that 


or 


To, = 


A) 
0 be 
C2 


0 
bs 


C3 


Ta: = abl tal dri sit tal ht 
cos Ed’ 0 -—sin Ed’ 
(eae 1 0 
sin Ed’ 0 cos Ed’ 
cos Bd’ sin Bd’ 0O 
Tsa = |—sin Bd’ cos Bd’ 0 
0 0 1 
1 0 0 
eee Se cosZo sin Zo 
0 —sinZo cosZo 
cos Eio 0 -—sin Zio 
dnt ae 1 0 
sin Kio O cos Eio 
cosB —sinB O 
Tin =| sin B cosB 0 
0 0) 1 
cos & Onesinwt 
La 1 0 


—sinH 0 cosE 


1 
Pat = TavT vt a 0 
0 
1 
Tot = fipr 0 
0 
cos Ed’ cz sin Ed’ 
0 be 
—sin Ed’ c.cos Ed’ 


Ja v0 
be bs 
C2 C3 
0 O 
be bs 
(7) C3 
c3 sin Ed’ 
bs 


c3 cos Ed’ 


(102) 


(103) 


Cre -C; 
so that 
Toe = TraT ai 
A,cos Bd’ + B,sin Bd’ x x 
= |—A,sin Bd’ + B,cos Bd’ « x} (104) 
Cy be ig 
where the x’s represent quantities which are of no inter- 


est here, since, by equating (103) and (104), 
Ed’ = sin-t (—C;) 


B 
Bao =tane! (=) ; 
Ay 


That is, in order to determine the output quantities Ed’ 
and Bd’ in terms of the inputs & and B, and the dis- 
turbances Eio and Zo, we need only to insert the values 
of Ai, By; and C, into (105) and (106) as indicated. 

Consequently, the evaluation of Ta: is considerably 
simplified because we seek only to find the first column 
of the matrix, made up of Ai, Bi, and C,. 

The transformation matrix Tg; may be evaluated by 
noting that 


(105) 


(106) 


Tar = Cee oa 
or 
Tit = TraTat 
A, 7 Aan 
=! B,cosZyo5 —CisinZ) “# x |. 
Bysin Zo + Cy cos Zo 4% .% 


That is, it is sufficient, for our purpose, to evaluate only 
the first column of 7;;, as follows: 


cos E cos B —sinB cosBsinE 
Tic = TinTne = | cos EZ sin B cos B sin B sin E 
—sin E 0 cos E 
Tre = TriTiz 
cos Hiocos BcosE + sin HiosinE x «x 
= cos F£ sin B 2 lait 
sin Hiocos Ecos B — sin Ecos Kio x «x 
so that 
A, = cos Eio cos B cos E + sin Eio sin E, 
B, = cos Zo cos E sin B 
+ sin Zo(sin Eio cos E cos B — sin E cos Eio), 
C,; = — sin Zo cos Esin B 


+ cos Zo(sin Eio cos E cos B — sin E cos Eto). 
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Therefore, from (105) 
Ed' = sin“ [sin Zo cos E sin B 
+ cos Zo(sin E cos Eio — sin Eio cos Ecos B)| (107) 
and from (106), 


January 


The set (xyz); represented by S; on Fig. 15 defines the 
reference space, while the set (xyz)a, represented by Sa 
on Fig. 15, designates the rotating body. S, and 5S, 
represent auxiliary frames. 


Bol = tiene” 


cos Zo cos E sin B + sin Zo (sin Eio cos E cos B — sin E cos Eio) 


(108) 


cos Lio cos B cos LE + sin Eio sin L 


Comparison of (107) with (21) in Part 1, and (108) with 
(17) in Part 1 shows that the technique of implicit solu- 
tion given in Section ITI is valid. 


APPENDIX II 


EULERIAN ANGLES AND VECTOR DIFFERENTIATION 


In Part 1, it was demonstrated that a vector could be 
transformed, via a series of planar rotations, from one 
set of three-dimensional coordinates to another. Gold- 
stein*! points out that three such rotations serve as three 
independent parameters which can be used to specify 
the attitude of a rigid body in space with respect to some 
fixed reference frame. Classically, three such angles, 
suitably chosen, are known as “Eulerian angles.” How- 
ever, as Goldstein points out,” there is no general agree- 
ment as to how the Eulerian angles are defined. Conse- 
quently, for the convenience of this paper, we have 
generalized this concept and designated as “Euler 
angles” any three successive rotations which uniquely 
and independently define the attitude of a rigid body in 
space. 

In particular, we will consider here three angles which 
are of special interest in naval and aircraft systems. 
These three angles uniquely determine the attitude of 
ship, airplane, etc., with respect to some arbitrary ref- 
erence, and, in addition, two of these angles are directly 
measurable with a standard vertical gyro.?’** These 
angles are represented by the symbols Co, Ezo, and Zo, 
where 


Co = Heading angle, *4 
Eio = Pitch angle, 


Zo 


I 


Roll angle. 


Fig. 15 shows a space-flow diagram utilizing those 
three Euler angles to define the attitude of a body in 
space. 


31 Goldstein, op. cit., pp. 107-109. 

% Thid., p. 108. 

33 See Part 3, to be published. 

34 These definitions are taken from OP-1700, except that Co, in 
OP-1700, designates direction from true north. This restriction is not 
required here and is not implied. Furthermore, as previously stated, 
the sign of these angles is given by the right-hand convection. 


() Co @ Eio © Zo (::) 
Fig. 15—Space-flow diagram for Euler angles. 
It is of interest to express the angular rate of the body 


in terms of the rates of change of the three Euler angles. 
To do this, we note that 


Wis = TarTrpW iy + Td;Wor + Wea 


where 
DZo 0 
Wi, Se 0 W.. = |DEio 
G4 
Win = 
DCo 
cos Kio O —sin Exo 
Lp = 0 1 0 
sin Kio 0 cos Eio 
1 0 0 
Lap = 40 cos Za". sin Zo 


0 —sinZo cosZo 


Performing the indicated operations and combining re- 
sults in the following three scalar equations, represent- 
ing the body axis components of the angular rate of the 
body with respect to S; as expressed in terms of the 
Euler angles Co, E1o, and Zo, 


W. = DZo — DCo sin Eio, (109) 
W, = DCo cos Eio sin Zo + DEio cos Zo, (110) 
W,. = DCo cos Eio cos Zo + DEio sin Zo. (111) 


Given the angular rate of a platform with respect to 
inertial space, as might be measured by three rate gyros 
for example, it is of interest to express the time rates of 
change of the three Euler angles Co, Eio, and Zo in 
terms of W,, W,, and W,. This may be done as follows: 
add the product of sin Zo (110) to the product of cos Zo 


= “at” > oe 


‘d 
7 


wy 
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(111) to yield (112); subtract the product of cos Zo (110) 
from the product of sin Zo (111) to yield (113); combine 
(109) and (112) to produce (114): 


W, cos Zo + W, sin Zo 


DCo = (112) 
cos E10 

DEio = W, cos Zo — W, sin Zo, (113) 

DZo = W,-+ (W.cosZo + W,sin Zo) tan EFio. (114) 


When dealing with coordinate systems which rotate 
with respect to each other, as in the cases considered in 
this paper, it is most important to exercise extreme care 
in performing differentiation and integration operations 
on vectors. For example, let V(¢) be a vector quantity 
which varies with time, and D=d/dt( ). Note that the 
operation D[V(t)] is ambiguous, since the time rate of 
change of V(#) varies in accordance with the coordinate 
frame to which this motion is referred; to make the dif- 
ferentiation operation meaningful, it is necessary to 
indicate, by means of a subscript to the operator D, for 
example, from what set of coordinates D[V(£) ] is to be 
measured. That is 


D[V()] = 
whereas 
D.[V(t)] = 
Ve DV, 
| tajaka|| Vy | | = | tafaka|| DV, 
V, DV, 


= time rate of change of V(¢) with respect to Sa, 


where i, ja, and k, are the unit vectors of S,, and Vz, 
V,, and V;, are the x, y, and z components of V(é) in Sa, 
respectively. 

Frequently, it is necessary or convenient to determine 
the time rate of change® of a vector with respect to, say, 
S,, in terms of its time rate of change with respect to 
Sp. The rule which expresses these relations is sometimes 
referred to as the Law of Coriolis, and may be written 
as*6 


DV = DV + WaXV (115) 
where 
V = any vector 
Dav = cs (V)a, 
DV = 2 (V)s, 
dt 
Was = angular rate of S; with respect to Sa. 


* 


% These statements are equally true for differentiation with re- 


to any other independent variable. 


36 Goldstein, op. cit., pp. 132-140, 


} o 
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Fig. 16 illustrates the geometric significance of (115), 
Let S, be the reference space. The rotation of S, with 
respect to S,, denoted by W,3s, is indicated in Fig. 16 by 
the vertical arrow. If V is a vector in .S;, then D;V is a 
vector along V and rotates with respect to S,. Since .S; 
is rotating with respect to S,, V changes direction with 
time when viewed from S,. This change is normal to V, 
and is given by W.XV. If V is a position vector, then 
the magnitude of the velocity of V due to W.s is given 
by | Wo |V| sin a, where a is the angle between Wz 
and (V). The direction of Was XV is normal to both Was 
and to V. Therefore, D,V consists of two components: 
one along V, proportional to the rate at which V is 
changing in length, and the other normal to V, propor- 
tional to the rate with which V is rotating with respect 
COnSa: 


=R (t) AFTER ROTATION 


"RAt = DERAt aan + At) 


R(t + At) RW) 


| Roll ¥,, | sin a@At 
(NORMAL TO R, W 


ab) 


>| 
*o 
° 


Fig. 16—Differentiation of rotating vector. 


Frequently, vectors exist which are products of a 
vector and a scalar. For example, let 


Vd) = SHU 
so that 
D.V(t) = D[SOU] 


D[S@|O@ + SOD. [TW]. 


(116) 


That is, D,[S(t)]=D[S(t)], since the rate of change of 
a scalar is independent of the reference frame from 
which this change is measured. 

As an example of operations of this kind, consider the 
angular acceleration of a body in terms of its Euler 
angles, Referring to Fig. 15, this angular acceleration 
may be written as follows: 


Aim D Wa 


where the components of Wiaa are given by (109)—(111) 
above. Expanding A iaa we have 
a ae —1 ——d 
Ly = DT aT aH a Cedl eae eee ea) 


Each of the terms of (117) may be expanded by the ap- 
plication of (115) and (116) as follows: 


Or 
= 
op) 
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Di (rs GAWD a D( ay) Wes 


0) 0 0 cos Lio O —sin Eio 0 
+ DZo| 0 -—sin Zo cos Zo 0 1 0 0 
0 —cosZo —sin Zo sin Hio O cos Hio DCo 
I 0 0 —sin Hio O —cos Kio 0 
+ DEio| 0 cos Zo sin Zo 0 0) 0 0 
0 —sin Zo cos Zo cos Lio 0 —sin Kio DCo 
1 0 0 cos Hio 0 —sin Lio 0 
+ 0 cos Zo sin Zo 0 1 0 0 
Q —sin Zo cos Zo sin Hio O cos Eio D*Co 
— D?Co sin Eio — DCoDEio cos Eio 
=| D’Co cos Eio sin Zo + DZoDCo cos Eio cos Zo — DEioDCo sin Exo sin Zo |. (118) 


D°Co cos Eio cos Zo — DZoDCo cos Eio sin Zo — DEioDCo sin Eio cos Zo 


Di(TarW or) a D(Tar)W or ata De Wae = Wis x Worl 


0 0 0 0 
= DZo| 0 —sinZo cos Zo DEio 
0 —cosZo —sin Zo 0 
1 0 0 0 — DCo sin Eio 0 
+ 10 cosZo sin Zo D? Kio | + DEio xX | DEio 
0 -—sinZo cosZo|| 0 DCo cos Eto 0 
— DEioDCo cos Eio 
=|. D? Kio cos: Zo — DEioDZo sin Z) — DEioDCo sin Eio sin Zo |. (119 
— D*FEio sin Zo — DEioDZo cos Zo) — DEioDCo sin Eio cos Zo 
D*Zo W, BDZo D*Zo 
ras mg ke eee wy 
DWra= DaWra+t+WaxXWa=| 0 |4+}WyiX} 0 |=|+DZoW, |. (120) 
0 W, 0 — DZoW, 


Combining Equations (110), (111), and (118)-(120) As = (D*Co cos Eio cos Zo — D?Eio sin Zo) 


Byes —(2DCoDEio sin Eio cos Zo 
A 
pane, : + 2DCoDZo cos Eio sin Zo 
ep eed 2 
“) + 2DEioDZo cos Zo). (123a) 
math The following observations may be made concerning 
oat (121a), (122a) and (123a). 
A, = (D?Zo — D°Co sin Eio) 1) 1 4 FeGrat homicht aueeeen 
: : n each case, the first term on the right side of the 
— (2DCoDEio cos Eio), (121a) equation is formally similar to (109)-(111), re- 
Ae = (D°Co cos Eio sin Zo + D?Eio cos Zo) spectively. 
Ho Eo DZe tos Eis eostZo 2) In each case, the second term on the right hand 
Stee Na ne 5 side of the equation has the familiar form of the 
— oDEio sin Eio sin Zo 


“Coriolis” acceleration—that is, it is proportional 
— 2DEioDZo sin Zo), (122a) to twice the product of two angular velocities. | 
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Consequently, we write (121a)—(123a) as follows: 


A, = A,+ (ACA), (121b) 
A, = A, + (ACA), (122b) 
Az; = A, + (ACA), (123b) 
where 
A, = D*Zo — D?Co sin Eio 
A, = D°Co cos Eio sin Zo + D*Eio cos Zo 


A, = D°Co cos Eio cos Zo — D*Eio sin Ao 
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and we read ACA as Angular Coriolis Acceleration, with 
the x, y, and z subscripts denoting the x, y, and z com- 
ponents, respectively so that 


(ACA), = — 2(DCoDEio cos Eio) 
(ACA), = 2(— DCoDEio sin Eio sin Zo 

+ DCoDZo cos Eto cos Zo — DEioDZo sin Zo) 
(ACA), = 2(— DCoDEio sin Eio cos Zo 

— DCoDZo cos Eio sin Zo — DEioDZo cos Zo). 


I 


Root Locus Properties and Sensitivity 


Relations in Control Systems” 
HANOCH URt+ 


Summary—the differential properties of root loci including pole 
sensitivity, angle of slope, and curvature at ordinary and irregular 
points are investigated in a unified manner. A relation between the 
sensitivity function and pole sensitivity is established. The sensi- 
tivity is shown to determine variations in the transfer function due to 
large (not only infinitesimal) variations in K. Additional properties 
of loci which are developed include loci of a variable pole position 
and the existence of asymptotes for open-loop transfer functions with 
no poles or zeros at infinity. 

The locus is treated as a transformation of a line (the real axis) 
in the K plane to the s plane, and properties of analytic functions are 
used to simplify calculations and results. It is shown that the properties 
obtained can be extended to the general root locus of a nonreal K. 


INTRODUCTION 


N the study of linear systems it is interesting to 

investigate the effects of variations of certain ele- 

ments, or parameters, on the network function. It 
has been noted by Bode! that one can write the network 
function in a bilinear form, 


ce A(s) + KB(s) 


= ——— (1) 
C(s) + KD(s) 


T(s) 


* Manuscript received by the PGAC, December 17, 1958; re 
vised manuscript received, June 23, 1959. This work was performed 
in connection with contract No, AF-18(603)-105 with the Air Force 
Office of Science Research of the Air Research and Development 
Command while the author was at the Microwave Res. Inst., Poly- 
technic Inst. of Brooklyn, Brooklyn, N. Y. 

+ Radio Corporation of America, Camden, N. J. : 
1H. W. Bode, “Network Analysis and Feedback Amplifiers De- 
sign,” D. Van Nostrand Co., Inc., New York, N. Y., p. 223; 1945. 


, 


t 


7 


where K is the variable element. A, B, C, and D are 
functions of all other fixed elements or parameters, and 
the complex variable s. In problems involving control 
systems, a commonly encountered special case occurs 
when A=0 and B=D. 

Although the transmittance 7(s) is written formally 
as a function of s alone, it should be noted that it is 
actually a function of s and K. Bode investigated the 
differential changes in JT due to small changes in K and 
introduced the term “sensitivity” commonly defined 
now as the logarithmic derivative of T with respect to 


the logarithm of K. 4 
dT 
reap nee A ae He it 
Sz = aan 2 
aie Sipedes ir 2 
K 


In recent years, with the increased utilization of the 
pole-zero approach, it has become increasingly impor- 
tant to examine variations in position of the poles and 
zeros of the network function due to changes in network 
parameters. 

This work will consider the effect that variation of K 
will have on the poles and zeros of 7, both in the “large” 
and in the “small.” The effects in the large, or the 
properties of root loci, are better known and more in- 
vestigated than the effects in the small, z.e., differential 
properties which will be discussed here in greater detail. 
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Both the numerator and denominator of 7 are of the 
form 


H(s, K) = q(s) + Kp(s) (3) 
where K is the variable parameter and p(s) and q(s) are 
polynomials. The roots of H(s, K)=0, which charac- 
terize H, therefore, depend upon the value of K. 
H(s, K) =0 constitutes an implicit relationship between 
sand K, and as K takes on all positive values, a number 
of curves in the s plane are obtained. These curves are 
the loci of the roots of H, 1.e., the “root loci.” 

The root locus was introduced by Evans? in 1948. He 
used it mainly for the investigation of control systems 
with a transfer function 7’, of the form 


Ciara K 20) 


: (4) 
1+G 4s) + Kp(s) 

where G=Kp/q is the forward gain in a unity feedback 
control system. For such systems, only the denominator 
is of interest. As a result, there arose the generally ac- 
cepted, but unfortunately confusing terminology, in 
which the zeros of p and q were called open-loop “poles” 
and “zeros,” respectively, and those of H, “closed-loop 
poles.” Root locus techniques are useful both as an 
analytical tool and as an aid to the rapid evaluation of 
control system characteristics. Rules helpful in drawing 
of root loci quickly were established.’ 

An alternate approach is obtained by viewing the 
relationship between s and K as a conformal mapping? or 
transformation. The root locus then becomes a mapping 
of the real axis of the K plane into the s plane (the defini- 
tion can be extended to other lines). The characteriza- 
tion of J in terms of its critical points (poles and zeros) 
has led to the investigation of the differential relation- 
ship between K and the s,>~® where the s; are the solu- 
tions of the equation H=0. 

The relationship, known as the pole-zero sensitivity, 
takes the place of the previously defined sensitivity. 


2W. R. Evans, “Graphical analysis of control systems” Trans. 
AIEE, vol. 67, pt. 1, pp. 547-551; 1948. 

G. Truxal, “Automatic Feedback Systems Synthesis,” Mc- 
Grav Hil Book Co., Inc., New York, N. Y., ch. 4; 1955. 

4F. M. Reza, “Some mathematical properties of root loci for 
control systems design,” Trans. ATEE, Commun. and Electronics, 
vol. 75, pt. 1, pp. 103-108; March, 1956. 

5 The roots of H=0 are denoted as s; to distinguish them from the 
variable s which may, of course, assume any value. It also empha- 
sizes the fact that H=0 may be satisfied by several values of s. 

6 A. Papoulis, “Displacement of the zeros of the impedance Z(p) 
due to incremental variations in the network elements,” Proc. IRE, 
vol. 43, pp. 79-82; January, 1955. 

Teas Truxal and I. Horowitz, “Sensitivity Considerations in 
Active Synthesis,” presented at Midwest Conf. on Circuit Theory; 


1956. 
E. J. Angelo, Jr., “Design of Feedback Systems,” Microwave 
Res. ee , Polytechnic Inst. of Brooklyn, Brooklyn, N. Y., Rept. No. 
R-449 
as & in, “Geometrical Properties of Root Loci,” M.S. thesis, 
_ Purdue University, Lafayette, Ind.; August, 1956. 
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THE RELATIONSHIP BETWEEN SENSITIVITY 
AND POLE SENSITIVITY 


In order to establish the relationship between sensi- 
tivity and pole sensitivity, let 


ell+s) 


LAS) te aie are (5) 
II w+) 


and let the z;, p; and g be functions of some parameter 
x (temperature, gain of some element, etc.). 
The logarithmic derivative of 7’ with respect to x is 


Oz Op: 
dT 1 dg Ox s Ox 
—=ds|— = +> -> (6) 
d: 20x tae Seeea ives 1 2: 
On dividing both sides by dx/x 
dT og 
dx Ox m  ~02; 1 n 
Sz? = = =r x = 
‘e g ~ Ox s+ 2; 2 
2 x 
OD: 1 
x n* . (7) 
Oz s +P: 


Eq. (7) gives the relationship between the classical defi- 
nition of sensitivity and the pole and zero sensitivities. 
It confirms the intuitive feeling that a small value of 
S;” would be achieved, for a given value of s, if the p; 
and g; in the neighborhood of s were associated with a 
low pole-zero sensitivity. The pole-zero sensitivities can 
be defined in several ways, but in light of (7), 


BEDE OD: 
z -—) xe— 
Ox 
AHN 02, 
x = a 
ax (8) 


is one apparently useful choice. 

This definition is not necessarily the most useful for 
all situations. Occasionally, an alternate definition may 
be necessary, meaningful or useful (see Appendix II). 


THE POLE SENSITIVITY 


In the euatiels of the root locus, the question of its 
behavior “in the small”—around a particular point s;— 
is of interest. This investigation starts with the examina- 
tion of ordinary points where the locus has a definite 
derivative and mapping between s; and K is conformal. 


7 “e me ° . . . 
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The root locus is defined by 


at ly 
H(s, K) = q(s) + Kp(s) = 0; (9) 
hence 
q(s) : 
K = — = — [J (s; — 5,4 10) 
6) Ls, (10) 
where A;=1 if s; is a root of gand A;= —1 if s;is a root 
of p. 
The logarithmic derivative of (10) is 
dK 
K dp Ai 
= 11 
ds; g p X Sj — Sj ( ) 


The reciprocal of (11), the pole sensitivities with respect 
to K, is 


EOF ELSON 2 a at Ca 
ak Ai qb’ — py 
K Se 
K 
gq + Kp 
it follows from (10) and (12) 
K K 
ary Saba sie Rg ee PS 
qg + Kp 0H 
Os 


where k; is the residue expressed in (4) at the point 
s;.6810 The same result can be obtained directly from 


(9) by 


oH 
ds; OK K 
fe, eae Est ieaust (14) 
dK 0H dq a) 
K Os; ds ds 


Forms (12) of S% are more convenient than others since 
they require only the knowledge of the “open-loop” 
transfer function and need not be in factored form. 
Evaluating the pole sensitivity by the residue of the 
closed transfer function will call for the complete deter- 
mination of the closed-loop transfer function, In many 
cases this will not be called for in the design as interest 
may be restricted, for example, only in the “control 
poles.” 

The slope of the root locus is that of As; (Fig. 1), but 
since dK /K is a positive number, one can find the angle 
of the slope by finding the argument of Sx 


-10R, Y, Huang, “The sensitivity of the poles of linear closed-loop 
systems,” Trans. AIEE, Applications and Industry, vol. 77, pt. 2, 
pp. 182-187; September, 1958. This paper came to the author's at- 
tention after the present paper had been written. 


d 
4 


A;sin 0; 
Se ase dw 


dy Im Sx? l; se 

dx Re Sx? A; cos 86; ae Oo 
oe : Lan 
yr OOO" 
yr: eae: Ewa 

‘(4 — ee ( bie Ne 
pi ay ey (15) 
re 


A; 
2, (x — wi)? + (y — yi)? 


where 6; and J; are the argument and magnitude of 
sj—S; (Fig. 2). It may be easier, in some cases, to com- 
pute dy/dx from the slope of ds/dK obtained by dif- 
ferentiating —K=p/q. 


GRAPHICAL EVALUATION OF Sk 


Sk, the pole sensitivity at any value s [point A 
Fig. 3(b) ], can be evaluated graphically by summation 


Fig. 1—Root locus slope. 


Fig. 2—Argument and magnitudes of (s;—s;). 


2Oi=7+Aa 


(b) 
Fig. 3—(a) Definition of AK. (b) Evaluation of Aa. 
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of angles and without evaluation of K, by the following 
procedure: Choose a point C such that AC is perpendicu- 
lar to the locus and evaluate the phasor AC =As. Evalu- 
ate the value of the phase of K at C by the summing 
angles!! and let this angle be w+Aa. S% is then As 
over jAa. 

This simple procedure is based upon the independence 
of the derivative of an analytic function on the direction 
of differentiation. Hence AK is chosen to be A’C’ which 
is mapped into AC. From the geometry one notes 
[Fig. 3(a) ] that AK =jAaK and hence 


As As As 
jha PAS AK 
4 jK K 


SENSITIVITY TO OPEN-LOoP POLES AND ZEROS 


Let a; be a root of q 


q(s) = gi(s)(s + a); (17) 


: oe ’,.\In this system F,(0, 0) =O and curvature is 
then ARE i q (4) + [X ply) 
x 
1 Puy 
as a a aie (24) 
Do, oe SEES oleae i are : a 
Oa; oH oH OH s;-+ a; 
ae ae: aa) The center of curvature is on a line normal to the 
4 3 4 curve. Because of the particular choice of coordinates a 
Upon comparing with (9) and (13) (a; is a “pole”), positive sign of p means that the center of curvature lies 
on the positive v axis and vice versa. 
Os; 1 K is expressed as a function of the variable z=u+j 
ead aye oe : (19) and Fis defined as * 
0a; Sj + ay -_ 
A 
Similarly if @; is a root of p (“zero”) F = Im log K. (25) 
Os; 1 (Sometimes it is more convenient to define F=Im K.) 
Pee = Sk? Eas j (20) As K takes on positive real values, F is equal to zero and 
E : ‘ therefore serves as the real function of (23). 
THE TOTAL SENSITIVITY 1 i (Im log K)uu Im [(log K) ux] od 
0s; dK OS; p F, (Im log K)» Im [og K)e] 
$7 = SA) | da;. (21) 
va is Oa. But for an analytic function (such as log K) 
K 
) d 
dK da; eee 
ds; = Sx?| — — >) A; | (22) ou dz 
iG Sp Ge 
0 0 d 
See Sig ta (27) 
THE CURVATURE OF THE Locus AT ORDINARY POINTS dv Ojv dz 
The curvature of a curve, which is defined by means it pees 
of an implicit relationship F(x, y) =0 is? 
1 Im (log K) 22 Im (log K) zz (28) 
1 Truxal, op. cit., p. 226. f Imj (log K). Re (log K), 
2 R. Courant, “Differential and Integral Calculus,” Interscience Be 


Publishers, Inc., New York, N. Y., vol. 2, p. 125; 1952. 
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1 Fale = 2h aly Layla 


23 
p ed PR) a 


The subscripts denote partial differentiation with re- 
spect to the indicated variable. 


It is convenient to wtrornce & new ar ceen a 
with its origin at s; and oriented so that the tangent to 


the locus coincides with the uw axis (Fig. 4). 


Fig. 4—u, v coordinate system. 


1960 
sin 2¢; 
iad: : 
12 
p cos ¢, 
ra" 


2(u — u:)(v — 2) 
"(u — us)? + (v — 0)? 


Pe 


“yea (29) 
uU— Uj; 
A; 
u (a — u:)? + (v — 0,)2 
where 
$5=0;—ea (30) 


and a@ is the angle between the wz axis and the «x axis. 

The transformation from s to gz need not be carried 
out by evaluating (28). Consider instead the general 
linear transformation 


z= (s — soe. (31) 
Then, for any function f, 
d df ds d 
Pete Os, (32) 
Csea COSMO zmeds 
a? ep yas\ |. af a's... d* 
BEER ae cy 
dz? ds? \ dz ds dz? ds” 


However, since all orders of derivatives of s with respect 
to z higher than the first are zero, 


ay d"f (=) 
dz” ve ds” \ dz i 


ay, 


mp cael 


ds” 


(34) 


Examples are given in Appendix I. 


PROPERTIES OF THE Root LOCUS AT 
IRREGULAR POINTS 


The previous discussion is for the case where the locus 
has a definite derivative. Such points are characterized 


by 


(35) 


At points where these restrictions do not hold, the root 
locus has no definite derivative and a multiple point 
occurs. Condition (35) implies that S% does not exist 
at an irregular point. 

As an illustration, consider a point where (log K)’=0 
on the real axis. If both p and gq have real coefficients 
then Im(log K)’=0 everywhere on the real axis and 
only Re(log K)’=0 need be investigated. 


N= hy 


Rees se Ee 0h (36 
(x — xi)? + (y — 94)? gO 


Re (log K)’ = >> 
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If this summation is split into a term ).>’ corresponding 
to critical points on the real axis and a term ).” cor- 
responding to critical points in the upper half plane, the 
following expression is obtained: 


1 2(% = Xi) 
O= DY A +" A; 
(@—#)? + ye 


The factor 2 in }>”” is necessary to account for the sym- 
metry about the real axis. Eq. (37) is the well-known 
expression for the “breakaway” points. Another example 
appears in Appendix I. 


(37) 
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THE SENSITIVITY AND SLOPE AT 
IRREGULAR POINTS 


Let 
K JK MSC : es (38) 
— Ay eo ee 
: ds 2 ase j 
and let the first 7 —1 derivatives vanish. Then 
1 dK 
AK = K — Ky =— As” (39) 
n! ds” ; 
Noe lv Ke 
or ie Dee 
we AK n n!} Praha: 
As = ie = INI 40 
/ * bf ay (40) 
ds” ds” 


As As/XYAK is finite it may be substituted for Sx? at 
irregular points. 

For any AK>0, ds has n values separated by 27/n 
which represent the outgoing branches in Fig. 5. 2 
branches with the same separation as before are again 
obtained if AK <0 but they are displaced from the pre- 
vious set by 27/2n, that is, the two sets intersect (Fig. 
5). Eq. (40) gives the slope at multiple points. This re- 
sult is of a different form than that obtained by 
Papoulis.® 


(a) (b) 
Fig. 5—n outgoing branches. (a) »=3. (b) n=2. 


Tue CURVATURE OF ONE OF THE Locus BRANCHES 
AT A DOUBLE POINT 


A set of coordinates uw, v is introduced as before (Fig. 
4) at the double point with each of the branches tangent 
to an axis. In the vicinity of a double point, F can be 
factored into two functions each describing one of the 
branches and having no singularity. 
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F = m(u, v)M(u, v) = 0 (41) 
because of choice of coordinates 
m(0, 0) = m,(0, 0) = M(0, 0) = M,(0, 0) = 0. 
The branch tangent to the wu axis is described by 
F 
M = — (42) 
m 


1 1 1 
Mee 7(—) + 2F, (—) + Pusu (-). (43) 
MS uu Mu ™m™ 


This expression is indeterminate at ~=v=0. After ap- 
plication of l’Hospital’s Theorem, 


1 
M,,,(0, 0) = 3 F (0, 0). (44) 
Similarly, 17,(0, 0) = F.»(0, 0). Hence the curvature is 


eee (45) 


Evaluation of the derivatives of F is carried out as be- 
fore and examples appear again in Appendix I. 


On CURVATURE, POLE SENSITIVITY AND 
THEIR RELATIONSHIP 


At first glance it may seem that the curvature and 
sensitivity are somehow related. Example 1, Appendix I, 
would reveal, however, that the curvature in that exam- 
ple is a constant, while the pole sensitivity is (in absolute 
value) equal to cot ¢—a continuous variable. The in- 
dependent nature of the curvature and the pole sensi- 
tivity can be understood by noting that the pole sensi- 
tivity depends on first derivatives while the curvature 
depends also on the second derivative and, in general, 
the two derivatives are independent. 

From another point of view, these two quantities are 
of an entirely different nature: drawing on mechanics, 
the curvature can be considered a “kinematic” property, 
while the pole sensitivity is a dynamic property. 


EXTENSION TO GENERALIZED Root Loci 


Conventionally, it is assumed that K takes on only 
positive or negative values. This case is the one treated 
in detail in this paper. The results, however, can be ex- 
tended easily to cases where the root locus is a mapping 
of some other line. Most common are the loci which are 
mappings either of the radii or of the concentric circles. 
The radii correspond to values of K with constant phase. 
The circles correspond to values of K for which |K | 
=constant. The mappings of these curves are loci which 
correspond to constant real or constant imaginary parts 
of log K. F is then’derived from log K and becomes 


F =Imlog K —y=0 (46a) 
F = Re log K — log Ky = 0. (46b) 
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Since only the derivatives of Ff appear in the formulas: 
derived and since the derivative of a constant is zero, , 
all results obtained remain valid. 

If K varies in such a way that either its real or imagi- 
nary parts are constant, the loci of interest are mapping 
of lines in Fig. 6(b). Fis then derived from K itself. 

The most common root locus corresponds to the real 
axis which is a grid line in both Fig. 6(a) and 6(b). 
Hence, F can be derived for this case from either K or) 
log K as is done in the examples. 


SENSITIVITY CAUSED BY LARGE VARIATIONS 
IN PARAMETERS 


Consider the transmission 7 of the graph shown in 
Pig yie 


ack 
T = —— = [7(s, k). (47) 
1 — Kb 
Let K=Ko+AK and T(Ko) =7». Then, 
T i} (48) 
Ti AK 
K 
1 — Kob 
But, 
i Spe ) (49) 
fESER Eee ee 
Hence, 
eo} AK n 
T =\To 4) (= s) (50) 
0 VK 


eran. 


(a) (b) 


Fig. 6—(a) Loci for K of constant magnitude and phase. (b) 
Loci for K with constant real and imaginary parts. 


Fig. 7—Typical flow graph with one feedback loop, 
no direct transmission. 
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or 


ik Peg IN GaN. 
In — = — Sey: (51) 
To i 942 IK 
These formulas are interesting since they show that as 
a result of the special form of 7, the sensitivity deter- 
mines all coefficients of the power series expansion of T. 


THE Root Locus WHEN ONE OF THE CRITICAL 
Points BECOMES VARIABLE 


In many cases, the position of one of the critical 
points (a pole or a zero) is either not known exactly or 
is a variable of design. Thus, it may be of advantage to 
know the locus of the poles when this movement is con- 
sidered. With notations similar to that of (17), one can 
show that the locus is described by the roots of 


q(s) + Kp.(s)(s + a; + Aa,) 


= q(s) + Kp(s) + Aa:Kp,(s). (52) 


K here is considered constant and Aa;,, the variable. 
Eq. (52) indicates that the new poles are the poles of 
the closed-loop system while the zeros are those of p;— 
in other words, the zeros of the original system less the 
movable zero. As an illustration, consider the second- 
order system of Fig. 8. The loop is closed and a certain 
gain of value K is chosen so that the closed-loop poles 
have the positions A and B. As there was only one zero 
in the original system, the a locus is that of a system 
with two poles at A and B and no zeros. The locus is 
indicated by the dashed Jine and the real axis in Fig. 8. 
A similar rule applies to a variable open-loop pole. 

As another example, consider a commonly used trans- 
fer function [Fig. 9(a) | 


K 


a 53 
s(s + a) a 


around which the loop is closed and the selected gain K 
corresponds to two closed-loop poles at A and B. The 
position of the closed loop poles for different values of 
a@ is given according to the rules by the locus shown in 
Fig. 9(b) by a dashed line. According to the rules given 
before, there are “poles” at A and B and there is a zero 
at C, and the “gain” is +Aa. 


Tue Locus WHEN Two ConjJuGATE PoLEs MOVE 


One has now to consider the roots of 
(s? + b;s + Ab;s to; + Aci) qi(s) + K p(s) 


= g(s) + Kp(s) + Abisgi(s) + Acigi(s) = 0. (54) 


Two cases will be considered: 
i) The two poles move on a vertical line and remain 
conjugate.'® In this case, Ab;=0 and the locus has 


13 This requirement is equivalent to stating that the sum of the 
_two poles remains constant. 
ed 
» 
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“poles” at the poles of the closed loop and the “zeros” 
are the open-loop poles which do not move. 

2) The two conjugate poles move on a circle centered 
at the origin, as in Fig. 10. This assumes that the prod- 
uct of their values is constant and hence Ac;=0, only 
Ab; varies. Examination of (54) shows that while the 
“poles” are the same as before, the zeros are given by q; 
and, in addition, there is a zero at the origin. By change 
of variable, it is easy to prove that if the variation is on 
a circle centered anywhere on the real axis, not just at 
the origin, the added zero has to be located at its center. 


EQUIVALENT LocI 


Theorem: Vhe root locus remains unchanged, if, in- 
stead of the original poles, one uses any set of points on 
the locus corresponding to any desired K=K;,. In this 
case, K has to be replaced by K—K,. 

Proof: The root locus is characterized by the zeros of 
(9) but 


qt Kp = (¢+ Kip) + (K — Kip = 0. (55) 


The right member describes the new locus mentioned in 
the theorem. 

Theorem: Interchange of open-loop poles and zeros 
would create an identical locus, which, however, has an 
opposite sense and in which K is replaced by 1/K. In 


(b) 


Fig. 9—(a) Root locus for system with two real poles. (b) Position of 
poles for various values of a. 


Fig. 10—Locus change for variation of conjugate poles on a circle. 
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the normal case of excess of poles over zeros, the loci ap- 
proach the asymptotes when K (the new one) is very 
small. For proof, divide (9) throughout by K. 


Loct witH No OpEn-Loorp POLES 
OR ZEROS AT INFINITY 


If the order of p(s) and q(s) is the same, then for 
K=-—1, the order of (9) would be reduced (assuming 
unity to be the leading coefficient in and q). Hence 
some roots will disappear or become points at infinity. 
Each point at infinity has one entering branch and one 
outgoing branch. The statement follows immediately, if 
in (9) —1+a is substituted for K: 


H(s, a) = (q— p) tap = 0. 


Here # is of higher order than g—# and hence the points 
at infinity occur for a—0. Let the difference in order be 
n. There would be x asymptotes as a—>-+0, and another 
set of m asymptotes as a——0O. The two sets intersect 
each other. (This is the same as in Fig. 5, except that the 
multiple point is now at infinity.) The rules for centroid, 
number of asymptotes, etc., are the usual ones, with 
q(s) — p(s) replacing qg. Fig. 11 describes two such cases. 


(56) 


AK=-I+a 


(a) (b) 


Fig. 11—Root loci for a system with an equal 
number of poles and zeros. 


CONCLUSION 


The usefulness of any analytical technique depends, 
to a large extent, on the familiarity of the user with that 
technique. This familiarity readily enables the user to 
extract more information about the properties of his 
problem, and guides him to the correct path of analysis 
without time-consuming false starts. Such familiarity 
requires the understanding of general theorems applica- 
ble to all problems where the technique applies. It is felt, 
therefore, that the utility of root locus analysis will be 
enhanced by the results given in this paper. These re- 
sults can be classified in several ways: properties in the 
“large” and “small.” The properties in the large which 
appear in the latter part of the paper are directed toward 
some particular cases. The case when there are no poles 
or zeros at infinity may occur quite often, yet the basic 
known textbook rules regarding asymptotes would be 
misleading here. The ability to obtain the critical points 
of a system in which an open-loop pole or pole-pair is 
varying is felt to be a useful design tool. The properties 
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in the small, such as slope, curvature, and location of 
multiple points can be very helpful in drawing the locus 
as in Example 2, Appendix I. 

The rest of the paper considers sensitivity relations 
which should be of particular value in design of practical 
systems where elements are either unknown or variable, 
and these tolerances have to be considered in trying to 
meet specifications. These ideas, and the deliberate vari- 
ation of elements to vary properties of systems in a de- 
sired way to make them adaptive, can be made a basis 
for new design procedures. 


APPENDIX I 
Example 1 (Fig. 12) 


ie i ee ee (57) 
$ S 
dK 1 a’K 2 
VES Gh STEY oe 
a?K 6 
dst st 8) 
Fig. 12—System discussed in Example 1. 
Double points occur at dK /ds=0 or s= +1. 
dK 
saa Pa 2je%% sin @ = 2 sin e—?(@tt/2) (59) 


for K <2 or on the circular portion of the root locus. 


dK 
dy ds ( ty é6 
— = ————_ =— cot¢=t ——). 
da dK oa are a Sy 
1S — 
ds 
The pole sensitivity will be 
ds s?.-- 1 2 74+1 
Sf Caine Sue gga 
dK s 1— s? st— 1 
1 
Sa 
S 
= 5 
' (61) 
iS! —_ -—— 
i) 
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and on the circular portion 


D ei? a eit 
Se aa, Sak Se oe?) 
K=e Saas ge'* cot d. (62) 
The curvature is, using (34), 
1 (Ps 8 Im (— 2e-%%¢) ef 2-7) 
=. =— = — 1, (63 
p i Im (— 2 sin ge*#e7(¢-7/2)) Sy 
The curvature at s= —1(¢=7) is 
1 a°K 
raken Im e13(r/2) 
1 3 dz’ 1 (—1)4 
=~ = =—1 (64) 
p ha 3 —2 
Imj Im (i J ener 
dz? ~ (—1)3 
The pole sensitivity at s= —1 (a double point) is infinite 


and hence (40) is evaluated. Here n=2 and, according 
to (40), 


AS... A 2 : 
TE // ayaa (65) 
3 


All these values can be verified otherwise. 

The graphical procedure for obtaining Sx? is tried at 
@ = 135°(Fig. 13). AC is chosen to be 0.20 and one finds 
that 


> 6: = 7 —- 12 —=r-—0.21 (66) 
Fig. 13—Graphical procedure for finding Sx. 
Hence 
As —0.20¢°:7577 
Sx? = — = — ———_ = 0.95 fei? 157, (67) 
jha 40.21 


According to (62) the exact value is 1 instead of 0.95. 


Example 2 (Fig. 14) 
1 
K= 
(s + 2)(s — 2)(s? + 16) 


It will be more convenient to consider w= —1/x which 
‘corresponds to the same root locus. 


(68) 
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Fig. 14—System discussed in Example 2. 


uw = s*+t 125? — 64 (69) 
du 
— = 453+ 24s = 4s(s? + 6) 
ds 
oe 12s? + 24 
he 10% 
ds? 
EN (70) 
— = 24s. 
ds* 


Multiple points occur at du/ds=0 which yield s=0, 
s=+j/6. These values of s correspond to real values 
of w and hence are double points of the root locus. The 
curvature at s=j/6 is 


1 1 


Im 24/6 
p 3 Imj(-12X6+24) V6 
ds/dz=1 in this case. 


(71) 


APPENDIX II 
On the Meaning of Sensitivity Functions 


The accepted definition of sensitivity utilizing log- 
arithmic derivatives might be ambiguous in some cases, 
especially when JT is small or zero. For a balanced bridge 
circuit, the sensitivity (with respect to any element) of 
Eo/Ein is infinite while an infinitesimal change in an 
element value corresponds to an infinitesimal change 
in Jord, 

Another example involves the sensitivity of a circuit 
to a resistance R. R is composed as a series combination 
of two resistances, R,; and Re, where R,=10R2. Then 
Sh, =11SR, while an ordinary derivative would make 
the total sensitivity equal to that of the sum. 
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Time Lag Systems—A Bibliography 


N. H. CHOKSYt 


Summary—In a recent issue of these TRANSACTIONS, Weiss! 
has given an excellent annotated bibliography on the subject of 
transportation lag. He was kind enough to refer to a previous 
bibliography which appeared in this author’s thesis (his refer- 
ence number [Ch5]). Since the latter is not generally available, the 
bibliography therein is given here. The items which have already 
appeared in Weiss’ listing have been omitted; items which have 
appeared since the thesis was written have been included. 

While most of the items here were found by library searches, 
acknowledgment must be made to Bellman’s? bibliographies as a 
source. 

If there are items which have not been mentioned either here 
or in Weiss’ paper, the author of the present paper would appreciate 
being informed about it by either a complete reference or a copy 
of the paper (or papers) in question. 

A short introduction on the mathematical characterization of 
time lag systems is given before the bibliography. 


INTRODUCTION 


INCE transient analysis is often of interest, we shall 
S consider here only the homogeneous forms of the 
governing equations. 

Linear physical systems’ are characterized by linear 
differential equations or by a set of such equations (and, 
as is well known, one nth order differential equation can 
always be written as a set of » first-order differential 
equations). Here, for convenience, we shall consider one 
equation. This, in general, has the form: 


apy (t) + avy Y(t) + + + + + any(t) 


= Day =0 (1) 
i=0 
where the a; are real or complex-valued constants, or 
may be functions of time. 

Associated with such an equation there is the so-called 
characteristic equation which is a polynomial equation 
in a complex variable if the a; are constants, or, indeed, 
may itself be a differential equation if the a; are poly- 
nomial functions of time. Again, here we shall restrict 
ourselves to the former case. Then, the characteristic 
equation can be written as: 


aos” + as" 1 +---+a, = dD ier), (2) 


i=0 


The study of the behavior of the physical system by an 
analysis of this equation fills the literature. 


* Manuscript received by the PGAC, August 17, 1959. 

+ The Johns Hopkins University, Baltimore, Md. 

1R. Weiss, “Transportation lag—an annotated bibliography,” 
ane TRANS. ON AUTOMATIC CONTROL, vol. AC-4, pp. 56-68; May, 

2 Tbid., references [Be2], [Be3], [Be4]. 

3 The word “physical” is used here in the widest possible sense to 
include all phenomena in the real world which may be characterized 
mathematically. 


Now, a differential equation as expressed in (1) tacitly 
implies that the future behavior of the system is fully 
determined if we have complete knowledge of its state, 
i.e., of the values of y(é), (¢=0, 1,---, ) at any 
given time ¢. In many systems, particularly in the fields 
of economics, population studies, biophysics, physi- 
ology, etc., this implication is not justifiable. In such 
systems, the future behavior is dependent not only on 
its present state, but also in some measure on its im- 
mediate past, and, indeed, in some cases on its im- 
mediate future. 

The author has always preferred to use the term 
“transportation lag” only when the dependency on the 
past (or future) is through the dependent variable only 
and not through any of its derivatives. The words “time 
lag” are reserved for the more general case. 

In the transportation lag case, if no derivatives are 
present, the corresponding equation becomes a dif- 
ference equation. However, of more interest is the time 
lag case where the corresponding equation is called a 
difference-differential equation. Systems characterized by 
such equations are variously called time-lag systems, 
retarded systems, hereditary systems, etc. These equa- 
tions have the form: 


n 


m 
YL ay + LV by" C-—n)=0 (3) 
7=0 k=0 
where a; and 0; are constants, tT; are real constants (posi- 
tive for dependency on the past, negative for de- 
pendency on the future) or may, in complicated cases, 
be functions of time (we shall not consider this latter 
possibility here), and m may be greater than, equal to, 
or less than a. If for a certain k, 7,=0, then the corre- 
sponding b,=0. 

The characteristic equation in such cases is no longer 
a polynomial equation, but a transcendental equation 
of the sort usually termed exponential polynomial 
equations: 


n ™m 
Dd assn-? + DY dyes” = 0. (4) 
i=0 k=0 
This equation plays a part, analogous to that of (2), in 
the study of time lag systems. 

This brief introduction can be closed by simply stat- 
ing that instead of difference-differential equations, 
time lag systems are also characterized by integral equa- 
tions of the form: 


ro aol! f nt ae ee 


This is called the renewal equation. 
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The following bibliography lists papers concerning 
time lag systems in a variety of ways: some papers dis- 
cuss physical systems and show that a difference- 
differential equation characterizes them; some go further 
and try to obtain solutions to such equations directly 
in the time domain; some examine the solutions for 
asymptotic behavior; some deal with the stability of the 
systems through the study of their characteristic ex- 
ponential polynomial equations, and so on. As the titles 
indicate, a large assortment of problems are’ considered, 
in the hope of covering all possible fields of interest. 
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Correspondence 


Properties of Root Locus 
Asymptotes* 


As a function of K, the roots of the 
rational polynomial P(s)+KZ(s) describe 
what are known as root loci in the complex 
s plane. The behavior of these roots is im- 
portant in understanding the behavior of 
linear systems as a function of the variable 
K. The variable K may be either the gain of 
some part of the system or some control 
parameter in the system. 

Based on the location of the roots of 
P(s) and Z(s) in the s plane, the root loci of 
P(s)+KZ(s) can be sketched by applying a 
number of standard rules.! A new rule is 
presented here which indicates the final side 
from which the locus approaches its asymp- 
tote as K and |s| become very large. 

The roots of 


P(s) = II @—s,,) 
ia 
are known as poles, whereas the roots of 


Zi Ve 


t=l 


are known as zeros. It is well known that, 
as K goes from zero to large positive num- 
bers, the roots proceed along the loci from 
the poles to the zeros. If there are more poles 
than zeros (n>m) in the finite s plane, the 
loci approach (n—m) asymptotes radiating 
from the centroid of the finite pole-zero pat- 
tern. The centroid is located at 


n ™ 
ee es 
t=1 


t=1 


, (1) 


C= 
n—-m™m 


and the angles of the asymptotes are located 
at 


Masai 


n—™ 


ok = a; k=0,1,°---,n—m —1. (2) 


The side from which the locus approaches 
the asymptote is determined as follows: 

1) The locus lies exactly along the 
asymptote if the pole-zero pattern is sym- 
metric about the asymptote line extended 
through the centroid. An analytical test for 
this is the substitution of 


s = s'ei*y (3) 


into the polynomials P(s) and Z(s). The 
locus lies along the asymptote if, and only 
if, the coefficients of the new polynomials in 
s’ are real. If this test fails, 
2) Find the rational 
P(s) and Z(s) in the form 


polynomial of 


* Received by the PGAC, July 20, 1959. This 
paper (Jet Propulsion Laboratory External Publica- 
tion No. 587) presents the results of one phase of re- 
search carried out at the Jet Propulsion Lab., Calif. 
Inst. Tech. under Contract No. DA-04-495-Ord 18, 
sponsored by the Dept. of the Army, Ordnance Corps, 
E LJ. F. Truxal, “Control System Synthesis, ch. 
4, “Root-Locus Methods,” McGraw-Hill Book Co., 
Inc., New York, N. Y.; 1955. 
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Z(s) = (s—c)™ 

+ a(s— P14 set ay 
P(s) = (s —c)™ 

cig O1S GIRS tein Orne (4) 


where ¢ is the centroid. 
The coefficients aq, bg may be found di- 
rectly from the pole-zero configuration by 


m—I+1 m—o+2 


GK ad)t Dine), 


tel j=tl k=j+l 


“(Su = Mba — 6) +++, — 2). (6) 


An identical expression gives bg, except that 
bz uses the poles where (5) has zeros. These 
coefficients may be found from the coeffi- 
cients of s rather than of (s—c) by 


‘ S : 
LEDD Pe Wd saerens 
7=0 m—p 
where ; is the coefficient of s”~7 in Z(s). The 
expression for bg is the same, with &; re- 
placed by the coefficients of s”~7 in P(s). 
3) Determine the first nonzero number 


p=1,2,3--- (7) 
4) Divide the s plane into 2p equal sec- 


tors about the centroid starting with the 
positive real-axis, as in Fig. 1. 


dp = bp — Ap; 


s- PLANE 


Ist 


—po 


Fig. 1—Division of the s plane into 2p sectors. If 
asymptote lies in odd sector; asymptote ap- 
proaches from right, if dp <0; from left, if dp>0. 
If asymptote lies in even sector: asymptote ap- 
proaches from left, if dp <0; from right, if dp >0 


5) If the asymptote lies in an odd sector, 
the locus approaches from the right if 
d»<0O and from the left if d,>0. If the 
asymptote lies in an even sector, the oppo- 
site is true: the locus approaches from the 
right if d,>0 and from the left if dp <0. i 

6) If the asymptote falls on a sector line, 
determine the nonzero number 

q-1 


dq = bq — % — ms Digi (8) 
t=p 
of least q such that the asymptote does not 
lie on one of the new lines which divide the 
s plane into 2g sectors. Apply steps 4) and 
5), with p replaced by q. 
In all cases, d; will be zero. The constants 
a. and by can be obtained either from the 
polynomial or from the pole-zero pattern by 
use of the formula 


_ m=1 mM 1 m 
a, = yD, oS O%2 02; TP 2 ra, Be,” (9) 
il jai+l ; jal 
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where 
az, = Re (sz, — ¢) 
Bz, = Im (s;,). 


The formula for be is identical to that for a», 
except that the b: formula uses the poles. 
Similar expressions exist for ap and by for p 
greater than 2. The numbers ag and bz may 
also be obtained from the coefficients of s 
rather than (s—c): 


a2 = ke + (m — Adie + m(m — 1)c, (10) 


where &; is the coefficient of s™~7 of Z(s) and 
c is the centroid. An identical expression 
exists for be, except that k; is the coefficient 
of s"7 of P(s). 

Two examples are given to illustrate the 
use of the asymptote rule: 

1) Consider the polynomial 


P(s) + KZ(s) = s+ 753 + 125? + K(s + 1) 
= s*(s + 4)(s +3) + K(s + 1). 


The root locus of this polynomial has a cen- 
troid at 


4/23) (21 
eS te = 2 +70 


There are three asymptotes at angles of 60°, 
180°, and 300°, as shown in Fig. 2. From the 


Lae 
s*(s+4)(s+3) 


5- PLANE 


Fig. 2—Root locus for Example 1. 


polynomial written in the form 
(ie 2) oS 2) — Os 2)- 
+4(s+2) +84 K[(s+ 2) - 1] 
it is determined that 
a, = 0 
bp = — 6. 


Since b2<a2 the locus approaches the 60° 
asymptote from the right and the 300° 
asymptote from the left. An identical result 
is obtained for ag and bz by using (9) directly 
on the pole-zero pattern. 
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2) Consider the polynomial 
P(s) + KZ(s) 
= 5 + 454 + 653 + 452 + K(s? + 4s + 8) 
= $%(52-b 25 2)(5 4-2). 4 Ks? + 45 +8), 


The root locus of this polynomial has a 
centroid at the origin and three asymptotes 
at 60°, 180°, and 300°, as shown in Fig. 3. 


jw 


a? 


_K(s*+45+8) 


=-| 
s*(s+2Ns*+2542) 


5-PLANE 


NOTE THAT LOCUS CROSSES 
ASYMPTOTE 


Fig. 3—Root locus for Example 2. Note that 
locus crosses asymptote. 


From this polynomial, 
a = 8 
bs =— (oe 


The 60° asymptote is approached from 
the right since a2>be. Note that the locus 
crosses the asymptote, a fact which might 
not be evident without an actual detailed 
plot of the root locus. 

The asymptotic behavior of the root 
locus is obtained by locating the root of the 
polynomial P(s)+KZ(s) for large values of 
K and |s|. Equating the polynomial 
P(s)+KZ(s) to zero produces the relation 
that 


—K= aa = (s se c)n-m + di(s = c)n-m-1 
a ae (11) 


The coefficients d, are 
dy = by —- a= 0 


q-1 


dg = bg — Aq — 2, Ady; (12) 
i=l 


where b,=0 for g>n, ag=0 for g>m. 
From the identity 


@-orm=-K[ | as) 


extraction of the (7—m)th root, and trans- 
position of the ¢ term to the right-hand side, 


the following equation is obtained: 
S=C + Kil Gm) i [(2k+1)/ (nm) Jr 


IK —1/(n—m) 
[- | : (14) 
s—oc)rm™ 


Substitution of (11) into (14) and applica- 
tion of the binomial theorem allow s to be 
written 


rater 
(n — m) 
[als — o)-* + dls —7+--:] 
1 i+” —m) se 

2! (n—m)?~ Cee 


s=cet Kil rm) of [@k+1/ (nm) Jr \1 a 


+ 


false Pte tA) 


For very large |s—c|, (145) shows that 
= —(s—c)r™., (16) 


Under such conditions, (15) may be ap- 
proximated by keeping only the first two 
terms within the braces. For large |s—c], 
the inner brackets become d,(s —c)~?, where 
d, is the first nonzero coefficient. Then 


Ss = 6+ KU) ei [2k+1)/(n—m) |r 


BE ices ication icel, (17) 
n— m 
Note that the first two terms of (17) consti- 
tute the vector from the origin to points 
along the asymptotes, and the third term is 
a first-order correction term possibly dis- 
placing the root off the asymptote. 
The asymptote vector is 


c + Km) gf [(2k+ Di (n—m) Jr =C¢ +. pei?k, (18a) 


while the first-order correction is 


KG-p)/(n-m) gti [(2k+1)/ (n-m) Jt (1p) 


n—m 


= dei#-»),  (18b) 


where ¢; is the kth asymptote angle. This 
approximation for the root of the polynomial 
represents the offset due to the centroid, the 
asymptote, and the lowest-order correction 
to the asymptote. Note that there are (x —m) 
asymptotes at angles 7(2k+1)/(m—m), that 
the asymptote grows as K’@™) and that 
the correction decreases as K(~p)/(»—m), The 
side from which the locus approaches the 
asymptote can be determined by considering 
the angle of the correction in relation to the 
angle of the asymptote. Substitution of 
(18a) and (18b) into (17) gives 


s= ct pele — AdyeiPe-P), (19) 


Divide the complex plane into 2p sectors by 
lines through the centroid ¢ at angles of rx/p, 
and consider points about an asymptote. 
If ¢, is an odd sector, that is, 
Qrr 2r+1 
ae < ok (“= T; 


then 
ge — 2rx > dk(1 — p) > be — (2r + 1)m. (20) 


Obviously, in these regions, for —\d,>0 the 
locus must lie to the right of the asymptote. 
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If —Adp<0, $4 (1 —P) isincreased by 7 radians 
and the locus lies to the left. Examination 
shows that this condition reverses at multi- 
ples of 1/p radians. If equal sectors of angle 
a/p are constructed in the s plane about the 
centroid, the angle condition is equivalent 
to determining the sector in which the 
asymptote lies. This argument then produces 
the first part of the rule stated here. 

If the lowest-order correction lies along 
the asymptote, it will not give an indication 
of the approach side. Likewise, powers of this 
correction term will not produce an indica- 
tion. This case occurs when the asymptote 
being considered lies on a sector line. 

The next lowest-order correction results 
from the next nonzero dy. The angle condi- 
tion requires a new set of sectors to be con- 
structed in the s plane. If the asymptote 
still lies on a sector line, the next nonzero dp 
must be considered, and so forth, until one 
is found which produces a sector which in- 
cludes the asymptote. Products and cross- 
products of the lower-order d,’s and the one 
under consideration produce either correc- 
tions along the asymptote or higher-order 
corrections than the one being considered. 
Thus, these powers and cross-powers may be 
neglected. 

If all nonzero d; (i=pe2-+ ©) require 
sector angles which lie along a certain 
asymptote, then from (15) there can be no 
displacement from any term, and thus the 
locus must coincide with the asymptote. 
This condition is fulfilled if the pole-zero 
pattern is symmetric about the asymptote 
line through the centroid. This then com- 
pletes the rest of the rule. 

Most pole-zero patterns are such that the 
loci in the upper half-plane all converge to 
the asymptotes from the upper side or all 
from the lower side. This case occurs when 
the coefficients a2 and be are nonequal. 

The coefficients a2 and by can be deter- 
mined from Z(s) and P(s) expanded about 
the centroid c: 


Z(s —c) = (s—c)™+ a,(s —c)™ 

“Gals OC) Ate ene 

P(s = c) = (s — 6)" + 3s = c)et 
FP hele OT? 


or from the pole-zero pattern with the use 
of (9). 

In such cases, these coefficients a: and by 
determine the approach to all asymptotes 
other than those at +90°. For those asymp- 
totes not equal to +90°, the loci in the upper 
half-plane converge from the upper side of 
the asymptote if b2 >a: and from the lower 
side of the asymptote if bs<a», 

A method of determining the final side to 
which the root locus converges to its asymp- 
tote has been presented. It is possible to 
make this determination analytically from 
the polynomial or visually from the pole-— 
zero pattern. 

C. S. LorENs — 

R. C. TitswortH 
Jet Propulsion Lab. 

Calif. Inst. Tech, 
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fessor of mathematics 

at Stanford University, Stanford, Calif., 
1948 to 1952; and visiting Professor of Engi- 
neering at the University of California, Los 
Angeles, in 1956. 

At the present Dr. Bellman is a mathe- 
matician at The RAND Corporation in San- 
ta Monica, Calif., and a consultant for 
the AEC, BOOZ, and Hughes Aircraft Co. 

Dr. Bellman has also been active in war- 
time and defense activities as a USAF in- 
structor in a preradar course (1942-1943), 
as an instructor in Specialized Training 
Program (1943-1944), as a mathematician 
in a sonar project at the U. S. Naval Radio 
and Sound Laboratory, Point Loma, San 
Diego (1944), as a mathematician on the 
atomic bomb project in the special engineer- 
ing division, U. S. Army, Los Alamos, 
N. M. (1944-1946), and asa member of Proj- 
ect Matterhorn, which was concerned with 
thermonuclear weapons, at Princeton Uni- 

versity, 1951-1952. 

He has published three books and over 
two hundred and fifty papers. 

Dr. Bellman is a member of the Council 
of the American Mathematical Society, a 
member of the Committee on Applied Math- 
ematics of the American Mathematical So- 
ciety, and was chairman of the Society for 
Industrial and Applied Mathematics, Los 
Angeles Chapter, 1956. 


R. BELLMAN 


Norbert T. Bold (S’55) was born in Chi- 
cago, Ill., on September 3, 1929. He received 
the B.E.E. degree from Marquette Uni- 

, versity, Milwaukee, 
Wis., in 1951 and the 
M.S. and Ph.D. de- 
grees from North- 
western University, 
Evanston, Ill, in 
1956 and 1958, re- 
spectively. 

From 1951 to 
1954, he served in the 
U. S. Navy as engi- 
neering officer aboard 
an aircraft carrier in 


sia 


, N. T. Botp 


the Korean Theater of Operations and a 
minesweeper in the Atlantic Fleet. 

While at Northwestern University, he 
was employed as a graduate teaching assist- 
ant in the electrical engineering department 
for three years. The last year of his educa- 
tion was completed with the assistance of a 
James M. Barker fellowship. 

Dr. Bold has been a student member of 
the AIEE. He is also a member of Tau Beta 
Pi, Eta Kappa Nu, Pi Mu Epsilon, and 
Sigma Xi. 


H. Chestnut (SM’57) was born on No- 
vember 25, 1917 in New York, N. Y. He 
received the B.S.E.E. and M.S.E.E. degrees 
in 1939 and 1940, re- 
spectively, from the 
Massachusetts Insti- 
tute of Technology, 
Cambridge. 

He has been with 
the General Electric 
Company since 1940, 
when he joined the 
Test Course. In 1942, 
he supervised the 
G.E. Electrical B Ad- 
vanced Engineering 
Program. Upon com- 
pletion of this program in 1943, he joined the 
Aeronautics and Ordnance Systems Division. 
From 1943 to 1956, he had assignments in 
AOS in Systems Development, Advanced 
Development, Guided Missiles, and Marine 
Systems sections. In 1954, he was appointed 
Ordnance Engineer, and in 1955, Project 
Engineer F-104, Since 1956, he has been 
Control Systems Engineer in the General 
Engineering Laboratory of General Electric, 
Schenectady, N. Y. 

Mr. Chestnut has been guest instructor 
in the Extension Division of Union College, 
Schenectady, a member of the Research and 
Development Board Committees, the chair- 
man of the AIEE Feedback Control Sys- 
tems Committee and the AIEE Subcom- 
mittee on Component Specifications, Vice- 
Chairman of the North American Control 
Council, and the first President of the Inter- 
national Federation of Automatic Control. 


H. CHESTNUT 


Nasli H. Choksy (A’50-M’55), for a 
photograph and biography, please see page 
117 of the December, 1958, issue of these 
‘TRANSACTIONS. 


Lester A. Gimpelson (S’55) was born on 
July 15, 1935, in New York, N. Y. In 1959, 
he received the S.B. and S.M. degrees in 
electrical engineering 
from the Massachu- 
setts Institute of 
Technology, Cam- 
bridge, where he is 
continuing his gradu- 
ate study. 

While on M.I.T.’s 
Cooperative Course 
program he was em- 
ployed by the Air 
Force Cambridge Re- 
search Center, Bed- 
ford, Mass., where he 
was engaged in the design of electronic and 
transistor circuits. He was also associated 
with the Research Laboratories of the 
British Thomson-Houston Company at 
Rugby, England. Since 1957 he has been a 
Staff member at M.I.T., first as a research 
assistant, and at present as a teaching as- 
sistant. 

Mr. Gimpelson is a member of Tau Beta 
Pi, Eta Kappa Nu, and Sigma Xi. 


L. A. GIMPELSON 


William A. Janos (M’59) was born on 
November 9, 1926, in Easton, Pa. He served 
in the U. S. Army from 1945 to 1947, and he 
graduated from Rut- 
gers University in 
1951 with a B.S. de- 
gree in physics. In 
1954 he obtained his 
M.A. degree in phys- 
ics, and in 1958 he 
received a Ph.D. in 
physics, both from 
the University of 
California, Berkeley. 
During this time he 
was the recipient of a 
University appointed 
teaching assistantship in the Physics De- 
partment and a Convair Scholarship Award. 

He has been employed by Convair and 
Convair Astronautics, San Diego, Calif., 
from 1951 until the present except for leaves 
of absence to obtain his advanced degrees. 
His work has covered a variety of fields in- 
cluding applied analysis and spectral theory 
related to analytical dynamics, wave propa- 
gation and diffraction, variational techniques 
in least time trajectories for thrust propelled 
flight, control system analysis and synthesis, 
noise theory and optimal linear estimation, 

He is presently conducting analytical re- 
search in digital-analog communication, 
computation and control systems, and initi- 
ating studies in information theory of con- 
tinuous sources. 

Dr. Janos is a member of the American 
Physical Society. 


W. A. JANOS 
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Robert Kalaba was born on September 
21, 1926 in Mt. Vernon, N. Y. He received 
both a B.A, in February, 1948 and a Ph.D. 
in February, 1958 
from New York Uni- 
versity, New York, 
IN@RY2 

He served in the 
U. S. Navy during 
1945-1946 and was 
a teaching fellow 
in mathematics at 
N. Y. U. from 1948 to 
1951. Since 1951 he 
has been a mathema- 
tician in the Elec- 
tronics Department 
of The RAND Corporation, Santa Monica, 
Calif. Over forty of his published papers 
have dealt with mathematical aspects of a 
variety of problems in physics, engineering, 
and operations research. Wave propagation 
and automatic control are among his current 
interests, 

He is a lecturer in engineering at Uni- 
versity of California at Los Angeles, and a 
member of Phi Beta Kappa, the American 
Mathematical Society, the Mathematical 
Association of America and the Society for 
Industrial and Applied Mathematics. 


R. KALABA 


Allen S. Lange, for a photograph and 
biography, please see page 71 of the May, 
1959, issue of these TRANSACTIONS. 
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Gordon J. Murphy (M’55) was born in 
Milwaukee, Wis., on February 16, 1927. He 
received the B.S. degree in electrical engi- 
neering from the Milwaukee School of En- 
gineering, Milwaukee, Wis., in 1949. After 
having done graduate work in an evening 
program at the University of Wisconsin, 
Milwaukee, from 1949-1952, he received the 
MS. degree in electrical engineering in 1952. 
He received the Ph.D. degree, also in elec- 


trical engineering, from the University of 
Minnesota, Minneapolis, in 1956. 

He was an assistant professor of electrical 
engineering at the 
Milwaukee School of 
Engineering until 
1951, when he ac- 
cepted a position as 
project engineer on 
inertial guidance sys- 
tems with the AC 
Spark Plug Division 
of General Motors in 
Milwaukee. In 1952, 
he accepted a posi- 
tion as instructor in 
electrical engineering 
at the University of Minnesota. There he de- 
veloped and taught courses in automatic 
control while continuing his graduate studies 
as a part-time student. In 1956 he was ap- 
pointed assistant professor of electrical en- 
gineering at the University of Minnesota, 
and he held this position until 1957, when he 
became associate professor of electrical en- 
gineering in the Technological Institute of 
Northwestern University, Evanston, Ill. He 
has been engaged since then in teaching and 
in research in the fields of statistical control 
theory, sampled-data theory, and adaptive 
control. 

Dr. Murphy is the author of two text- 
books on automatic control. 

He is a member of Sigma Xi, Eta Kappa 
Nu, AIEE, and ASEE. 


G. J. MurpHy 
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Joseph Otterman (S’52—M’56) was born 
in Warsaw, Poland, on April 12, 1925. He 
attended the Hebrew Institute of Tech- 

nology (Technion), 
=. -Haifa;. Israel, and 
graduated in 1947 
with the degree of 
Ingenieur in electrical 
engineering. 

From 1948 to 1950, 
he served in the Is- 
raeli Army, ending 
his military service as 
a lieutenant in the ca- 
pacity of Engineering 
Reconnaissance Offi- 
cer of the Command 
of the South. From 1950 to 1951, he worked 
in the Israeli Government Scientific Insti- 


J. OTTERMAN 


January ' 


tute on developing instrumentation systems | 
for ballistic measurements. 

He came to the University of Michigan, 
Ann Arbor in 1951, and obtained the M.S.E. 
and Ph.D. degrees in electrical engineering | 
from the University in in 1952 and 1955, 
respectively. 

From 1955 to 1959, he was associated 
with the Engineering Research Institute, 
University of Michigan, working on rocket 
investigation of the upper atmosphere within 
the International Geophysical Year pro- 
gram, on navigation problems, and on prob- 
lems of simulation by analog and digital 
techniques. In 1958 he proposed an experi- 
ment for determining upper atmosphere 
density and winds, utilizing a falling sphere 
and xyz accelerometers, and subsequently 
was engaged in investigation of this experi- 
ment under the sponsorship of the Air Force 
Cambridge Research Center, Bedford, Mass. 
He is currently with International Tele- 
phone and Telegraph Labs., Nutley, N. J., 
working at Bell Telephone Labs., Whippany, 
N. J., on topology and monitoring of com- 
munication systems. 

Dr. Otterman is a member of the Ameri- 
can Physical Society, Association for Com- 
puting Machinery, and Sigma Xi. 


Hanoch Ur was born in Bucharest, Ro- 
mania on October 16, 1931. He received the 
B.S. and Dipl.Ing. degrees from the Tech- 
nion-Israel Institute 
of Technology, Haifa, 
in 1954 and 1955, re- 
spectively. From 
1955 to 1956, he was 
a research assistant 
at the Computing 
Laboratory of Pur- 
due University, La- 
Fayette, Ind., and 
received the M.S.E.E. 
degree from Purdue 
in 1956. After a short 
period as a logical de- 
signer with Univac, St. Paul, Minn., he 
joined the Polytechnic Institute of Brook- 
lyn, Brooklyn, N. Y., as an instructor and_ 
graduate student. Recently he joined RCA 
in Camden, N. J., where he is working on 
high-speed computers. 

Mr. Ur is a member of Sigma Xi. 


H. Ur 


1960 


PGAC News 


Moscow CONFERENCE 


Preparations for the first Congress of the 
International Federation of Automatic Con- 
trol (IFAC) to be held in Moscow next June 
are being completed. In September, 1959, 
there was an IFAC General Assembly 
Meeting at Chicago, and Professor A. M. 
Letov, President-Elect of IFAC and Vice- 
Chairman of the U.S.S.R. National Com- 
mittee of Automatic Control, gave his re- 
port on the “Preparation for the First IFAC 
Congress.” With the approval of William E. 
Vannah, editor of Control Engineering and 
publicity chairman of the Congress, parts of 
Professor Letov’s report are printed here to 
illustrate the progress that is being made 
and the difficulties that are encountered in 
the preparation of an international technical 
meeting. 


I. INTRODUCTION 


1) The First IFAC Congress is to be held 
in Moscow, June 27 to July 8, 1960, in ac- 
cordance with the decision of the IFAC Gen- 
eral Assembly held in Paris in September, 
1957. 

2) The scientific program of the Con- 
gress was worked out by the U.S.S.R. Na- 
tional Committee of Automatic Control. 
This program was approved by the IFAC 
Executive Council in Zurich in March, 1958. 
An outline of the program was published in 
the IFAC Bulletin, as well as in a special 
invitation distributed by Dr. Ruppel among 
IFAC national member organizations. 

3) The U.S.S.R. National Committee 
continued the preparation for the Congress 
in conformity with this scientific program. 
The U.S.S.R. National Committee worked 
out a preliminary schedule of preparation 
for the Congress, of which the Executive 
Council was informed at the meeting in 
Rome in March, 1959. 

4) This report outlines the picture of the 
future Congress, as well as proposes solutions 
of some organizational problems. 


_IJ. SELECTION OF PAPERS 


In accordance with the decision of the 
IFAC Executive Council passed in Rome in 
March, 1959, the U.S.S.R. National Com- 
mittee, acting as the IFAC Papers Com- 
mittee, in the spring of 1959 began its work 


B of selecting papers for the Congress. 


The U.S.S.R. National Committee has 
been happy to see very high activity dis- 
played in Congress preparation by the ma- 
jority of the IFAC member organizations. 


Four hundred and nineteen abstracts of 


papers have been received from nineteen 


countries. 
The First IFAC Congress must demon- 


strate in the best way possible the scientific 


- 
: 


achievements of scientists and engineers that 
the Federation unites. Therefore the prob- 
lem of selecting papers is at present one of 
the most responsible and most important 


problems. The difficulties in solving this 


problem lie in the fact that the U.S.S.R. 
National Committee has received more than 


twice the number of abstracts of papers 
that can be presented and discussed at the 
Congress technical sessions. 

When considering these abstracts, the 
U.S.S.R. National Committee came to the 
conclusion that almost all of them touch 
upon very interesting problems of science 
and engineering, the elucidation of which 
would be most desirable. 

By June 27, 1959, approximately 220 
invitations had been sent to authors to sub- 
mit the full texts of their papers to the 
U.S.S.R. National Committee. After June 
27, 1959, the U.S.S.R. National Committee 
received 172 more abstracts. These abstracts 
were considered at the Committee’s meeting 
on September 1, when a decision to send 
additional invitations was passed. A total of 
321 invitations have now been sent. 

But the selection of papers is not com- 
pleted at this stage. The fact is that each 
abstract may be a good basis for an author 
to write a very interesting paper, and we ex- 
pect all the authors invited to make use of 
this opportunity. However, there may be 
some cases in which the full paper proves 
to be of insignificant scientific interest, 
which will give reason for its rejection. For 
example, Soviet papers will be very strictly 
selected. We expect to select 60 to 70 full 
papers. We also passed a decision that each 
author be allowed to present only one paper 
at the Congress, just as the U.S.A. Auto- 
matic Control Committee has decided. 

If we take into account that 10 to 15 
per cent of the authors cannot attend the 
Congress for some reason or other, we shall 
have between 210 and 215 papers to be pre- 
sented. These papers have been roughly 
classified into groups according to the prob- 
lems treated, which enabled us to draft a 
program for the Congress sections. 


Section I—Theory (7 groups of 141 
papers), 

Section I]—Components (6 groups of 
71 papers), 


Section I1I—Application (6 groups of 
118 papers). 


Of the eleven days (the duration of the 
Congress), two days will be devoted to 
plenary sessions, three full days to visiting 
plants, and six days to presentation of 
papers and discussions. At present the 
U.S.S.R. National Committee is working 
out a detailed program of the Congress. 


IIJ. NOMINATION OF SECTION CHAIRMEN 


After the final selection of papers to be 
presented has been completed in October, 
1959, it will be advisable to appoint chair- 
men and vice-chairmen of the three sections 
from among the Congress participants, in 
conformance with the division of papers into 
groups. 

Outstanding specialists in the fields cov- 
ering scientific problems which refer to these 
groups must be invited to serve as section 
chairmen and vice-chairmen. The U.S.S.R. 
National Committee shall nominate possible 
candidates to serve as chairmen and vice- 
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chairmen, work out a draft classification of 
scientific problems, and submit this infor- 
mation to the IFAC Executive Council for 
approval, 


IV. PREPRINTS, PUBLICATION, AND 
DISTRIBUTION 


The U.S.S.R. National Committee has 
organized seven scientific-technical com- 
mittees: 


Theory—Prof. J. Z. Tsypkin; 

Components—Prof. B. S. Sotskov; 

Automatic Electric Drives—Prof. M. G. 
Chilikin; 

Automation of Continuous Processes— 
E. P. Stefani; 

Automation of Technological Processes— 
Academician V. I. Dikuskin; 

Automation of Agricultural Production 
Processes—\. A. Budzko; 

Terminology—Professor M. A. Gavrilov 


to select abstracts of papers to be presented 
at the Congress, and they will be responsible 
for the preparation of the papers selected 
for publication; namely, translation of 
foreign papers into Russian, scientific edit- 
ing, etc. 


V. ORGANIZATION OF SIMULTANEOUS 
‘TRANSLATION 


During the plenary sessions and the sec- 
tion sessions, the U.S.S.R. National Com- 
mittee is arranging to have simultaneous 
translation made from English, French, and 
German into Russian, and from Russian 
into English only. It is very complicated to 
have simultaneous translation made from 
French and German into English. There- 
fore, the U.S.S.R. National Committee re- 
lies on the friendly assistance of Congress 
participants, who can speak various foreign 
languages, to organize simultaneous trans- 
lation from French and German into Eng- 
lish, and from Russian and English into 
French and German. It may happen that 
during the section sessions, the number of 
which exceeds three, some groups of the 
Congress participants will have no simul- 
taneous translation available. The U.S.S.R. 
National Committee asks all national mem- 
ber organizations to inform the Congress 
Secretariat (Moscow I-53, Kalanchevskaja 
ul. 15 A) of all participants wishing to assist 
in a friendly way in organizing simultaneous 
translation at the Congress sessions. 


VI. ConNGRESS COMPOSITION 


The U.S.S.R. National Committee and 
the “Intourist” Agency have received a 
number of letters whose authors request to 
be allowed to attend the First IFAC Con- 
gress. Taking this into account, the U.S.S.R. 
National Committee draws the attention 
of the IFAC General Assembly, Executive 
Council, and member organizations to the 
necessity of very carefully forming the dele- 
gations to the Congress. The work of the 
Congress must not be interfered with a great 
number of curious tourists. Therefore, we 
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urge that each delegation consist of Congress 
speakers, specialists engaged in research 
work or teaching, and others whom the re- 
spective national organizations may con- 
firm. 

The national member organizations 
should now begin to form their delegations 
to the Congress, and then send the U.S.S.R. 
National Committee lists of delegates in- 
dicating the name and surname, place of 
work, permanent address. This information 
will be required for the U.S.S.R. National 
Committee to provide the delegates with 
hotel accommodations and other facilities. 
The information should reach the U.S.S.R. 
National Committee not later than Febru- 
ary, 1960. 


Secretariat, Press Department, Excursion 
Office, Chief of the Congress Program 


The U.S.S.R. National Committee is or- 
ganizing the Congress Secretariat, the Press 
Department, the Excursion Office, and ap- 
pointing the chief of the Congress Program. 


Exhibition of Scientific and Technical Lit- 
erature 


In accordance with inquiries received 
from various people, institutions and firms, 
the U.S.S.R. National Committee agrees in 


principle with the proposal to organize an 
exhibition of scientific and technical litera- 
ure during the Congress. 


VII. ConGrEss REGISTRATION 


In order to reimburse partially the ex- 
penditure the U.S.S.R. National Committee 
is to bear, and in conformity with inter- 
national traditions, it is necessary to estab- 
lish Congress membership dues of $20 before 
March 31 and $30 afterward for people 
who are not presenting papers. Authors may 
register for $15 before that date, and for 
$20 later. 

Registration paid in advance of March 31 
will entitle registrants to a full set of pre- 
prints in English. Butterworths Scientific 
Publications, the Proceedings publisher, will 
distribute the preprints some time in April. 

A. M. LEtov 

President-Elec of IFAC 
Vice-Chairman of U.S.S.R. National 
Committee of Automatic Control 


Ir You ARE GOING TO 
IFAC’S ConGRESS... 


For information on arrangements for at- 
tendance at the First IFAC Congress, Mos- 
cow, U.S.S.R., June 27—July 8, 1960, write to: 


January 


Secretary, American Automatic Control 
Council 

(Editor, Control Engineering) 

330 West 42nd Street 

New York 36, N. Y. 


Inform him of your plans to attend the Con- 
gress by the first week in February. AACC 
needs your full name, title, place of work, 
and address so that it may assemble for 
IFAC a complete list of people who will 
attend from the U.S.A. 


1960 NATIONAL CONTROL 
CONFERENCE 


The 1959 National Control Conference 
was sponsored by the IRE with official par- 
ticipation by the AIEE, ISA, and ASME. 
The 1960 National Control Conference will 
be held in Boston, Mass., during September, © 
and it will be sponsored by the ASME with — 
participation of the IRE, AIEE, and the 
ISA. Harvey A. Miller, Manager, Electronic 
Development Division, Taylor Instrument 
Companies, 95 Ames Street, Rochester 1, 
N. Y., will accept papers for the IRE-PGAC 
members who wish to participate in the 
Conference. Further details about the Con- 
ference can be obtained from Mr. Miller. 


INSTITUTIONAL LISTINGS 


The IRE Professional Group on Automatic Control is grateful for the 
assistance given by the firms listed below, and invites application for In- 
stitutional Listings from other firms interested in the field of Automatic 
Control. 


THE RAMO-WOOLDRIDGE CORPORATION 
P.O. Box 45215, Airport Station, Los Angeles 45, Calif. 


cee 


PHILCO CORP., GOVERNMENT & INDUSTRIAL DIV., 4700 Wissahickon Ave., Philadelphia 44, Pa. 
Transac S-2000 All Transistor, Large-Scale Data-Processing Systems; Transac Computers 


The charge for an Institutional Listing is $75.00 per issue or $125.00 for two 
consecutive issues. Applications for Institutional Listings and checks (made 
out to the Institute of Radio Engineers, Inc.) should be sent to Mr. L. G. 
Cumming, Technical Secretary, Institute of Radio Engineers, Inc., 1 East 
79th Street, New York 21, N. Y. 


